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Dear  Lou: 

The  primary  emphasis  of  this  contract  was  to  calcu¬ 
late  the  parameters  appropriate  to  a  quasi-linear 
description  of  the  oceanic  internal  wave  field.  The 
main  results  are  described  in  Dewitt's  thesis (1).  There 
have  been  two  published  papers  (2,3)  on  the  early  part 
of  this  work.  A  paper  describing  the  later  calcula¬ 
tions  (4)  is  enclosed  and  will  be  submitted  for  publica¬ 
tion  shortly.  We  also  began  work  on  an  eikonal  descrip¬ 
tion  of  high  vertical  wave  number  internal  waves.  That 
work  is  not  yet  fully  finished.  In  addition  we  have 
been  trying  to  understand  general  features  of  nonlinear 
systems  with  the  hope  that  we  can  have  a  better 
phenowenolgical  picture  of  weakly  turbulent  systems  such 
as  the  oceanic  internal  wave  field.  That  resulted  in 
three  published  papers (5,8,7)  plus  another  paper  that 
has  been  submitted  for  publication (8) . 

Our  calculations  on  the  Internal  wave  field  provide 
better  numerical  values  for  lifetimes  and  other  parame¬ 
ters  than  those  existing  calculations  that  relied  on 
weak  interaction  theory.  However  the  effect  of  the 
large  scale  internal  waves  on  the  smaller  scale  waves  in 
the  induced  diffusion  region  (advection  and  doppler 
shifts)  are  so  strong  that  our  results  are  not  complete¬ 
ly  reliable.  In  regions  of  wave  vector  space  where 
these  effects  are  less  important  but  still  strong  enough 
to  invalidate  weak  interaction  theory  we  have  results 
that  we  believe  to  be  reliable.  In  particular  in  the 
region  of  vertical  wave  length  from  100  to  500  meters 
and  frequencies  larger  than  four  times  the  inertial  fre¬ 
quency  our  decay  rates  are  two  to  five  times  greater 
than  weak  interaction  rates  in  spite  of  the  claim  that 
weak  interaction  theory  is  valid  for  these  waves.  The 
reason  for  the  failure  of  weak  interaction  theory  is 
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that  all  members  of  a  resonant  triad  must  be  long  lived, 
not  just  the  member  whose  lifetime  is  being  calculated. 
In  most  other  regions  our  decay  rates  are  smaller  than 
or  equal  to  weak  interaction  rates.  Our  eikonal  ap¬ 
proach  is  an  attempt  to  handle  the  region  where  our 
present  results  are  not  completely  reliable. 
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SELF-CONSISTENT  EFFECTIVE  MEDIUM  PARAMETERS 
FOR  NONLINEAR  RANDOM  OCEANIC  INTERNAL  WAVES  ' 

Robert  Jon  DeWitt,  Ph.D. 

Department  of  Physics 

University  of  Illinois  at  Urbana-Champaign,  1982 

In  this  thesis  we  study  methods  for  obtaining 
self-consistent  parameters  that  describe  the  lifetimes  and 
frequencies  of  random  internal  waves  in  the  ocean.  We  discuss 
previous  treatments  of  the  problem  and  indicate  that  these 
methods  are  inadequate  in  large  part  because  they  are  not 
self-consistent.  We  perform  a  self-consistent  calculation  of 
oceanic  parameters  and  compare  our  results  to  previous  numbers. 
We  end  by  discussing  possible  future  improved  self-consistent 


calculations. 
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Self-consistent  effective-medium  theory  of  random 
internal  waves 

By  R.  J.  DeWITT  and  JON  WRIGHT 

Department  of  Physios,  University  of  Illinois  at  Urbana-Champaign, 

Urbana,  Illinois  61801 

(Received  6  June  1080  and  in  revised  form  11  May  1081) 

In  the  first  part  of  this  paper  we  introduce  a  path-integral  formalism  for  the  internal- 
wave  field  of  the  ocean.  The  intent  is  to  show  that  this  type  of  formalism  may  be  useful 
in  suggesting  improvements  to  current  calculations,  as  it  provides  a  framework  for 
applying  a  wide  variety  of  approximations  that  have  been  and  are  currently  being 
developed  in  other  areas  of  physics.  We  demonstrate  the  method  by  deriving  equations 
for  a  self-consistent  field  approach  (also  known  as  the  direct-interaction  approxi¬ 
mation).  The  experience  in  other  areas  of  physics  is  that  the  self- consistent  field 
approximation  is  more  reliable  than  lowest -order  perturbation  theory.  The  mid  result 
of  the  DIA  is  the  determination  of  an  effective  linear  model  for  the  description  of 
internal  waves  in  the  deep-ocean  environment.  In  the  second  part  of  the  paper  we 
obtain  H&aselmann’s  source  function  by  a  prescribed  limiting  process  and  are  able  to 
indicate  possible  improvements  in  related  calculations  by  comparing  the  limiting 
aasumptiens  with  numerically  oomputed  values. 


1.  Introduction 

In  the  past  two  decades  the  study  of  the  transfer  of  energy  in  the  internal-wave  field 
via  three- wave  (nonlinear)  interactions  has  become  an  area  of  active  research.  The 
possibility  of  three-wave  retonant  interactions  was  first  mentioned  by  Phillips  (1977). 
Subsequently  a  doltsmann-like  equation  giving  energy  transfer  rates  in  the  Resonant- 
Interaction  Approximation  (RIA)  was  obtained  by  Hassolmann  (1906,  1967)  using 
multiple-tiine-scaJe  methods.  Inherent  in  bis  calculation  was  the  assumption  that 
transfer  times  were  long  compared  with  typical  oscillation  periods,  a  reasonable 
assumption  given  the  widespread  belief  that  intemal-wave  interactions  were  weak. 
Later  synthesis  of  an  analytic  form  for  the  internal-wave  spectrum  by  Garrett  ft 
Hunk  (1979,  1975)  raised  hopes  that  results  of  numerical  calculations  using  Hsssel- 
maim’s  transfer  equation  might  contain  some  measure  at  reality.  Unfortunately, 
subsequent  calculations  fay  Others  (1976),  MoComae  (1977)  and  Pomphrey,  Heim  ft 
Watson  (1960)  predieted  vary  short  interaction  Hum  in  a  rather  large  region  of  the 
spsetn>m,antlmt  the  rewrite  were  iatsrnally  inoonsist  snt.  What  was  needed,  then,  was 
a  systematic  way  to  relax  the  strict  assumptions  inherent  in  the  dsrivation  of 
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most  of  the  emphasis  has  been  cm  critical  dynamics  and  folly  developed  turbulence. 
We  propose  to  formulate  the  dynamics  of  both  internal  waves  and  surface  waves  using 
these  methods. 

We  should  emphasize  at  the  outset  that  we  make  no  claim  that  the  formalism 
developed  herein  provides  a  complete  picture  of  the  very  complicated  dynamics  of 
internal  waves.  Rather,  our  interest  is  in  demonstrating  that  the  formalism  provides 
a  much  more  general  framework  in  which  to  study  internal  waves  than  do  the  earlier 
schemes.  In  fact  the  reason  for  mentioning  surface  waves  at  all  is  to  demonstrate  that 
the  formalism  is  easily  generalized  from  that  obtained  for  the  simple  three-wave 
interacting  system  used  to  model  internal  waves.  It  is  in  part  this  flexibility  that  makes 
this  method  powerful. 

The  aim  of  our  program  will  be  to  calculate  the  linear  response  functions  (Green’s 
functions)  and  the  two-time  correlation  functions.  The  ocean  will  be  assumed  to  be  a 
stationary  random  system.  The  time-development  of  perturbations  to  the  wave  held 
is  described  by  the  linear  response  functions. 

In  the  theory  of  turbulence  there  is  a  1  ng  history  of  attempts  at  setting  up  such  a 
formalism.  The  early  work  of  Kraichnan  (1858)  and  Wyld  (1881)  was  followed  by  a 
1873  paper  by  Martin,  Siggia  &  Rose  (1873)  (henceforth  MSR)  which  introduced  the 
idea  of  doubling  the  number  of  variables  in  the  problem,  leading  to  the  current 
formulation.  There  have  been  a  number  of  papers  which  have  rewritten  the  MSR 
formalism  in  various  ways,  see  for  example  Phythian  (1877)  and  Langouche,  Roekaerts 
&  Tirapegui  (1878).  We  choose  to  use  the  path-integral  formalism  of  Phythian  to  set 
up  the  problem.  The  reason  for  this  is  that  there  is  currently  a  great  deal  of  activity  in 
various  areas  of  physics  using  path  integrals.  This  will  enable  us  easily  to  make  use  of 
some  of  their  techniques.  In  this  paper  we  will  discuss  perturbative  and  self-consistent 
evaluations  of  the  path  integral.  In  fact  for  the  approximation  we  propose,  we  could 
use  directly  the  original  work  of  Kraichnan  (1868)  who  proposed  the  direct-interaction 
approximation  which  has  been  extensively  studied.  We  refer  the  reader  to  the  reviews 
of  leeKe  (1873)  and  Oreag  (1877).  For  a  similar  approximation  see  Phythian  (1868), 
Herring  (1885)  and  Edwards  (1864).  It  is  essentially  a  self -consistent  field  approxi¬ 
mation,  and,  in  the  language  of  diagrams,  vertex  corrections  are  ignored.  The  semi- 
classical  or  WKB  approximation  is  another  standard  approximation  to  path  integrals. 
We  will  not  discuss  that  here,  but  it  does  appear  that  it  is  the  appropriate  approxi¬ 
mation  for  part  of  the  interral-wave  Add. 

In  order  to  show  how  such  a  formalism  might  shed  more  light  on  previous  calcu¬ 
lations  we  show  that  a  particular  limiting  case  of  our  model  precisely  corresponds  to 
the  results  of  Hsssclmann  and  consequently  indicates  how  one  might  improve  an  the 
calculations  of  Others  and  MoComas.  We  also  show  that  the  formalism  splits  Hassel- 
mann’i  source  function  into  two  pteoes,  which  have  natural  interpretations  as  effective 
damping  and  driving  foress.  This  finally  leads  aa  to  interpret  our  model  of  the  internal- 
wave  field  ss  a  collection  of  uneoapied,  damped  oectilatoi*  in  the  presence  of  random 
driving  fbrcM,  in  the  spirit  of  effective-medium  theorise  found  in  other  areas  of  physios. 

In  j  2  w»  review  the  ierivztion  of  the  path  integral  for  stationary  random  pro  oames. 
In  ft  we  cfisoum  the  oorrelstion  functions  and  the  linear  response  functions  and  the 
direst  interaction  spprnilmsrinn  In  |4  we  interpret  the  not  Maser  problem  as  an 
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the  important  physios  that  is  involved.  The  final  section  is  devoted  to  a  summary 
and  discussion  of  some  unresolved  problems. 


2.  Path  integral 

In  this  section  we  will  sketch  the  derivation  of  the  path-integral  form  of  the 
characteristic  functional  describing  the  random  field.  For  some  of  the  more  technical 
mathematical  points  we  refer  to  the  literature,  Phythian  (197?)  and  Langouche  et  at. 
(1979). 

As  far  as  is  possible  we  will  use  the  notation  of  01  ben  ( 1 976)  for  the  mode  description. 
The  equations  to  be  solved  are  of  the  form 

A*(k,  t) + «Q*(k)  A*(k,  I) + v«(k)  A*(k,  t) 

+ £  JVfc,  *k,  *(k  -  k,  -  k.)  BzftX  ASK *) 

-/*(k ,t),  s-±l.  (£.1) 

Here  we  have  introduced  the  terms  v9(k)  A‘(k,  t)  and/*(k,  t)  as  external  damping  and 
driving  forces  on  the  internal-wave  field.  Hie  motivation  for  the  particular  choice  of 
these  terms  is  twofold.  The  first  is  that  of  simplicity  -  we  wish  to  demonstrate  the 
application  of  the  formalism  to  this  problem  without  introducing  complexities  that 
might  tend  to  obscure  the  general  ideas.  The  second  motivation  can  only  be  seen 
a  posteriori.  When  we  subsequently  write  down  the  DyBon’s  equations  for  this  problem 
it  rill  be  seen  that  new  terms  appear  which  have  the  same  formal  structure  as  the  two 
‘  external  ’  terms ;  this  leads  us  to  interpret  the  new  terms  as  ‘  corrections  ’  to  the  initial 
terms,  i.e.  as  ‘effective’  damping  and  driving  terms.  For  most  modes / and  vt  are 
negligible  and  can  be  ignored.  They  are  just  a  convenient  way  of  parametrizing  the 
coupling  of  inertial  waves  to  the  rest  of  the  world. 

The  amplitudes,  mode  frequencies  and  couplings  satisfy  the  relations: 

A*(k,  t)  -  [A-(  -  k,  t)]*,  0*(k)  -  fl-(  -  k),  (2.2),  (2.3) 

Bm.  “  (2-4) 

The  index  » is  either  +*  or  —  «  with  0+  >  0.  The  label  a  is  an  additional  mode  index, 
which  for  internal  waves  is  equal  to  unity,  If  an  appears  without  the  label,  it  is  taken 
to  be  Il+.We  make  the  ad  hoc  assumption  that  the  driving  foroes/g(k,  f)  are  stationary 
and  random  with  a  Gaussian  distribution  function, 

</»<M/*(k',0>  -  *<k+kU-,  jyk,f-n.  (2.6) 

This  particular  assumption  allows  us  to  obtain  closed-form  solutions  for  the  path 
integrals  that  will  appear  later.  Other  models  for  the  interaction  of  the  external  forces 
with  the  internal  waves  could  be  used,  such  as  random  parametric  couplings  in  (2.1), 
hut  such  models  are  somewhat  more  difficult  to  treat. 

The  probability  distributions  for  /  are  specified  fat  terms  of  its  characteristic 
ftmetfcnal  C(<f), 

m  -  (expjij d*  2  j_^  **-(  -  k,  l)/1k,  1 J) ,  (2.6a) 

C{f)  W-0).  <*«*) 
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The  equations  are  much  simpler  when  expressed  in  terms  of  their  Fourier  transforms, 


A‘( k,<o)  -  ^J^^k,!)*’ 

^*(k,  (O)  «  jt«p(k,t)dt, 

4(k,«)  =  f*  e<*It(k ,t)dt. 

J  -00 


With  this  convention,  equation  (2.6)  becomes 

C(4>)  =  expj-^sJ<PkJ*^"*(-k,  -w)J,8(k,w)^(k,w)|.  (2.8) 

We  are  interested  in  finding  the  probability  distribution  for  the  amplitudes  A.  To 
that  end  we  introduce  the  characteristic  functional  for  those  variables, 

m  -  <exp  »J  Jd*k  j  it  £-»{  -  k,  f)  .d*(k,  l)> 
or  in  terms  of  the  Fourier-transformed  amplitudes 

Z(|)  -  <oxp[2wi2  Jd*kJ<fo»£-*(-k,  -w)4*(k,w)]>.  (2.9) 

The  angle  brackets  refer  to  averaging  over  the  distribution  for  /.  They  are  given 
formally  as 

<X<J)>-j9fW(f)X<J)  (2.10) 

where  W(J)  is  the  probability  distribution  function  for  /  and  is  just  the  Fourier 
transform  of  C.  The  integral  is  a  functional  integral  in  that  we  are  to  integrate  over  all 
possible  functions /(k,  f)  for  each  mode  k.  One  way  to  achieve  this  is  to  discretize  k 
mull  and  integrate,  ^,nn#|g.  (2.11) 

*  a 

For  the  purposes  of  this  paper  it  is  not  necessary  to  be  specific  and  Pbythian  (1977)  can 
be  oraumlted  for  more  details.  The  idea  now  h  to  replace  1T(/)  by  its  Fourier  transform 
Ctf)  and  to  change  variables  from /to  A  wing  (2.1): 

<*>  (2.12) 
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1977;  Langouche  et  al.  1979)  and  is  given  formally  by 

J£  -  exp  j  +  (2.14) 

It  has  been  shown  (Langouche  el  al.  1979;  Dominicis  &  Peliti  1978)  that,  if  diagrammatic 
perturbation  theory  is  used,  the  Jacobian  cancels  ‘tadpole’  diagrams  and  both  the 
tadpoles  and  the  Jacobian  can  be  omitted.  See  the  appendix  for  a  discussion  of  this 
point. 

We  now  extend  the  definition  of  the  angle  brackets  to  include  <j>  functionals  as  well 
as  A‘(k,  w)  functionals, 

<X(A, $))  -j®AS>4§£ X(A, <t>) 0(4) 

xexp{-2m^jd*kdto4~‘(-^>  —  to)f*(k,u)}.  (2.1 5) 

The  generating  functional  introduced  in  (2.9)  is  now  extended  to  include  <j>, 

Z(£,y)  -  <exp{2ff»2  fd*kfdto(i-(-k,  -to)A>(k,to)  +  V-(-k,  -«)0*(k,  «)]}>.  (2.16) 

As  we  are  presumably  interested  in  correlation  functions  of  A‘(k,  to)  which  can  be 
obtained  by  differentiation  with  respect  to  £,  it  seems  unnecessary  to  introduce  y. 
However,  the  field  4  (and  hence  y)  plays  a  crucial  role  in  defining  a  response  function 
and  in  simplifying  the  perturbation  series.  To  demonstrate  this  we  imagine  replacing 
/-*•/+  e  in  the  equation  of  motion.  This  shows  up  in  (2.15)  as  an  additional  factor  in 

the  integral  exp{+2ff»  £  Jd*kJdw0-*(-k,  -«)e*(k,w)}. 

If  we  expand  this  to  first  order  in  e  and  consider  the  change  in  (A)  we  see  that 

(u-tffi-h,)  -  *'«'*•<*.<■•)<**.»».  I*-*’) 

which  can  also  be  written 

(i£f= iff?))  “  »'<^'(k,«)^(k',<')>.  (2.18) 


Thus  the  term  on  the  right-hand  side  of  (2.18)  gives  the  linear  response  of  A*(k,t)  to 
a  known  perturbation.  It  is  linear  only  in  the  sense  that  e  is  treated  to  first  order.  Hie 
rest  of  the  dynamics  is  nonlinear.  Response  mid  correlation  functions  are  calculated 
by  taking  the  appropriate  functional  derivatives  of  Z,  for  instanoe 


<Atk,w)^(k',«')> 


1  /  1  \«  »  * 

f  te)  jfh i--.)^p=r 


*(£,?)  I  •  (*•!») 


The  standard  trick  used  to  develop  perturbation  theory  is  to  notice  that,  if  the  theory 
were  Hnsar,  all  of  the  integrals  in  (>.16)  sonld  he  obtained  exactly.  In  the  actual  case 
the  factor  preventing  exact  analytical  integration  is  given  by 


sxp{-2>rf£  /#k#kt#M(»-*i-»i)Ja(k-l‘i-k^ 

x*SS£A*<ki,ml)A*Qt»'et)4-+(-k,  -a.)}. 
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If  we  recognize  that  A’( k,  w)  can  be  replaced  by 


and  each  tj>  by 


^*(k1w)^2w.Jg-t(_k  _w) 


£(f,i/)  =  exp  |  ^ j  J *{«! + w8 + ws)  +  lc,  +  kg)  rfw*  dw8  <f*k3 

*  ,,E.  2n‘1'a‘  ^-(-kn  -«i)  ^-(-kg.-o*)  t£~‘'(  -  kg,  -oj  Z°(5,  ^ 

(2.20) 

The  term  Z0(£,  y)  is  obtained  from  (2. 15)  and  (2. 1 8)  by  finding  all  of  the  integrals  now 
that  the  nonlinear  piece  has  been  removed.  The  result  is 

*■«.»>  - 

X  exp  {2*i  Jd*k  S  £-( -k,  -a>) %(k, «)} .  (2.21) 


The  new  functions  D  are 

D*(k,  w)  =  K#(k) + »n*(k)  -  io),  D+(k,  w)  -  [2>-(  -  k,  -  «)]*.  (2.22) 

We  can  now  write  down  a  diagram  expansion  for  the  perturbation  series.  We  do 
that  in  the  next  section. 

In  principle  quantities  of  interest  can  be  calculated  directly  from  (2.15)  without 
resorting  to  the  perturbation  theory  that  follows  from  (2.20).  Currently  there  is  a 
major  effort  aimed  at  extracting  non -perturbative  information  directly  from  path 
integrals  describing  a  wide  variety  of  physical  phenomena.  As  those  methods  tend  to 
be  difficult  to  implement,  we  propose  to  first  test  the  SCFA  (DIA).  There  are  a  number 
of  quantities  of  interest  for  the  internal-wave  field.  However,  there  is  relatively  meagre 
information  on  higher  moments  so  in  this  paper  we  will  be  primarily  interested  in  the 
two-point  functions, 


1  V"f(k,w)fi(k+k’)H<o+<u')  -  <A*(k,wM*'(k>')>, 

*V(k,t)«»(k+k')  -  <A*(k,iM'(k,0)>, 

and 

k,w)a«(k+k')i(«+w')  -  »<A*(k, w)^'(k',w')>, 

fr'lk.W k+k')  -  «A*(k,l)^(k',0)>. 

The  definitions  of  V  and  0  lead  to  some  important  symmetry  properties : 

£7-f(k,«)-t£7-^'(-k,-w)]*  C7*(k,*>)  -  Cf'*(-k,  -«), 

0*(k,w)  -  [0— '(-k,  -«)]*. 


(2.281 

(2.24) 
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3.  Perturbation  theory 

In  tills  section  we  will  exhibit  the  general  structure  of  the  two -point  functions  and 
diseuss  briefly  the  diagrammatic  perturbation  theory.  We  begin  by  listing  the  free 
Green’s  functions  and  correlation  functions,  and  we  will  then  point  out  the  general 
features  that  the  exact  functions  will  have.  The  free  correlation  function  is  defined  by 

1 UU  k,  w)3*(k+k')i(o+«')V^ 

~  k  fe)  sppfe  -«')  z°  Lo 


2w  D~*{  -  k,  -u) 


The  free-responee  function  is  given  by 


)5=(£^) 

„5gE iA 

k)+t«-Q 


i*(k+k')«(w+w')d. 


±0y&U)Pik+V)9i»+m’)i^ 

-^nk,«)^(k',«')> 

id£-(‘-k,-«)Jrv<-k',  -^)  wL. 

If  if*  is  independent  of  »,  the  corresponding  time-dependent  functions  are 


I.I)aod(8.S), 
on  «.  To  the 


at  Mast  Hn  start-time  behaviour.  If  the 
%  sinmtelanl*  1 1  "  liiuw  tlnSi  tmbairlmir  wffl 

aft*  HK»  iwlP'lfllW  wv»l*“IWwr  W1U 
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i  s,  k,  to 

£H(k,  u) 

-F(  k)  j,  k,w— ► 

D-'  (-k,  -w)  0*(k,  w) 

Figure  1.  Lines  representing  response  function  and  correlation  function. 


Figure  2.  Lowest -order  vertex  function.  The  line  with  the  incoming  arrow  is  the  first  argument 

x«*(k,-k,  +  k,)x  i(w, -<j, +  (■>,). 

— -o— 

Figure  3.  Lowest- order  graph  contributing  to  the  correlation  function,  U. 

Because  of  its  symmetry  there  is  a  factor  of  l  associated  with  it. 

- 1 -  -  - > -  ♦ 


Figure  4.  Dyson  functions  for  the  two-point  functions. 

two  other  lines  that  either  have  arrows  pointing  away  from  the  vertex  or  no  arrows. 
In  figure  2  we  show  a  vertex  which  has  associated  with  it  the  factor 

+  ^(wi  -  +  W*). 

The  first  pair  corresponds  to  the  arrow  in  the  diagram .  B  is  symmetric  in  the  other  pairs. 

(c)  All  internal a,  k  and  «’s  are  Bummed  over. 

(d)  For  Ath-order  perturbation  theory,  draw  all  topologically  distinct  diagrams. 
Oraphs  possessing  a  symmetry  have  a  symmetry  factor;  for  example  the  graph  of 
figure  3  gets  a  factor  of  J. 

(e)  Each  diagram  gets  a  factor  (2w)-*w+w. 

(/)  Each  external  A  field  gets  a  factor  of  ( - »).  (The  extra  i  here  is  because  the 
response  function  is  see  (3.2).) 

We  mention  that,  if  vertices  with  more  lines  are  present,  the  only  rule  that  is  changed 
is  (e).  Each  vertex  will  still  have  only  one  line  with  an  arrow  pointing  into  it. 

The  Dyson  equations  have  been  given  by  MSR.  For  completeness  we  show  them  in 
figures  4—0.  Note  that,  although  there  is  only  one  kind  of  bare  vertex,  there  are 
three  kinds  of  renormalized  vertices  which  are  shown  to  lowest  order  in  figure  9. 
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PiotJRK  5.  Equation  for  Zt  and  Z)  in  terms  of  renormaUsed  vertices. 


IXD4  -  >*  * 


£»  -  £>*  — - 

Fiodbk  6.  First  terms  in  expansion  of  vertices. 


Finally,  it  can  easily  be  shown  that  any  graph  containing  a  dosed  loop  whose  lines  ate 
entirely  response  functions  (arrows)  is  identically  sera. 

'hie  Dyson  equations  ate  easily  solved  to  give 

«K(k,to)  -  (D-«(k,a.)»,i^-<Ef(k,W)]-‘  (3.5) 

(note  that  this  is  a  2  x  2  matrix  equation).  The  correlation  function  is  given  by 

IK  -  0"»(k,  o)  [FJk,  m)  + £f*(k,  a»)]  (?*'*»(  -  k,  -«)  (3.6) 

(implies  summation  over  «„  «*). 

If  vertex  corrections  are  ignored  (direct  interaction  approximation)  then  we  have 
£f(k,  w)  =  -  S  j  d*k1d*ktd6J1dw,3*(k+k1-k1)5(w+w1-w1) 

x  Big%  *&&&**'(** «»)]  t  -  U«*’(k»  «,)]  (3.7) 

and 

£f(k,«)-£  S  J  *1%  tPk,  ifcii  «fa»i  3*(k + k, + 1^ )  i{u + &>, + w, ) 

(3.8) 

ZS,  has  the  same  symmetry  properties  as  U  aud  ffe,  has  the  sane  as  6. 


U  and  <3  can  be  ignored.  U  and  O  are  assumed  to  have  the  form  U*  -  DM,-,.  The 
corresponding  statement  fori  is  £'  « 


•  r :•> •  ■  :  *r-v.- ■ 
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Figure  7.  Graphs  for  £  for  four-wave  interactions. 

We  now  comment  briefly  on  the  application  of  these  techniques  to  surface  wares. 

The  principal  difference  is  that  the  equations  involve  four-wave  coupling,  that  is 

there  is  a  term  „  .  .  , 

BxAxAxA 

in  the  equation  of  motion,  (1.1).  The  term  involving  only  two  A’b  is  still  there,  but  it 
is  not  possible  to  have  three-wave  resonanoes.  This  means  that  in  (3.12)  the  limit 
»>->0  gives  aero  and  the  answer  must  come  from  finite  v  which  presumably  gives  a 
small  contribution.  In  any  case  if  we  replace  the  BAA  terms  by  BAAA  in  (2.1)  we 
obtain  a  very  similar  theory .  The  graphs  for  £  are  shown  in  figures  7  and  they  are  to  be 
solved  self-oonsistently  with  the  diagrams  of  figure  4.  The  number  associated  with  a 
vertex  is  now  3!  instead  of  21.  There  are  some  additional  (2ir)  factors.  The  net  effect  is 
to  change  rule  (e),  (2!)*(2n)-i<*+*‘>-> (3!)* (2*)-**'+*, 

and  rule  (b)  gets  a  factor  of  t  for  each  vertex.  These  two  changes  can  be  stated  in  a  way 
that  is  easily  generalized.  If  there  are  NA  A  fields  at  a  vertex,  the  vertex  has  associated 
with  it  (i^+*)^!.  Each  dosed  loop  has  associated  with  it  an  integral  fdw/Zn. 


4.  Linear  effective  medium 

In  this  section  we  shall  elaborate  on  the  effective  medium  picture  and  discuss  the 
direct-interaction  approximation  and  the  relationship  to  the  theory  due  to  Haasel- 
mann.  We  will  first  remind  the  reader  of  a  few  simple  properties  of  linear  systems,  and 
compare  these  with  our  effective  medium.  The  simplest  situation  assumes  no  memory 
and  is  given  by  ^ + iQ>A, + yA,  m  ^  (4  , ) 

which  is  (2.1)  with  .8  =  0  and  the  mode  index  k  ignored.  If  we  define 

00.0  -  (A*(t)A(t')>,  (4.2) 

then  in  the  steady  Btate  U  depends  only  oat—t'  and  is  given  by 

C(t)  „  f"  —  —  (4  31 

w  2»r  (v-iQ+w)  (i>+iQ—w)‘  [  > 

If  we  calculate  dU(t,  t)/dt  we  obtain 

iC(t, t)/dt  -  - 2i»0(O)  +  J-0,  (4.4) 

where  J  is  given  by  (2.5)  and  (2.7).  This  just  represents  the  belanoe  between  the  driving 
forces  and  the  dissipative  forces.  We  now  compare  with  the  correlation  function 
obtained  in  the  previous  section, 

°'{t)  " J (Ji+5® L+VQ-^T-ft^)]-  (45) 
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We  hare  suppressed  the  mode  index  lc.  If  £1(  £t  are  approximately  independent  of 
10,  then  it  is  clear  that  real  is  a  frequency  shift  and  is  an  effective  visooeity 
representing  loss  to  the  other  modes.  This  particular  interpretation  may  be  somewhat 
misleading,  however.  Consider  the  case  where  the  mode  in  question  is  Doppler-shifted 
because  of  passage  through  an  inhomogeneous  medium  (large-scale  wave).  In  this 
case  we  would  formally  say  that  one  mode  has  been  damped  while  another  mode  of 
neighbouring  wavenumber  has  grown.  In  reality  of  course  we  are  just  seeing  the 
original  wave  with  a  slightly  different  frequency. 

In  addition  to  the  ‘effective  damping ’  term  there  is  also  an  effective  driving  foroe 
due  to  the  other  modes  and  its  correlation  function  is  The  correlation  function  for 
the  total  driving  foroe  is  thus  2^+.F#. 

Now  in  general  2^  and  2^  will  depend  on  u.  In  this  case  it  is  easy  to  see  that  the 
equation  giving  the  evolution  of  the  amplitudes  is  given  by 


i(t)+r  t(t-r)A(t’)dt'-f(t),  £(t)  .  o,  <<o. 

(4.6) 

J  “• 

If  we  write 

r+(w)  -  J"  e«  P(i)  dt,  r-(«)  -  e“  t*(t), 

(4.7) 

then  we  see  that  F^w)  is  analytic  in  w  for  fto  >  0.  The  corresponding  response 

functions  are 

0(),“  *irj_„r^(w)-»a.’ 

(4.8) 

which  leads  to  an  obvious  interpretation  ofZf(eo), 

T+(«)  -  v,+*0-£J(a»). 

(4.8) 

The  fact  that  Q(t)  m  Ofor  t  <  0  (causality)  implies  that  there  are  no  sevos  of 
in  the  upper  half-plane.  Because  of  these  analytical  properties  of  O  and  I\  Kramer’s 
Krtinig  dispersion  relations  can  be  used  to  calculate  the  real  part  of  O  and  £  ones  the 
imaginary  part  is  known  (or  vice  versa). 

The  equilibrium  condition  relating  the  energy,  dissipation  and  driving  forces  is 
obtained  by  evaluating 


at  t  m  0. 


F(u) _ do> _ 

2tt  [p-(  —  hj)  —  uo]  (T+(<ti)  —  uo] 


5.  Relationship  to  previous  results 

In  this  section  we  obtain  the  results  of  Hasselmann  as  a  special  case  of  our  model 
equation.  The  limiting  assumptions  made  in  obtaining  the  results  provide  a  more 
transparent  way  to  check  the  validity  of  the  assumptions  inherent  in  Hasselmann’s 
equation;  in  the  end  of  this  section  we  indicate  where  improvements  could  be  made  in 
his  scheme. 

To  obtain  the  results  of  Hasselmann  one  performs  the  following  operations  on 
(9.7)  and  (3.8)  (using  only  £•  *  £“**). 


Vi. 
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(1)  Assume  that  Ft,  Z„  and  I,  are  independent  of  to  and  perform  the  Wj  and  w, 
integrations.  For  example,  defining  »>(k)  m  vt(k) + >S1(k),  we  get 

Jdu,  da,  *(«  - -  «,)  M^i)* + (Q(k,)  -  a^)*)-1  [Kk,) + *Q(k,)  -  *wJ-» 

=  j^-)[a»-0(kl)-n(k#)-t(^(kl)  +  Kk.))] 

x[(<u-Q(k1)-Q(k,))*+(>»(k1)  +  Kki))*rl.  (8.1) 

(2)  Take  the  limit  Ft,  vt,  Slt  £,->0  (in  such  a  way  that  6 (k,  0) «  F/2v  is  main¬ 
tained).  Then  (5.1)  becomes 

<«•*) 


where  &  represents  a  principal  part  integral.  When  the  actual  structure  of  the  B 
coefficients  is  considered  it  is  easily  seen  that  the  delta  function  contributes  to  v(k) 
and  the  principal  part  to  a  frequency  shift.  The  Dirac  delta  function  has  come  from 
the  limit 


v 

p*-t-(A£l)* 


nt(AQ) 


(0.3) 


and  we  call  the  function  on  the  left-hand  side  the  frequency  filter  for  the  problem. 
Such  a  function  was  earlier  obtained  by  Holloway  A  Henderson  (1977)  and  Holloway 
(1978, 1979),  where  he  also  noted  the  RIA  as  a  limiting  case.  This  is  the  first  indication 
of  a  possible  improvement  in  the  calculations  of  Olbers  and  McComas.  What  they  did, 
in  essence,  was  to  assume  that  v  was  small  compared  with  typical  oscillation  fre¬ 
quencies  and  then  calculated  values  of  v  on  the  order  of  or  even  greater  than  the 
frequencies.  The  earlier  assumption  of  neglecting  the  diagonal  elements  of  U,  G 
corresponds  approximately  to  assuming  v/Q  <  1. 

Equations  (3.7)  and  (3.8)  must  be  solved  self-consistently,  however,  so  that  values 
of  v(k)  used  for  input  on  the  right-hand  side  are  the  same  as  those  obtained  on  the 
left-hand  side.  Unfortunately  solving  these  equations  self-consistently  does  involve 
integration  over  one  more  variable,  mid  so  may  require  an  inordinate  amount  of 
computer  time  (the  two-dimensional  integral  is  already  relatively  time-consuming). 

Having  now  carried  out  the  limiting  process,  we  are  ready  to  make  a  detailed  com¬ 
parison  with  the  source  function  of  Hasselmann,  which  he  writes  as 

SF  -  /d%ld*k,{2’+<^k-k1-k1)d(n-fl1-Qt)[iH»H-nih+»»4] 

+  2 T~  i*(k -k,+kt)J(Q-fl1  +  Q1)[*1*1+n*1- »**],  (6.4) 

where 

»(k)*(k+k')  -  <A+(k,f)A-(kM)>5^-). 

We  split  i?  into  two  pieces. 


n k)  -  fd%d*kt(Tm-kl~k,)HQ^at-Ot)nlnt 

+ST~S(k-k1+kt)t(Q-Ql+(it)»1*t}  (5.8) 

and 
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Figure  8.  Plots  of  o(k)/fl  «t.  m  for  various  frequencies.  The  values  of  fl  for  each  curve  are  as  fol¬ 
lows:  (A)  7-0  x  i(h»  s-‘,  (B)  7-8  x  10-»  s-»,  (0)  8-0  x  lO"4  s-»,  (JO)  11  x  10-*  *-»,  (*)  1-7  x  10~»  s~>, 
(F)  3-2  x  10-*  s-‘,  (0)  MxKHr1,  (H)  1-8  X  10-*s~l,  (J)  2-8  x  10"*  s-1.  This  labelling  applies  to 
all  curves  in  figures  8  and  8.  The  values  of  Q  correspond  to  equal  increment*  in  log  [(Jt{ +!£)/*{]. 


We  find  the  following  identifications, 


and 


H(k)^,(k)  (5(M) 


(5.6) 


^  ^(k)  -  2(Kk)  -  y,(k))  -  ^(k),  (5.7) 

where  0  ia  defined  in  (2.28),  S,  is  the  effective  force  due  to  all  other  modes,  and  F  is 
the  total  force,  F(k)  -  J,(k) + J^(k).  (5.8) 

From  equation  (4.4)  we  have,  in  the  stationary  state, 

F  -  +  2v0(k,O)  -  F^k)+L1( k) 


or 


7(k,0) 

m  *..«.>  -^(k)  s\(t)  ^k) 


(*.») 


«  mi 
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From  this  equation  we  see  that  the  calculations  of  Olbers  and  MoComae  give  the 
coupling  to  external  systems  in  terms  of  the  difference  between  an  effective  force  and 
an  effective  viscosity. 

Earlier  we  noted  the  breakdown  of  the  resonant-interaction  approximation  (RIA) 
because  of  the  large  values  of  v(k)  in  some  regions  of  the  internal-wave  spectrum.  We 
would  next  like  to  test  the  assumptions  that  E*  and  S,  are  independent  of  t».  We 
consider  first  the  case  of  Ev  In  order  to  ascertain  whether  the  non-resonant  structure 
of  £,(<u)  might  have  any  appreciable  effect  on  the  values  computed  in  the  right-hand 
side  of  (3.7)  we  choose  a  particularly  ample  form  of  2^: 

-  a(k)  +  («-Q(k))d<k),  (All) 

where  both  a0t)  and  b{ k)  are  assumed  real.  If  the  »  constant  assumption  is 
warranted  we  should  expect  to  compute  very  small  value*  for  b(k).  For  the  present 
we  still  assume  that  E,  -  ocmstant  and  use  the  RIA.  The  effect  offr(k)  is  to  reduce  the 
residue  of  the  pole  of  6  from  l  to  1/(1  +ib)  and  to  decrease  the  damping  from  v  -  a 
to  v  -  o/(l+ft*).  The  correlation  function  is  reduced  by  the  same  amount: 

We  wish  to  check  the  siw  of  b  numerically.  We  take  0(k,  0)  to  be  given  by  GM76,  a 
recently  revised  version  of  the  Garrett-Murk  spectrum  (see  MoOomas  4  Bretherton 
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1977).  The  coopting  coefficient*  used  in  the  calculation  am  derived  in  the  we  manner 
as  those  given  by  MoComas. 

We  have  computed  values  of  o(k)  and  b(k)  by  assuming  that  both  wen  initially 
zero  (i.e.  just  £,->0  as  before),  and  then  iterating  (8.7)  cnee.  Values  of  o(k)/Q(k) 
and  6*(k)  an  shown  in  figures  8  and  9  respectively.  We  first  note  the  previously 
mentioned  fact  that  the  values  of  o(k)  for  much  of  the  spectrum  violate 


o(k)/(0(k)/2»)  <  1. 


The  dashed  line  corresponds  to  o(k)/(0(k)/9ar)  —  1.  Ike  RIA  assumes  that  Q  ►  a; 
that  condition  is  clearly  violated  hare .  We  further  note  that  an  examination  of  figure  9 
reveals  values  of  6*(k)  that  are  extremely  large  for  much  of  the  spectrum.  This  dearly 
violates  the  assumption  that  Z1  is  constant  in  *».  We  oondude  that  any  reliable  self- 
consistent  calculation  must  indude  information  on  the  non -resonant  structure  of  Z^ 
It  should  be  noted,  however,  that  the  simple  linear  form  of  Zi,  chosen  here  to  test  the 
constant — Zj  assumption,  is  wholly  inappropriate  for  a  serious  calculation. 

Finally,  we  can  also  determine  the  validity  of  the  assumption  that  Z,  is  constant 
in  u.  When  we  actually  compute  Z^k,**)  we  find  that  Z^  actually  does  possess  an 
appreciable  structure  with  a  finite  width  in  a*.  However  it  is  not  sufficient  merely  to 
calculate  Z,  for  various  values  of  a,  since  frequency  filter  cuts-off  may  allow  only 
values  of  Z,  from  a  narrow  band  around  resonance  to  enter  the  calculation.  In  that 
case  Z,  may  be  considered  effectively  oonstant.  One  way  to  circumvent  this  is  to 
compare  the  functions 

and 

B(t)  - 


We  note  that  Z,  always  appears  in  the  equations  through  U{u),  which  contains 
information  on  the  width  of  the  frequency  filter  as  well  as  on  the  width  of  Z,.  On  the 
other  hand  JR(<)  contains  information  on  the  width  of  Z,  only.  We  envision  two  possible 
extreme  cases.  In  one  oaae  v  is  extremely  small,  so  that  the  frequency  filter  becomes  a 
resonance  delta  function  and  V(t)  assumes  the  form  earlier  obtained  by  assuming 
Z,  -  constant.  On  the  other  extreme,  v  might  be  so  large  that  the  width  of  the  fre¬ 
quency  filter  would  be  much  larger  than  the  width  of  Z,.  In  this  caee  we  would  expect 
the  shapes  of  V (!)  and  £(*)  to  be  identical.  The  calculations  that  we  have  performed 
indicate  that  both  extreme  situation  do  indeed  obtain  hr  different  regions  of  the 
spectrum.  Figure  10  shows  the  first  ease.  Here  Z,  is  oonstant  and  0(f) ace***-*1*1. 
If  Z,  were  actually  constant,  however,  R(t)  should  be  a  delta  function  in  t.  In  actuality 
we  see  that  B(i)  does  have  some  width,  indicating  that  Z,  is  not  constant  but  that  the 
variation  is  unimportant  in  the  calculation.  On  the  other  hand  figure  li  shows  the 
oese  where,  appropriately  sealed,  U{f)  end  )  are  identical,  indicating  large  values of 
r(k).  In  this  region,  then,  not  only  is  the  RIA  very  bad,  but  aho  the  assumption  that 
Z,  -  constant  must  be  dropped  if  the  frequency  filter  cOrreoticn  is  to  be  used. 

We  have  also  calculated  Z,(w)  and  from  it  the  response  fraction  to  verify  directly 
that  the  frequency  dependence  of  Zj  iehaportaat  for  some  modse. 

In  general,  than,  we  eohdode  the!  any  self-consistent  calculation  of  Z*  and  Z,  in 
the  region  whew  When  and  HeOeaea  found  inconsistencies  in  their  calculations  must 

necessarily  inriode(l)aftequanoy  filer fruMtioai  to  relax  the  RIA,  end  (t)non^aeonant 
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information  on  and  Unfortunately  this  greatly  increase*  the  difficulty  of  the 
numerical  solutions  of  the  problem,  since  it  requires  an  extra  integration  over  at, 
including  a  principal  part  integral  that  has  not  been  present  in  earlier  calculations. 
In  addition,  the  equations  must  be  solved  self-consistently,  implying  an  iteration 
process  on  an  already  time-consuming  calculation.  We  are  currently  attempting  a 
simplified  version  of  the  self-consistent  calculation. 


6.  WKB  approximation 

In  the  introduction  we  alluded  to  other  approximations  to  the  path  integral.  One 
of  the  failings  of  weak  interaction  theory  can  be  traced  to  the  interaction  of  two  small- 
scale  high-frequency  waves  with  a  large-scale  inertial  frequency  wave.  This  interaction 
has  been  called  induced  diffusion  by  McComas  ft  Bretherton  (1977).  It  is  appropriate 
to  think  of  it  as  the  effect  of  advectkm  of  the  small-scale  wave  by  the  inertial  wave. 
This  picture  suggests  studying  the  propagation  of  a  small-scale  wave  packet  in  the 
geometric  optics  approximation  including  the  effects  of  a  time-varying  background. 
The  aim  would  be  to  calculate  a  coherent  decay  time.  The  time  scales  of  the  response 
functions  in  the  induced  diffusion  region  are  dominated  by  processes  that  leave  the 
wave  packet  coherent,  although  somewhat  altered. 


7.  Summary 

We  have  written  this  paper  with  a  twofold  purpose.  In  the  first  part  of  the  paper  we 
have  introduced  a  formalism  for  solving  nonlinear  equations  of  the  type  encountered 
with  internal  waves.  From  experience  with  self-consistent  field  treatments  in  other 
areas  of  physics  we  very  well  might  expect  such  methods  to  be  helpful  in  extending 
calculations  beyond  some  of  their  inconsistencies.  To  this  end  we  have  taken  a 
simplified  steady-state  model  of  the  oceanic  internal  wave  field  and  obtained  via  path 
integral  methods  the  corresponding  diagrammatic  perturbation  theory  from  which 
Dyson’s  equations  can  be  obtained .  In  the  process  we  have  suggested  as  interpretation 
of  the  ocean  as  a  collection  of  uncoupled,  damped  oscillators  driven  by  random  forces 
and  each  having  a  memory.  In  this  case  we  used  the  direct  interaction  approximation 
of  Kraichnan,  which  ’Is  easily  implemented  in  the  diagrammatic  language.  We  should 
note  that  one  problem  with  this  approach  is  that  it  is  not  currently  known  what  kind 
of  approximations  guarantee  the  positivity  of  quantities  like  the  effective  viscosity 
and  also  guarantee  consistency  in  the  Green’s  function  sum  rule.  In  turbulence 
problems  it  has  been  possible  to  construct  models  of  fictitious  systems  that  satisfy  the 
DIA  ex*Hly,  thus  guaranteeing  the  appropriate  positivity.  However  this  approach 
is  of  very  limited  usefulness  and  it  is  not  known  how  to  extend  it  to  other  approxi¬ 
mations. 

In  the  second  part  of  the  paper  we  have  obtained  the  previous  reeuhs  of  Hsssslminn 
via  a  prescribed  limiting  process.  The  particular  assumptions  inherent  in  the  process 
were  checked  snd  found  to  break  down  in  the  region  where  Olbtrt  and  McComas  found 
inootufaient  results  in  their  calculations.  We  suggested  that  the  oakmlations  might  be 
improved  by  solving  equations  for  £,  and  £|  self-consistently,  but  that  the  imple¬ 
mentation  of  the  procedure  aright  be  quite  difficult. 

In  con  clarion,  we  fed  that  the  formalism  of  the  type  discussed  here  should  be 
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Figure  12.  Tadpole  graph  with  response  function  integrated. 


Figure  13.  Tadpole  graph  with  correlation  function  integrated. 


studied  further  for  its  possible  application  to  the  internal  wave  problem.  The  prelimi¬ 
nary  calculations  indicate  ways  to  correct  for  deficiencies  in  earlier  calculations,  but 
further  vork  is  needed  to  formulate  positively  constraints  and  related  constraints. 

This  work  was  supported  in  part  by  the  Navy  SBC,  Johns  Hopkins  University, 
APL  601103,  and  by  the  Office  of  Naval  Research,  N00014-80-C-0840. 


Appendix 


We  discuss  briefly  the  role  of  the  Jacobian  (equation  (2.14))  in  cancelling  some 
tadpole  graphs.  We  will  illustrate  the  cancellation  by  considering  the  expectation 
value  of  A(k,  «*>),  which  should  be  zero.  There  are  two  graphs  contributing  to  {A  (k,  o»)> 
and  they  are  shown  in  figures  12  and  13.  If  we  work  to  lowest  order  in  B,  the  Jacobian 
gives  the  contribution 


<A«(k,»))jMoM»  -  2fd*1dM,g  (Al) 


where  (AA}t  is  given  by  (3.1).  Using  the  rules  given  in  the  text  it  is  easy  to  calculate 
the  contribution  of  figure  12  and  verify  that  it  cancels  (A  1).  The  graph  shown  in 
figure  13  vanishes  for  another  reason.  It  is  easy  to  see  that  it  is  proportional  to 


and 


D“(k  «  0,  o> > 


0)* 

BSZaO. 


When  these  graphs  are  inserted  into  more  complicated  graphs  they  will  oontinue  to 
give  no  contribution.  Similarly  more  complicated  tadpoles  (vertex  corrections  of 
figures  12  and  13)  will  continue  to  vanish. 
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ABSTRACT 


V 

In  this  paper  we  apply  a  formalism  introduced  in  a  previous  paper  to  write 
down  a  self-consistent  set  of  equations  for  the  functions  that  describe  the 
near-equilibrium  time  behavior  of  random  oceanic  internal  waves.  These  equations 
are  based  on  the  Direct  Interaction  Approximation.  The  self-consistent  equa¬ 
tions  are  solved  numerically  (using  the  Garrett-Munk  spectrum  as  input)  and  the 
results  are  compared  to  parameters  obtained  in  the  weak  interaction  approxima¬ 
tion  (W1A) .  The  formalism  points  out  that  an  extra  parameter  that  is  implicitly 
vanishingly  small  in  the  MIA  has  a  significant  effect  on  decay  rates  when 
computed  self-consistently.  Me  end  by  mentioning  possible  future  self-consistent 
calculations  that  would  improve  upon  our  own. 


1.  INTRODUCTION 

In  the  past  few  years  the  study  of  transfer  of  energy  in  the  internal 
wave  field  of  the  ocean  has  become  an  active  area  of  research.  Olbers  (1976) , 
Me Comas  (1977),  Pomphrey,  Meiss  and  Watson  (1980), and  DeWitt  and  Wright  (1982), 
have  done  calculations  in  the  resonant  interaction  approximation  (RIA)  of 
lifetimes  and  action  transfer  rates. 

These  calculations  assumed  that  the  interactions  were  weak  and  the 
lifetimes  of  particular  modes  were  long  compared  to  the  Intrinsic  frequencies. 
Because  of  the  diversity  of  scales  in  the  ocean,  this  approximation  is  invalid. 
The  large  scale,  near  inertial  frequency  waves  interact  strongly  with  the 
smaller  scale  waves.  Attempts  at  solving  this  part  of  the  problem  have  been 
made  by  Meiss  and  Watson  (1982)  and  Henyey  and  Pomphrey  (1982) .  In  this  paper 
we  present  a  calculation  using  the  formalism  presented  in  DeWitt  and  Wright 
(1982)  and  DeWitt  (1982).  The  idea  is  to  use  the  Direct  Interaction  Approx¬ 
imation  (DIA).  In  fact  we  will  use  additional  approximations  to  the  DIA.  The 
details  and  further  references  can  be  found  in  DeWitt  and  Wright  (1982)  and 
DeWitt  (1982). 

The  DIA  equations  give  an  effective  linear  theory  with  memory  terms  and 
forcing  terms.  The  memory  function  and  forcing  function  are  determined  by  a 
self  consistent  calculation  that  will  be  described  in  Section  II.  In  this 
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section  we  write  down  the  effective  linear  theory  as  that  will  enable  the 
reader  to  better  understand  what  is  being  calculated. 

Our  formulation  of  the  nonlinear  problem  is  the  same  as  that  of  McComas 
(1977)  or  Meiss ,  Pomphrey  and  Watson  (1979).  We  very  briefly  review  the 
equations.  A  Lagrangian  is  introduced  with  the  variable  being  the  displace¬ 
ment  field  £j(r,t) 


L(Q  -  2  ^j) 


'Pg?j6j3 


3P 


e  .  1 


32P 


3r^  +  2  5j*k  Sr^  +  *  *  * 


+  2  ar^rj  +  !  '  *  j 


£  is  then  expanded  in  terms  of  the  normal  modes 

V?-'>  ‘  Z/A  <  M  **' 

s  ■  ± 


A(t)  is  the  mode  amplitude  and  ”  £*r  +  (2^  t 


(1.1) 


(1.2) 


is  the  dispersion  relation 


< 


The  equations  of  motion  for  the  amplitudes  A  are  given  in  the  next  section.  For 
details  on  the  coefficients  Z  and  the  couplings  B  of  the  next  section  see  DeWltt 
(1982)  or  Olbers  (1976).  In  the  DIA  approximation  all  of  the  Greens  functions  an^ 
correlation  functions  are  identical  to  those  of  the  following  linear  system  of 
equations: 


7T  + 


(k.t.f)  A"(lc,t’) 


Ck,t,t’)  A+  Ck,t')^  dt'  -  f“  (k. 


t) 


75-  Oc.t)  +  y  <>.t,t')  A’  (k,f)  +  r_+(k,t,t')  A+  (k,f)^  dt* 


f+  (k,  t) 


(1.3) 


^(t, t ' )  is  a  function  determined  by  the  DIA  equations.  Casuality  demands 
r(t,t')  ■  0  for  t'>t.  We  are  treating  a  stationary  problem  so  T  depends  only 
on  t-t'.  f(k, t)  is  a  random  forcing  function  on  the  k**1  mode  due  to  all  of 
the  other  modes  and  any  external  random  forces.  In  the  DIA  approximation 
the  f 's  are  Gaussian  random  variables  with  zero  mean  (provided  that  the  external  forces 
are  also  Gaussian  with  zero  mean).  Their  correlation  function  is  given  by 


Fj(k,t-t’)  5s>_g,  +  l\  8'  (k,t-t*)  -  <f*(k,t)  f8'  (k’,t’)>  .  (1.4) 

Fq  is  the  correlation  function  for  the  external  part  of  the  forces  (1.4)  and 
calculated  from  the  DIA  equations.  It  is  convenient  to  work  with  the 
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Fourier  transform  in  the  time  domain,  then  F  has  the  following  representation 


pSS 


(k,<u)  - 


'»-*  [  e 


Oc)  +  i 


n_8(k) 


J  -  Ej  8'  (k,u>) 


(1.5) 


Here  v  is  the  damping  due  to  interaction  with  the  external 
environment  and  If  E^  were  Independent  of  u»  it  would  represent 

three  effects.  The  obvious  ones  are  a  finite  lifetime  and  frequency  shift.  In 
addition  however,  there  is  a  coupling  between  the  ±  modes.  Previous  calcula¬ 
tions  have  used  weak  interaction  theory  which  is  the  limit  of  the  DIA  equation 
in  which  E^  and  E^  are  Independent  of  b)  and  infinitesimal. 

In  the  following  section  we  introduce  the  correlation  function  for  the 

A's 


7~US8'  (k,u>)  <53(£+k)S(W)  -  <a£(0)>  aJ,'  (w')> 


(1.6) 


and  the  Green's  function  for  equations  (1.3)  , 

vQ  +  is£  -  iw  -  iEj8'  (k,o>) 


G88  (k,u)) 


-1 


(1.7) 


The  Inverse  is  a  matrix  inverse  of  a  2  x  2  matrix. 

The  only  additional  information  that  one  gets  from  the  DIA  equations 
other  than  the  correlation  and  Green's  functions  comes  from  an  Inspection  of 
the  calculation  of  E^  and  E^  where  it  is  possible  to  deduce  where  the  action 
in  a  mode  k  is  flowing  and  where  it  came  from.  We  do  not  present  any  of  those 
results  in  this  paper.  From  the  fact  the  DIA  equations  give  correlation  and 
Green's  functions  that  are  equivalent  to  a  linear  theory  in  which  different 


modes  are  uncoupled,  it  Is  clear  that  they  can't  give  a  completely  accurate 
picture  of  the  Internal  wavefleld.  The  effect  of  these  correlations  Is  lumped 
Into  finite  lifetime  and  frequency  shifts.  Thus  we  can  expect  some  improve¬ 
ment  on  weak  Interaction  theory,  but  the  Interaction  between  diverse  scales 
may  still  not  be  adequately  treated. 

We  find  substantial  differences  from  weak  interaction  rates  for  many  of 
the  internal  wave  modes.  That  Is  the  functions  E^(k)  and  Z^Ck)  are  quite 
different  from  those  of  weak  Interactions.  One  problem  that  immediately  arises 
in  this  case  is  the  validity  of  truncating  the  Lagrangian  in  Eq.  1.1  at  third 
order.  An  Eulerlan  formulation  would  remove  this  problem  but  introduce  others. 
In  this  paper  we  use  Eq.  (1.1)  and  ignore  higher  order  corrections  without 
attempting  to  justify  that  approximation. 

The  remainder  of  the  paper  is  organized  into  four  sections.  In  Section  II 
we  write  down  the  DIA  equations  and  our  approximation  to  them.  In  Section  III 
we  discuss  briefly  the  meaning  of  the  parameters  and  in  Section  IV  we  give  the 
results  of  our  calculations.  Section  V  gives  our  conclusions. 
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2.  THE  SELF-CONSISTENT  EQUATIONS 

In  this  section  we  write  down  the  specific  equations  that  will  be  used 
in  the  numerical  solution  of  the  problem.  The  derivation  of  the  more  general 
Direct  Interaction  Approximation  (DIA)  equations  is  given  in  an  earlier  paper. 
(Dewitt  and  Wright,  1982).  The  exact  DIA  equations  are 


£  j88  <k,u> 


■  I  £ , /Vw, 


»vsi 

$2’*2 


-8  82“®! 


6(k+k2-k1)6(aj+ti)2“UJ;L)Bk  _k  x 


(2.1) 


'*1  *2  -S‘P 


82®2 


\-k2-k  [iG  C-U  (k2,a)2)J 


and 


2**  (k,u)  ■  ~  V  f  d^k^d\2du^dw2 

1, . s’  J 


'1*  "1 
*2’s2 


“•  b1  s- 

6(k+k2+k1)<S(u>*w2+<»»1)Bk  kx  k  x 


(2.2) 


“•  *1  *2  Vl  *2  *2 

B-k  k*  -k2  U  <V“1>  0  {k2’“2> 
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where  6  is  given  by  Eq.  1.7  and  U  by 


t  ss. 

UM  <k,u>)  -  G  X(k,u>)  C5S|_s.F0(k,ui) 


8.8-  S's, 

+  S2X  X(k,ui)3  G  i-k,-<u] 


All  of  the  theoretical  information  about  the  interaction  of  the  internal 

®18283 

waves  is  contained  in  the  coupling  coefficients  B.  .  ,  which  in  turn  are  de- 
rived  from  first  principles  and  Include  some  approximations.  They  are  dynami¬ 
cally  defined  by  the  equations  of  motion  for  the  internal  wave  amplitudes  A: 


AJ(t)  +  inj(t)  +  vo(k)Aj(t)  +£  J d3k1d3k2:(k-jc1-jk2) 

1,2*^ 


8'*8l"82  8i  ®2 

Vk,-k2  \  \2 


f0(k,t) 


(2.4) 


The  goal  of  the  calculation  is  to  compute  the  matrices  and  in  a  self- 
consistent  manner.  All  of  the  Information  about  the  decay  and  correlation  of 
the  waves  is  contained  in  these  two  quantities. 

He  hasten  to  remind  the  reader  that  even  a  complete  solution  of  the  DIA 
equations  does  not  constitute  an  exact  solution  of  the  problem,  since  these 
equations  were  obtained  by  completely  ignoring  all  vertex  corrections  in  the 
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relevant  diagrams.  It  is  obviously  not  clear  a  priori  that  such  an  approach 
is  reasonable;  one  must  eventually  compute  at  least  the  leading  terms  in  the 
vertex  expansions  to  vindicate  the  procedure. 

As  was  pointed  out  in  the  previous  paper,  much  has  been  said  about  possi¬ 
ble  ways  to  solve  systems  of  equations  that  are  related  to  the  full  DIA  set  of 
equations.  However  the  only  numerical  attempts  at  solution  in  this  framework 
have  used  the  Weak  Interaction  Approximation  (WIA) ,  which  computes  and  Eg 

by  letting 


E1*E2 


(2.5) 


in  the  expressions  on  the  right  hand  side  of  equations  2.1  and  2.2.  The  only 

outcome  that  would  justify  this  approach  would  be  one  in  which  the  computed 

values  of  were  small  compared  with  typical  oscillation  frequencies.  As 

has  been  pointed  'nit  many  times  before,  this  is  not  the  case.  This  makes  the 

WIA  a  non-self-consistent  approach. 

Our  parameterization  of  the  DIA  equations  assumes  that  G(k,<*>)  can  be 

represented  by  two  poles,  one  at  u  -  ft^  -  iSj  and  a  second  at  ui  •  -f^-ia^. 

Symmetry  between  ±  k  is  assumed  as  well  as  between  ±  ft.  The  equations  are 

then  perameterlzed  by  the  pole  positions.  Also  we  require  G8®  (k,t  -  0+)  - 

<5  ,.  Finally,  there  is  the  qqestlon  of  where  to  evaluate  £-(&>) .  (Carnavale 

■  »“  ■  J. 

and  Fredericksen,  1982).  That  is  do  we  evaluate  it  at  the  pole  for  complex  u 
or  somewhere  on  the  real  axis.  We  chose  to  evaluate  E^(bi)  at  u  -  ft(k)  with  ft 
given  by  the  linear  dispersion  relation.  See  Section  3  for  further  comments. 

In  particular  we  choose 
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r1(k,u)  -  z^k.n) 


E2(k,ai)  -  £2(k,ft) 


(2.6) 


The  quantities  and  Z2  are  in  general  2X2  matrices  with  complex  compo¬ 
nents  and  so  superficially  appear  to  represent  16  real  parameters.  However 
using  the  constant-u)  assumption  and  the  general  symmetry  conditions 

H88'  (k,u)  -  [U"S*~S'  (-k,-u>)]* 


.ss'tk.u)  -  U8  8(-k,-u)> 


G88'  (k,a>)  -  [G-8*-8'  (-k.-a))]* 


(2.7) 


enables  us  to  reduce  these  to  7  independent  real  quantitites.  We  denote 


a(k)  =  Z7'(k,n) 

**  1 


c(k)  =  z^(k,n) 

(2  8) 

d(k)  5  Et+ok.n) 

^  4 


a(k)  5  E2  (k,n) 


The  first  symmetry  condition  implies  that  e(k)  is  real;  the  other  three  quan¬ 
tities  are  in  general  complex. 

By  invoking  a  further  assumption  we  can  reduce  this  set  of  parameters 
even  further.  We  procede  as  follows;  with  tbs  above  assumption  we  can  write 


the  Green’s  function  for  this  problem  as 


G(k,w)  ’  (w-n)  (u>-a) 


(w-a)  1 

\a>-0+a 


T) 


(2.9) 


In  writing  matrices  we  will  use  the  convention  that  s  ■  +  is  represented  by  the 
first  row  (or  column).  In  the  expression  for  G 


n  -  fi.-iaj 


=  -Or  -  iax  (2.10) 


flj  =  (fl-a^2  -  |c|2 


Here  £)r  Is  the  renormalized  "frequency";  and  a^  represent  the  real  and 
imaginary  parts  of  a  respectively.  He  note  that  in  the  weak  interaction  limit 
but  that  in  this  self-consistent  calculation  &r  may  be  either  purely  real 
or  purely  imaginary  depending  on  the  relative  values  of  (Q-a^)  and  |c|.  This 
will  be  explored  more  fully  later. 

Notice  that  in  this  particular  case  the  two  poles  of  the  Green's  function 
satisfy 


(mo)  -  -(me) 
(no)*  -  no 


(2.11) 


This  result  can  be  seen  directly  by  costparing  with  2.10,  but  more  importantly 
it  can  be  obtained  by  applying  the  symmetry  property  of  G  in  2.7  without 
knowing  the  analytical  expressions  for  n  sad  o.  This  is  useful  because  the 
sane  structure  appears  when  Green's  functions  with  more  poles  are  considered. 
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In  such  cases  it  is  not  generally  possible  to  obtain  closed  form  expressions 
for  the  poles. 

From  the  equation  for  G  we  can  compute  the  correlation  function 

U(w)  -  G(u)  S2  GT  (-u)  (2.12) 

where  the  superscript  T  denotes  transpose  of  the  matrix.  If  we  denote 

d(uj)  -  [(oj-n)  OiH-n)  (w-o)  (a*+c]_1  tJ(to)  (2.13) 

then  U(u>)  has  no  poles. 

Now  we  are  ready  to  impose  "experimental"  constraints  on  our  expressions. 
The  function 

U(k, t)  - y*du  e lwt  U(k,u>)  (2.14) 

Gives  the  two-time  correlation  of  the  internal  wave  amplitudes  and  so  the  quan¬ 
tity 


0**'(t-0)  -  3*0L,k  )  6  ,  (2.15) 

B  V  S,— S 

is  given  experimentally  by  the  Garrett-Munk  spectrum  (Garrett  and  Hunk  1972, 
1975).  We  take  this  as  an  additional  constraint  on  the  equations.  If  the 
expressions  for  U  are  made  to  satisfy  this  constraint  it  is  found  that 


d  »  -ic  (e/ax) 
3F  -  e/2aI 


(2.16) 


This  last  expression  Is  the  same  one  obtained  in  the  WIA.  We  see  that  Z2 
completely  determined  by  Ej  and  the  Garrett-Hunk  spectrum.  This  leaves  us 


U 


with  only  four  independent  reel  quantities  to  compute.  The  expressions  for  0 
under  this  approach  are  also  greatly  simplified: 

(w2  +Off)  (art-Jl-ap)^  +  a2  -  1  c  | 2  \ 

h  1  «•«> 

(w-ftha^)2  +  a2-|c(2  ^  (to2  +  no)  / 

The  great  benefit  of  this  approach  is  that  we  now  need  only  iterate  one  of  the 
two  DIA  equations  (the  equation  giving  £^) .  As  a  matter  of  fact  if  one  reviews 
the  steps  taken  to  derive  the  final  expression  for  6  and  D  it  is  not  difficult 
to  see  that  this  is  generally  true  whenever  £2  chosen  to  be  constant,  regard¬ 
less  of  the  form  of  £^.  After  the  equation  has  been  solved  self-consistently 
the  value  of  e  may  be  computed  from  the  £2  equation  to  determine  how  good  the 
ratio  e/2aj.  actually  gives  back  the  Garret -Hunk  spectrum.  If  we  use  the  equations 
in  Section  1,  it  is  easy  to  derive  the  following  balance  equation: 


Vk) 

^<k) 


-  2Vk) 


2Im£1(k)  - 


Z2(k) 


(2.18) 


If  there  is  no  forcing  or  dissipation  on  mode  k  due  to  noninternal  wave 
sources,  the  left  hand  side  is  zero.  In  this  case  if  the  right  hand  side 
is  nonzero  one  would  infer  some  inconsistency.  One  possibility  is  that  the 
Garrett-Munk  spectrum  ?S(k)  is  not  correct.  A  second  is  of  course  that  the 
theoretical  calculation  of  £^  and  £2  is  Invalid.  The  third  possibility  is  that 
there  is  either  an  effective  Fq  present  or  an  effective  Vq,  i.e.  mode  k 
might  be  driven  or  damped  by  direct  coupling  to  the  external  environment. 
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3.  PHYSICAL  SIGNIFICANCE  OF  THE  PARAMETERS 

Before  we  display  our  numerical  results  for  the  parameters  E^  and 
we  would  first  like  to  discuss  the  meaning  of  Che  parameters  introduced  in  the 
last  section.  To  do  this  we  recall  that  the  nonlinear  evolution  Equation  2.4 
may  be  equivalently  written  as  a  set  of  effective  linear  evolution  equations, 
which  for  the  constant  E  case  considered  in  this  paper  are 

A+(k,t)  +  (vQ+ifi-ia)  A+  (k,t)  +  iC*  A_(k,t) 

-  f+r k.t) 

(3.1) 

A~(k,t)  +  (uQ-iftfia*)  A"  (k.t)  -icA+(k,t) 

-  f"(k,t) 

The  functions  f+  and  f~  give  the  effective  driving  of  the  waves  due  to  all  other 
modes.  The  parameter  Vq  is  a  free  parameter  that  may  be  used  to  Inject  additional 
information  about  the  decay  of  internal  waves  due  to  interrction  with  external 
systems.  In  all  numerical  work  we  always  set  ■  0. 

The  parameter  a  represents  two  physical  quantities.  The  Imaginary  part 
of  a  represents  an  additional  effective  damping  of  the  internal  wave  mode 
due  to  losses  to  other  internal  waves  modes.  The  real  part  of  a  gives  the 
frequency  shift  of  the  wave.  The  parameter  a  is  not  new  to  this  type  of 
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analysis.  The  parameter  c,  however,  has  not  to  our  knowledge  appeared  in  this 

context  before.  We  see  that  the  formalism  compels  us  to  include  a  parameter 

that  represents  a  coupling  between  +  and  -  traveling  waves. 

The  pair  of  first  order  equations  for  the  amplitudes  may  be  written  as  a 

pair  of  second  order  equations  involving  only  one  type  of  amplitude  each. 

■f  «• 

If  we  ignore  the  f  ,  then  for  either  A  or  A  we  have 


A  (k,t)  +  2(vQ+aI>  A  (k,t) 


(3.2) 


+  [(fl-a^2  +  (Vg+a-j)2  -  | C | 2  ]  A  (k.t)  -  0 


The  solution  has  a  time  dependence 


A  ~  exp 


-(Vg+aj)  t  ?  \f\c\2  -  (H-ar) : 


(3.3) 


i  2 

If  |c|*'>((2-ar)  ,  then  we  refer  to  the  mode  as  overdamped.  In  the  underdamped 
case  the  damping  is  given  by 


V  "  V0  +  aI  '  (3.4) 

In  the  overdamped  case  the  damping  is 

V  "  V0  +  aI  4  Inri  (3.5) 
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We  will  define 


*  *1  +  |nri 

\  *  *!  ‘  ^r1 


In  the  overdamped  region 


and 


(3.9) 


(3.7) 


A  spike  disturbance  in  the  equilibrium  spectrum  will  contain  both  a  short  lived 
and  a  long  lived  component.  We  will  see  from  the  results  of  the  numerical 
calculations  that  it  is  that  measures  the  liftime  of  a  disturbance  for  a 
large  part  of  the  internal  wave  spectrum. 


4.  NUMERICAL  RESULTS 


In  Figure  (1)  we  show  a  contour  plot  of  log  |c/(Q-aD)|  vs.  log  (m)  on 
the  horizontal  axis  and  log(O)  on  the  vertical  axis.  Here  m  refers  to  the 
vertical  wavenumber.  The  line  of  critical  damping  has  the  value  zero; 
everything  to  the  left  of  this  line  is  in  the  underdamped  regime,  while 
everything  to  the  right  is  in  the  overdamped  regime.  It  is  instructive  to 
compare  this  plot  with  Figure  (2),  a  plot  of  log(a^(WI)/ft)  for  the  weak 
interaction  (WI)  case  using  the  same  wavevector  cutoffs  for  both.  Notice 
that  the  SC  overdamped  region  covers  the  entire  region  and  overlaps 

the  a1/n<l  region,  so  that  the  overdamping  takes  over  j  ust  at  the  place  where 
the  old  weak  interaction  calculation  begins  to  break  down.  What  is  interes¬ 
ting  is  the  extent  to  which  the  self-consistent  (SC)  calculation  reduces  the 
rates  in  the  large-m,  large-0  region.  Of  course  we  do  expect  the  SC  calcu¬ 
lation  to  reduce  the  damping  rates  in  this  region,  for  reasons  described 
below,  but  the  way  that  this  reduction  comes  about  is  rather  unexpected. 

Recall  that  the  rate  of  decay  of  a  spike  disturbance  introduced  into  the 
steady  state  ocean  is  given  by  and  c*s  in  the  SC  calculation.  (See  Eq.(36>) 
If  we  plot  log|tx^/ag |  in  Figure  (3)  we  see  that  in  essentially  the  entire 
overdaaped  region  cxL«agt  which  implies  that  ar(SC)r|0R| .  This  means  that 
the  decay  rates  are  much  smaller  than  one  would  expect  by  examining  the 
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scale  of  a^.  Unfortunately  the  same  device  that  causes  the  dancing  rates  to 
be  so  greatly  diminished  also  makes  it  more  difficult  to  obtain  accurate 
numbers  in  the  large-m  region.  A  small  change  in  the  parameters  (a,c)  in  this 
region  can  cause  a  disproportionately  large  change  in  (though  not  in  <Xg) . 
Since  it  is  (a,c)  that  are  computed  in  each  iteration,  factors  such  as  the 
integration  accuracy  make  it  difficult  to  control  such  changes  after  a  certain 
limit  is  reached.  We  have  solved  the  equations  by  an  iterative  scheme. 

After  many  iterations  the  results  do  not  change  much,  but  the  cancellations 
and  inaccuracies  in  the  integration  scheme  lead  to  some  changes  from  itera¬ 
tion  to  iteration.  In  the  next  two  figures  we  compare  the  from  two 
successive  iterations.  Figure  (A)  is  a  plot  of  the  percent  deviation  of 
the  input  values  of  compared  to  the  output  values  of  Sj. 

aICin)-ax(out) 


ACaj) 


100 


(A.l) 


ax(in)+ax(out) 

In  a  perfect  calculation,  of  course,  these  numbers  would  be  identical.  We 
see  in  fact  that  in  most  of  the  region  of  interest  we  have  convergenced  to 
within  a  few  percent,  and  it  is  only  for  very  large  values  of  m  that  we  reach 
values  of  8Z  to  10Z.  However,  if  we  compare  this  with  Figure  (5),  which  is  a 
plot  of 


logA  (c^)  =  log  100 


aLCin)-qLCout) 


o^ClnJ-to^Cout) 


(A. 2) 


we  see  that  the  error  Increases  in  more  quickly  than  it  do**  in  a^.  In  the 
large-m  region  this  means  that  we  only  know  ct^  within  a  factor  of  2  or  so. 

We  beaten  to  point  out  that  this  does  not  alter  our  conclusiona  about  the 
dramatic  effects  of  the  SC  calculation.  First  of  all  the  fact  that  the  values 
of  (a,c)  for  a  given  mode  do  not  depend  on  the  behaviour  of  modes  with  even 
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moderately  greater  values  of  (m,fl)  means  that  the  computed  values  in  most  of 
the  space  are  not  affected  by  the  uncertainty  in  the  large-m  extreme.  Second, 
we  will  see  shortly  that  compared  to  the  WIA  the  rates  for  the  SC  case  in  the 
large-m  region  are  reduced  by  orders  of  magnitude,  and  so  even  a  factor  of  2 
uncertainty  does  not  eliminate  the  qualitatively  new  behavior. 

Referring  back  to  (3.7)  we  see  that  a^«ctg  Implies  that,  for  example, 

"V 

A  (k,t)  -e  L  (4.3) 

so  that  we  may  ignore  and  consider  only  in  determining  the  long  time 
behavior  of  G.  For  this  reason  we  plot  log(aL/0)  in  Figure  (6)  with  the 
understanding  that  in  the  underdamped  regime.  It  is  clear  that  these 

rates  are  very  much  different  from  those  calculated  using  the  WIA.  The  rates 
for  most  of  the  spectrum  now  satisfy  aL/ft«(2ir)  \  but  there  still  exists  a 
substantial  region  for  which  a^/fl >  1,  mostly  in  modes  with  high  frequencies. 
This  agrees  with  our  general  expectations;  on  the  one  hand  we  know  from 
experience  that  the  internal  wave  field  appears  to  be  weakly  nonlinear, 
particularly  for  low  f requencies .  However  we  also  know  that  all  rates 
could  not  be  small,  otherwise  the  WIA  would  have  been  successful  in  the 
first  place.  We  explicitly  compare  the  WIA  decay  rates  with  SC  rates  in 
Figure  (7)  where  we  plot  log  (<x^/a^(WI)) .  For  small-m  and  saall-0  both 
calculations  give  roughly  the  same  results.  For  moderate  to  large  m  we  see 
substantially  reduced  rates  for  G.  Notice  though  that  in  the  small-m, 
large-0  region  the  SC  calculation  actually  increases  rates  by  a  great  deal. 

The  change  in  rates  for  the  moderate  to  large  a  region  is  easy  to 
understand  if  one  remembers  that  McComas  (1977)  and  McComas  and  Bretherton 
(1977)  identified  the  induced  diffusion  mechanism  to  be  the  dominant  energy 


transfer  mechanism  for  the  WIA  in  this  region.  This  means  that  energy  dif¬ 
fuses  from  a  given  large-wavevector  mode  to  a  nearby  mode  in  wavevector 

space  by  interacting  with  a  third,  small-wavevector  mode.  Since  in  the  WIA 
all  waves  are  assumed  to  have  Infinite  lifetimes  the  rate  of  energy 

exchange  is  only  governed  by  the  coupling  between  the  modes.  However  in 
the  SC  calculation  the  finite  lifetime  of  each  mode  is  built  into  the 
calculation.  This,  and  the  fact  that  the  decay  rates  depend  on  the  ampli¬ 
tudes  of  the  interacting  modes,  implies  that  the  rate  at  which  energy  is 
transferred  between  modes  decreases  as  the  waves  decay.  Therefore  we  would 
expect  the  lifetimes  of  the  modes  in  this  region  to  be  increased. 

We  would  now  like  to  consider  another  measure  of  the  consistency  of 
the  calculation.  Recall  that  we  have  found  a  correspondence  between 
e  and  a  given  by  (2.16).  Given  the  parameters  (a,c)  we  can  now  compute  the 
ratio  e/2a^.U(0) .  In  truth  the  degree  to  which  e/2a^U(0)«l  depends  both  on 
how  well  GM  76  represents  an  equilibrium  spectrum  and  on  the  reliability  of 
the  calculation.  It  is  not  possible  under  the  current  assumptions  to  differ¬ 
entiate  between  these  two  effects,  so  it  is  best  to  talk  about  the  degree  of 
equilibrium  of  GM  76  under  a  given  approximation.  (See  also  the  discussion 
at  the  end  of  section  2)  For  example.  Figure  (8)  shows  a  plot  of 
e/2a^U(0)  for  the  WIA.  We  see  that  for  large  frequencies  the  spectrum 
appears  to  be  in  a  high  degree  of  equilibrium,  while  the  lower  frequencies 
are  out  of  equilibrium  by  50Z  or  more.  We  compare  this  with  Figure  (9) 
which  shows  log(e/2a^U(0>)  for  the  SC  calculation.  We  note  immediately 
that  the  previous  discussion  concerning  the  accuracy  of  is  relevant  here, 
since  a  factor  of  2  in  the  high-m  region  can  greatly  affect  the  equilibrium 
balance.  The  equilibrium  nature  of  GM  76  may  be  difficult  to  determine  for 


large-m  in  this  SC  calculation.  However  in  general  we  notice  a  similar 
kind  of  deviation  of  6/2^13(0)  from  1,  on  the  order  of  50%  of  so,  although 
where  the  VIA  sees  depressed  values  the  SC  method  sees  increased  values. 

We  take  this  to  be  a  validation  of  the  computational  process,  given  the 
tentative  nature  of  the  GM  spectra,  which  Garrett  and  Munk  propose  as 
"straw  men"  to  be  knocked  down  by  the  next  generation  of  experiments. 

Next,  we  have  done  another  test  calculation  to  determine  how  much  the 
SC  calculation  might  affect  the  off-resonant  character  of  E^  and 
We  have  calculated  the  derivative  of  E^  at  In  Figure  (10)  we  plot 

log|^'+(n)n/E[+(fl')  |  (4.4) 

a  measure  of  the  rate  of  deviation  of  E^  from  its  value  at  if  this 

quantity  is  much  greater  than  1  then  the  constant  E^  assumption  breaks 
down.  If  we  compare  this  with  the  same  quantity  calculated  for  the  weak 
interaction  case ,  shown  in  Figure  (.11) ,  we  see  that  the  region  where  the 
relative  deviation  exceeds  1  is  greatly  diminished  in  the  SC  case.  It 
should  be  noted  in  passing  that  the  very  restrictive  nature  of  the  weak 
interaction  calculation  (only  allowing  modes  to  Interact  if  they  satisfy 
the  strict  frequency  resonance  conditions)  allows  a  new  class  of  interaction 
to  suddenly  "turn  on"  at  <d»2f,  and  this  accounts  for  the  large  derivatives 
around  that  frequency. 


5.  CONCLUSIONS 


In  this  paper  we  have  described  a  two  pronged  advance  in  the  study  of  the 
time  behavior  of  oceanic  internal  waves.  First,  we  have  applied  a  formalism 
which  allows  for  a  systematic  study  of  Internal  wave  interactions.  Second 
we  have  performed  the  first  self-consistent  calculation  of  oceanic  internal 
wave  parameters.  We  will  first  discuss  the  importance  of  the  formalism. 

Since  the  work  of  Olbers  (1974,  1976)  and  McComas  (1977)  showed  that  a 
naive  perturbation  theoretic  approach  to  Internal  waves  did  not  give  reliable 
decay  rates  the  possibility  of  improvement  on  their  calculation  has  been 
discussed,  Holloway  (1977,  1978  1979)  introduced  a  "frequency  filter" 
correction.  The  problem  is  that  by  unsystematically  introducing  physically 
motivated  corrections  one  leaves  open  the  possibility  that  other  corrections 
of  comparable  magnitude  have  been  ignored.  In  some  cases  the  combination  of 
all  of  these  possible  corrections  could  actually  yield  results  much  different 
from  those  obtained  from  keeping  only  one  correction.  We  believe  that  the  new 
diagonal  term,  c(k) ,  in  our  calculation  provides  an  example.  The  case  where 
a(k)«*flnite,  c(k)-0  for  all  k  corresponds  to  the  frequency  filter  correction. 
We  found,  however,  that  a  nonzero  c(k)  had  both  a  qualitative  and  a  quan¬ 
titative  effect  an  the  results.  First',  a  nonzero  e(k)  allows  for  the 
possibility  of  am  overdamped  regime,  and  in  fact  we  found  that  most  of  tbs 


Che  small  rates  caused  by  the  subtraction  of  aT  and  |  In  the  overdamped 
regime  give  us  a  much  different  picture  of  the  time  behavior  of  internal  waves 
than  would  be  given  simply  by  a^,  the  decay  rate  dictated  by  the  frequency 
filter  alone.  We  therefore  conclude  that  it  may  be  difficult  to  include  all 
significant  corrections  to  the  calculation  from  physical  motivations  alone, 
and  that  the  formalism  provides  a  systematic  way  to  indicate  the  relevant 
effects. 

The  numerical  results  of  this. paper  indicate  the  degree  to  which  the 
earlier  numerical  calculations  were  unreliable.  In  fact  we  conclude  that 
while  that  earlier  decay  rates  were  basically  correct  for  the  small  wave- 
vector — small  frequency  regime,  the  decay  rates  for  most  of  the  spectrum  were 
much  too  large.  Our  calculation  indicates  that,  except  for  high  frequencies, 
internal  waves  tend  to  decay  relatively  slowly.  We  point  out  that  the  fact 
that  some  region  of  the  spectrum,  is  still  predicted  to  decty  quickly  is 
another  validation  of  our  procedure.  Had  all  decay  rates  turned  out  to  be 
small  we  would  have  been  forced  into  the  impossible  conclusion  that  the  weak 
interaction  approximation  was  correct  after  all.  Any  future  improvements  on 
our  calculation  cannot  change  this  fact.  Further,  this  argument  tends  to  imply 
that  any  future  calculation  cannot  substantially  reduce  rates  in  the  high 
frequency  regime,  since  our  numbers  there  are  not  much  greater  then  the  point 
where  weak  interaction  theory  becomes  valid.  This  leads  us  to  believe  that  our 
rates  are  fairly  representative  of  reality,  or  at  least  represent  a  rough  lower 
limit  to  it. 

We  close  by  considering  hew  our  SC  calculation  might  be  improved  upon.  The 
moat  logical  next  step  is  to  study  the  effect  of  non-constant  and  on  the 
results.  This  might  be  dona  by  choosing  functional  forms  for  Z^Go)  that  allow 
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for  more  than  two  poles  In  the  Green's  function.  A  continued  fraction  repre¬ 
sentation  of  has  been  tried  in  other  contexts  and  would  probably  work  well 
here.  The  problem  with  going  to  a  larger  number  of  poles  is  that  it  is  no 
longer  possible  to  find  analytical  expressions  for  the  poles  and  so  some 
insight  might  be  lost.  Further  difficulties  arise  because  a  continued  fraction 
representation  allows  for  new  poles  to  appear  in  G  in  successive  iterations. 
Since  it  is  not  usually  possible  to  keep  all  of  these  poles  for  computational 
reasons,  one  must  determine  a  reasonable  way  to  discard  some  of  these  poles. 
Nevertheless  choosing  some  consistent  scheme,  like  keeping  the  slowest 
decaying  terms,  would  make  such  a  calculation  an  attractive  candidate  for  the 
next  level  of  calculations. 


This  research  is  supported  in  part  by  contract  number  N 0001 4-80-00840. 
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Figure  Captions 

1.  Contour  plot  of  log  |c/(  ft-a^  |  vs.  log(m)  on  the  horizontal  axis  and  log 

(Q)  on  the  vertical  axis.  The  contour  interval  is  0.4.  The  zero  line 

corresponds  to  critical  damping.  All  modes  to  the  left  of  this  line  are 
underdamped;  all  lines  to  the  right  are  overdamped.  Dashed  contours 
indicate  lines  with  negative  values.  The  horizontal  axis  is  marked  off  in 
equivalent  vertical  wavelenths  (in  meters)  and  the  vertical  axis  is  marked 
off  in  units  of  the  Coriolis  frequency  f. 

2.  Contour  plot  of  log  a^(Wl)/  ft  vs.  log(m)  on  the  horizontal  axis  and  log 

(ft)  on  the  vertical  axis.  The  contour  interval  is  0.2.  Dashed  contours 

indicate  lines  with  negative  values.  The  horizontal  axis  is  marked  off 
in  equivalent  vertical  wavelenths  (in  meters)  and  the  vertical  axis  is 
marked  off  in  units  of  the  Coriolis  frequency. 

3.  Contour  plot  of  log  lo^/o^l  vs.  log(m)  on  the  horizontal  axis  and  log  (ft) 
on  the  vertical  axis.  The  contour  Interval  is  0.4. 

4.  Contour  plot  of  A(a^),  the  percentage  deviation  of  the  input  value  of  a^ 
compared  to  the  output  value,  vs.  log(a)  on  the  horizontal  axis  and  log 
(ft)  on  the  vertical  axis. 

5.  Contour  plot  of  log  A(o^) ,  the  logarithm  of  the  percentage  deviation  of 
the  input  value  of  compared  to  the  output  value,  vs.  log(a)  on  the 
horizontal  axis  and  log  (a)  on  the  vertical  axis.  The  contour  interval  is 


6.  Contour  plot  of  log  jc^/Qj  vs.  log(m)  on  the  horizontal  axis  and  log  (ft) 
on  the  vertical  axis.  The  oentonr  interval  la  0.4. 


7. 


Contour  plot  of  log  ja^/a^(WI) |  va.  log(m)  on  the  horizontal  axis  and  log 
(ft)  on  the  vertical  axis.  The  contour  interval  is  0.5. 

8.  Contour  plot  of  e/ (2ax(WI)  U(0))  vs.  log(m)  on  the  horizontal  axis  and  log 
(ft)  on  the  vertical  axis.  The  contour  interval  is  0.4. 

9.  Contour  plot  of  log(e/ (28^(0))  vs.  log(m)  on  the  horizontal  axis  and  log 
(ft)  on  the  vertical  axis.  The  countour  interval  is  0.4. 

10.  Contour  plot  of  log  |r^(ft)ft/Z^+(ft)  |  vs.  log(m)  on  the  horizontal  axis 
and  log(ft)  on  the  vertical  axis,  for  the  self-consistent  calculation. 

The  contour  interval  is  0.4. 

11.  Contour  plot  of  log  | E^+(ft)ft/Z^+(ft) |  vs.  log(m)  on  the  horizontal  axis 
and  log  (ft)  on  the  vertical  axis,  for  the  weak  interactions  case.  The 
contour  interval  is  0.4. 
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Ha  have  atudled  the  1-dlnanslonal  Iteration  nap  associated  with 
the  even  quart  ic  polynomial  x^j  ■l+ax*  +  b  x*.  This  nap 
allows  a  awooth  transition  from  a  single  hunp  to  a  double  hunp. 
Bifurcations  and  higher-order  transitions  occur  as  wa  vary  the 
parameters  a,  b.  In  addition  to  the  usual  universal  bifurcation 
behavior  discovered  fay  Fsigenbauu,  we  find  a  new  universality 
class  of  bifurcations  which  la  associated  with  a  trlerltlcal 
point.  Trlerltlcal  points  serve  as  natural  boundaries  to 
Mlgenbauu  critical  lines.  For  the  quartic  nap,  the  trlerltlcal 
points  which  are  the  end-points  of  tha  original  Fsigenbaus  line 
are  (a,b)  -  (0,  -1.59490)  and  (-2.81403,  1.40701).  Associated 


as  well  as  two  Independent  exponents.  Tha  exponents  are  6i. 
7.2851  and  d£2'  -  2.8571.  At  thg  trlerltlcal  point,  we  can 
Introduce  a  universal  function  fT(x)  which  obeys 

fl^CffCx/Oj))  -  fj(x) 


with  the  scale  factor  cl  ■  -1.69031.  The  quartic  nap  has  a 
special  duality  transformation  (a,b)  -  (a',b*),  such  that  the 
two  mappings  axe  Intrinsically  related.  Tha  trlerltlcal  points 
which  are  dual  to  tha  above  pair  of  trlerltlcal  points  are 
located  at  (-3.18980,  2.54371)  and  at  (.95561,  -1.14981)  and  are 
Joined  fay  a  line  ddeh  la  the  dual  of  tha  original  fhlgenbaum 
Una.  There  are  an  Infinite  number  of  different  trlerltlcal 
points  Which  form  at  least  a  Cantor  ant. 


i.  umooocnoH 


decently,  Ihlgenbaum  has  shown  that  tha  bifurcation  sequence  la  a  single  hump 
Iterative  map  x^.*  «  f(x_)  with  a  quadratic  peak  obeys  a  universal  behavior.1*2 
One  can  understand  this  behavior  from  a  xenermallsatiam-growp  point  of  view. 

At  tha  limiting  point  of  a  bifurcation  sequence,  tha  2”  Iteraetlon  of  tha  map 
with  appropriate  sealing  approaches  a  universal  function  fp(x).  In  tha  neigh¬ 
borhood  of  this  universal  function,  It  appears  that  therm  Is  only  one  relevant 
direction,  dang  Which  tha  eigenvalue  6(-  4.6692)  la  larger  tha a  one.  Wo  have 


)  with  a  quadratic  peak  obeys  a  universal  behavior. 


la  a  single  hunp 

•rul  or.  »* 


recently  studied  a  norm  general  map  described  Ip  a  quartic  polynomial 
a**!  •  f(0  il  +  ir  +  b  x?.  An  Inpottsat  property  of  tha  quartic 
nd  f(x)  Son  describe  both  a  slngle-huap  and  a  Amifats  lamin  nan.  By  i 


a  a»b 
to  a  t 


exist  in  the  lteratlana  of  a 
deuble-hsnp  map  is  that  of  do 


of  a  doubly  ay able  opal*  la  tha  a^b 


a  eat  la dose  a  amaoth  traaaltlon  from  a 
▼lea  versa.  It  ia  easy  to  sea  that  a 
at  kinds  of  stable  of*laa  which  do  net 
-hoop  nip.  da  important  osmeept  in  a 
fifela  cycles. '  theme  ate  cycles  whoea 
al  pooh  and  a  side  peak,  in  tha  vieimlty 


stability  of  the  sue  cycle.  Away  froa  doubly-stable  cycles,  however,  the 
peaks  control  different  cycles  which  are  unrelated  in  x-apace.  Thus,  the 
regions  controlled  by  these  different  peaks  become  dynamically  coupled  at  a 
doubly-stable  cycle.  This  kind  of  dynamical  coupling  Is  a  new  phenomenon 
which  cannot  occur  in  a  single-hump  map.  Its  -presence  modifies  the  bifurca¬ 
tion  processes.  The  limiting  point  determined  by  a  bifurcation  sequence  of  2** 
doubly-stable  cycles  Is  associated  with  a  fixed  point  in  function  space  with 
two  relevant  directions  in  analogy  with  Felgenbaum's  discussion  of  the  qua¬ 
dratic  case.  This  limiting  point  also  serves  as  the  end-point  of  a  critical 
line  which  describes  the  usual  fhlgenbaun  behavior.  This  type  of  fixed  point 
is  known  as  a  trlcrltical  point  in  phase-transition  language.  The  trier it ical 
point  discovered  here  has  two  universal  exponents  6^  '-7. 2851,  **2.8571. 

At  the  fixed  point,  the  map  is  self-similar  near  the  peaks  after  2”  Iterations. 
Just  as  in^ the  Jhigenbaum's  case.  It  is  possible  to  introduce  a  universal 
function  ff(x)  at  a  trlcrltical  point  which  obeys  the  same  functional  equation 

Ojf^f^x/Of))  -  f£(x) 

but  with  a  different  scale  fac\.  Jr  a_  -  1.69031.  Our  f l£(x}.can  be  expressed  as 
a  power  series  in  x  .  We  have  discovered  an  Identity,  5^  ^  -  *;|,  relating  one 
of  the  exponents  to  the  scale  factor.  This  Identify  is  a  member  of  a  whole 
class  of  identities  which  exist  at  the  fixed  point  of  an  iterative  nap. 

II.  QOARTIC  MAPS 

(A)  Single  vs  multiple-hump  maps 

lie  study  the  one-dimensional  iterative  naps  generated  by  an  even  quart 1c 
polynomial 

*irl-l  “  £<*n5  5  1  +  *  +  b  *n  *  <2*1) 


figure! 


Iterative  asps  generated  by  the  even  quart lc  polynomial  (2.1).  the  meg 
im  M  fa  advent  aentapaait  to  a  fix)  with  (a»h)  ia  that  quadrant. 


Depending  on  Che  values  of  a  and  b,  f(x)  can  have  either  a  single  huap  or  a  3 
double  huap.  See  Fig.  1.  R>r  region  defined  bp  a  >  0,  b  >  0,  f(x)  can  at 
aost  give  rise  to  a  stable  fixed  point.  For  a  <  0,  b  <  0,  the  function  f(x) 
has  a  single  quadratic  huap,  and  Felgenbata  universality  holds.  In  one  para- 
aeter  single-peaked  aaps,  the  Unit  point  of  infinitely  bifurcated  2"  cycles 
is  a  single  point.  In  two  paraaeters  single-peaked  aaps  that  point  becoaes  a 
line  that  we  refer  to  as  the  Ihigenbaua  line.  For  the  regions  a  <  0,  b  >  0 
and  a  >  0,  b  <  0,  f(x)  has  either  two  peaks  and  one  valley  or  vice  versa.  In 
the  following,  we  shall  refer  to  the  peak  or  valley  at  x  *  0  as  the  central 
peak,  and  the  peaks  or  valleys  at  x  *  0  as  the  side  peaks.  The  existence  of 
both  central  and  side  peaks  inplies  that  f(x)  nay  develop  independent  itera¬ 
tion  regions  as  Indicated  in  regions  I  snd  II  in  Fig.  2. 

f(x) 


X 


Figure  2 

A  asp  for  a  <  0,  b  >  0.  Rota  that  regions  I  and  II  are 
independent  iteration  regions. 

In  region  I,  f(x)  describes  a  single-heap  asp  with  a  quadratic  peak  at  x  -  0. 
Re  shall  rater  to  region  I  as  the  region  controlled  by  the  central  peak.  In 
tegloe  n,  f(x)  describes  an  inverted  single-heap  asp  with  the  peak  at  x  - 
^a/b.  Rs  shall  refer  to  n  as  the  region  controlled  by  the  side  peak. 

Mace  the  tee  side  peaks  have  equal  height,  they  asp  into  the  saae  x  (-1- 
a'/db).  React,  they  cannot  both  be  aaabers  of  an  R-cyela.  Thus,  ne  have  at 
aaat  oaa  aeat  stable  region  controlled  by  the  side  peaks.  The  existence  of 
tan  independent  regions  any  also  appear  la  the  case  a  >  0,  b  <  0. 


(1)  Duality  treneforaatlon 

ter  the  study  of  cycles  and  their  sequences,  as  only  need  to  none 1 dot  the 
region  controlled Vy  the  austral  peak.  f*  cycles  Ueoclated  with  the  Vide 
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peaks  can  be  obtained  froa  those  with  the  central  peak  by  a  duality  transfor¬ 
mation  described  below. 


Consider  the  quartic  asp  (2.1).  We  can  rewrite  it  as 

*M-l*1  +  *4+b4»  b*° 

* 1  -ft ♦  bi!b ♦  •  «•» 

Thus ,  wa  can  factor  the  quartic  asp  as  the  product  of  two  quadratic  maps : 

*n  -fE  + A  •  <2*3) 


*afl  "  1  -  b  ^  • 

He  can  now  visualise  the  original  sapping  as 


*i  *  h  *  *2  *  h 


♦  x«  ♦  •••• 


(2.4) 


(2.5) 


It  is  easy  to  see  that  if  the  x-napping  has  a  stable  N-cyde,  so  does  the 
(-mapping  and  vice  versa.  The  mapping  froa  ^  to  £p^  is  sinply 

^1  -$S+  U  -^+  b  tit  *  <2*6> 

After  a  rescaling,  we  can  put  the  (-napping  (2.6)  into  the  standard  fora 
-  1  +  a*  (J  +  b*  4  <*-7> 

with 

-  (4b-a2)[8ab  +  (4b-a2)2!  (2#8) 

32b1 

and 

.  :  ...  2.2,5i 

(2.9) 


b.  .  I8*h  +  Cjtr-*2)2]3  . 
(8b)4 


It  is  easy  to  verify  that  the  side  peaks  of  x  are  napped  into  the  central  peak 
of  (,  kd  the  side  peaks  of  (  are  napped  into  the  central  peak  of  x.  The 
duality  relation  is  reciprocal:  The  duality  transformation  of  (a',b')  is  the  . 
original  (a,b). 

(C)  Moat-stable  cycles  and  doubly-stable  cycles 

Ibr  a  quartic  map,  the  aost-stable  N- cycles  (x^Xj...,^)  are  described  by 
f(xt)  -  x1+1,  -  Xj,  i  »  1,2,...,N  (2.10) 


(2.11) 


An  H- cycle  ia  post  stable,  if  either  the  central  peak  (x  -  0)  or  one  of  the 
side  peaks  (x  -  ±/-a/Zb)  is  a  neaber  of  the  cycle  x^  .  The  polynomials 
associated  with  a  stable  H- cycle  give  rise  to  a  line  in  the  a-b  plane.  In 
fig.  3,  the  solid  lines  represent  the  loci  of  aost-stable  2-cycles,  fractions 
whose  peramsters  satisfy  a  +  b  +  1  •  0  have  most-stable  2-cydes  associated 
with  the  central  peak.  The  curved  lines  are  the  lod  of  points  of  aost-stable 
2-cyde  associated  with  the  side  peaks.  It  is  assy  to  see  that,  with  the 
exception  of  the  origin,  the  most  stable  H-cyde  trajectories  associated  with 
the  central  peak  never  intersect  among  themselves,  neither  do  tht  most  stable 
trajectories  associated  with  the  side  peaks.  However,  the  trajectories  asso¬ 
ciated  with  the  central  peak  do  intersect  those  associated  with  the  side 
peaks.  Actually,  they  intersect  in  two  completely  different  ways.  At  inter- 


sections  A  mid  A'  In  fig.  3,  the  none  (Ublt  2-cydes  pass  through  both  tho 
control  peak  and  one  of  the  aide  peaks.  He  call  this  kind  of  cycle  a  doubly- 
stable  cycle.  Two  trajectories  which  Intersect  at  a  doubly-stable  cycle  are 
“dynamically  coupled"  at  the  Intersection.  At  Intersections  B  and  B' ,  the  2- 
cydes  are  controlled  by  different  peaks  Which  are  not  related  t  soma  regions 
of  x  are  attracted  to  one  cycle  and  some  to  the  other.  He  refer  to  this  kind 
of  Intersection  ss  "dynamically  independent".  Note  that  doubly  stable  points 
A  and  A'  are  dual  to  each  other.  The  concept  of  dynamical  coupling  la  Invar¬ 
iant  under  duality  transformation. 


figure  3 

The  solid  Unas  represent  the  most-stable  2  cycles,  the  crosses 
denote  trlcrltlcal  points,  and  the  dashed  lines  describe  the 
fhlgenbaum  Unas. 

111.  TKICXmCAL  POURS 

(A)  Bifurcation  along  b-axls 

If  we  vary  b  (negative),  beeping  a  small  and  negative,  ms  encounter  the  usual 
bifurcation  sequence,  studied  by  Kilgeabeun.  If  me  vary  b  (negative),  keeping 
jt  small  and  positive,  me  find  that  the  map  goes  through  three  different  most 
stable  2-cyclee  before  bifurcating  Into  a  4-cycle,  iheee  2-cycles  are  con¬ 
trolled  by  the  three  different  peaks  as  shews  la  fig.  1.  Iheee  phenomena 
persist  no  met ter  horn  small  the  value  of  a.  The  sans  phenomsna  repeat  them¬ 
selves  at  higher  2*  cycles.  Thus,  there  mre  qualitative  differences  between 
bifurcations  on  the  two  different  sides  of  the  b-asls.  When  a  ■  0  and  b  <  0, 
the  curve  f(x)  Is  single-lamped  nth  a  quartlc  maximum,.  Vs  find  a  ssqnsnrs  of 
bifurcations  which  always  occur  at  the  doubly-stable  2*  cycles s 
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The  2*  cycles  hsve  «  Halting  point  st  ba  •  -1.59490,  in  close  analogy  to  the 
quadratic  asp  studied  by  Feigenbaum.  However,  because  our  asp  has  a  quart lc 
(rather  than  quadratic)  naxlnun,  It  is  to  be  expected  .that  there  will  be 
different  critical  exponents,  as  we  now  discuss.  This  Halting  point  has 
several  Interesting  properties  which  we  shall  discuss  below. 

(B)  Critical  exponents  and  the  universal  function 

There  are  two  relevant  elgen-dlrectlons  at  the  critical  point  (0,  -1.59490). 
Along  a  relevant  elgen-dlrectlon  and  for  large  N,  we  have 

bg-tu with  |6f|  >  1  (3.1) 

«T 

where  6y  Is  the  critical  eigenvalue  *1988  this  elgen-dlrectlon.  One  elgen- 
dlrectlon  Is  along  the  b-axia  where  di1'  -  7.2851.  The  other  elgen-dlrectlon 
Is  (1,  -1.2347)  where  -  2.8571.  Note  that  the  Rsigenbaun  fixed  point  has 
only  one  relevant  direction  with  6«  ■  4.6692.  In  phase-transition  language,  a 
critical  point  with  two  relevant  directions  itoich  also  serves  as  the  end-point 
of  a  critical  line  Is  known  as  a  "trlcrltlcal"  point.  4s  we  shall  see,  our 
fixed  point  is  Indeed  a  trlcrltlcal  point.  See  Discussion  Section  below. 


it  the  critical  point  (0,  -1.59490),  the  function  Is  self-slallar  near  x  -  0 
to  within  a  proper  rescaling.  The  scale  factor  here  Is  Oj  ■  -1.6903.  Just  as 
1|  the  Balgenbaue  case,  the  trlcrltlcal  point  has  Its  own  universal  function 
fT(x)  obeying  the  sane  functional  equation 

^f*<2)(x/qr)  -  fj(x),  0y  ■  -1.6903  •  (3.2) 

In  the  neighborhood  of  x  -  0,  f£(x)  can  be  represented  approxlaately  by 
f^x)  -  1  -  1.83413  x4  +  0.01301  x8 

♦  0.31188  x12  -  0.062035  x16  +  ••••  .  (3.3) 


(C)  An  Identity 

The  second  exponent  ft!2'  and  the  scale  factor  cl,  are  related  by 

42)  -  4  .  (3.4) 

This  relation  Is  exact, ^ as  can  be  shown  as  follows:  Let  h(x)  be  the  eigen- 

(3.5) 

(3*6) 

(3.7) 


Cta  M  CSX 

vector  belonging  to  Jj1'.  If 

f(x)  -  +  «¥*>»  0  <  t  «  1  , 

then  iteration  under 

lU)  *  f(2)(xf(l)) 

leads  by  definition  to 

g(x)  -  fj<*>  ♦  42>«h(*)  ♦  0 («2)  . 
Since  h(0)  ■  0  (because  f(0)  -  1),  we  nay  taka 


(3.8) 


h(x)  ■  +  *4*^  +  ®6*^  +  ••••  * 

Direct  Iteration  of  (3.6)  gives  (c^  -  l/f*(l)) 

of^W)  -  -*»(l)fj(«)  + 

♦  [HD  x  £(£-)  +  h(i-)]  ,  (3.9) 

which  is  a  functional  equation  for  h(x).  It  is  easy  to  verify  that  the  right- 
hand  side  vanishes  at  x  -  0.  Only  the  last  tern  contributes  to  0(x2)  and  by 
equating  both  sides  of  (3.9)  one  finds 

42)  -  fj'd)/®*  •  (3.10) 

One  evaluates  4  (D  by  differentiating  (3.2)  with  respect  to  x  and  then 
setting  x  ■  0.  The  result  is 

fj'(l)  -  oj  .  (3.11) 

which  gives  (3.4).  Hie  other  eigenvalues  belonging  to  directions  out  of  the 
x*,  x  ,  x  ...  subspace  are  irrelevant.  Eq.  (3.4)  is  an  exaapla  of  a  general 
class  of  identities  which  we  discuss  elsewhere. 

(D)  The  trlcrltlcal  point  and  doubly-stable  cycles 

• 

Our  present  trlcrltlcal  point  (0,  -1.39490)  is  the  bifurcation  Unit  of  a 
sequence  of  doubly-stable  2”  cycles.  This  turns  out  to  be  a  general  property 
of  a  bifurcation  trlcrltlcal  point.  We  have  used  this  property  to  deteralne 
nsny  other  trlcrltlcal  points  which  do  not  lie  on  the  b-axis. 

IV.  DISCUSSION 

In  the  a  <  0,  b  <  6  region,  we  always  have  a  slngle-huap  function  f(x).  If  we 
increase  the  paraneters  -a,  -b,  we  encounter  an  infinite  sequence  of  bifurca¬ 
tions,  as  described  by  Ihigenbaun.  The  Uniting  points  of  these  infinite 
bifurcations  now  lie  on  a  line  in  the  a-b  plane  which  we  refer  to  as  a  Phlgen- 
beun  line.  The  trlcrltlcal  point  (a,b)  -  (0,  -1.59490)  .serves  as  a  natural 
boundary  to  the  original  Phigenbaum  line.  When  we  extend  the  original  Feigen- 
baun  line  in  the  other  direction,  we  find  that  it  terni nates  at  another  tri- 
crltlcal  point,  located  at  (-2.81402,  1.40701).  One  can  check  easily  that 
this  trlcrltlcal  point  is  also  the  Unit  of  a  sequence  of  doubly-stable  2s 
cycles  and  that  it  has  the  sans  critical  exponents  6^'  *®d 

Und«r  the  duality  transformation,  a  trlcrltlcal  point  naps  into  a  trlcrltlcal 
point  end  a  Ihigenbaun  critical  Una  naps  into  a  Ihigenbaun  line.  The  dual 
Images  of  the  previous  trlcrltlcal  points  are  located  at  (-3.18980,  2.54371) 
and  at  (0.95561,  -1.14981).  We  have  plotted  the  dual  transformed  Ihigenbaun 
line  in  Pig.  3. 

Since  bifurcation  occurs  after  each  stable  cycle,  and  since  there  are  an 
infinite  nuaber  of  different  stable  cycles,  there  are  an  infinite  nuabar  of 
Mlgenbaua  Unas  in  the  e-b  plane.  Hence,  there  must  also  be  an  infinite 
nuaber  of  trlcrltlcal  points  in  the  a-b  plus.  Indeed,  we  are  able  to  show 
that  between  the  critical  points  (0,  -1.59490)  and  (.95561,  -1.14981),  there 
are  an  infinite  nuaber  of  trlcrltlcal  points,  forming  a  Canter  set.  Ws  shall 
discuss  the  geometrical  meaning  of  the  Center  set  elsewhere. 
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Since  toe  seepage  to  controlled  bp  toe  else  of  toe 
•■to,  *e  expect  tost  toe  vafc£rdtfrr%iiilt  Oie  gaps 
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lag  M«|in».  Figures  6—8  give  actual  graphical 
repress stations  of  p  in  the  neighborhood  of  a  3-1 
transition. 

Mow,  lat  oa  determine  toe  condition  tor  toe  3  - 1 
chaotic  transition.  As  we  have  mentioned  earlier, 
the  3*1  transition  occurs  when  points  D  and  E 
coincide.  Since  D  is  f*iP),  and  since  point  E  maps 
into  the  condition  for  D  and  E  to  coin¬ 

cide  is  tor  /(D) -/(F)- F.  This  then  implies 

fV>)  •/*&).  (3.4) 

Note  that  the  points  B,D,  F  form  an  unstable  three 
cycle  at  the  transition.  In  terms  ot  Ft*)  ■/*(*), 
we  can  rewrite  the  above  equation  as  the  universal 
equation 

FW-F'lP)*  FiP).  (3.5) 


Even  though  the  3—1  aad  S— 1  chaotic  transitions 
arise  from  different  mechanisms,  they  obey  the 
same  universal  equation.  The  only  difference  is 
in  the  definition  of  F:  F  is/*  for  the  3-1  transi¬ 
tion,  and  is  /*  for  the  3  - 1  transition,  we  can 
gensraUae  the  shoe*  rtsautt  to  an  arbitrary  hmisi- 
tion.  tot /»Xx(l  ~x),  (he  3-1  chaotic  transition 
occnrvat 

1*3.856800  883  477766*..  . 

».  Btrteeresf  Aslseder.  W*  now  discuss  toe  un¬ 
iversality  el  the  *»  to  r-*  transttfcma.  In  Table  I, 
weeeetfaat  as*- «,  a.  ipphOaches  the' iqifiiite 
two-cyele Mfttrcntice pofcft. - ... 
from  atmve.  Ia  fact.  1.  shtlsflea 


1  i 

X,-X.+y 


WMMjMlffft., 


Vers >8,  the 


than  enqptwt  is ; 


FU)mf%d  ■» 

hoed  at  *«lPs 


hyfW*- 
Mrler  of 


•*  ■■!*'„■*}  V.'V-frf’)  S  *«!'  " 


we  arrive  at  an  invariant  function  G(y)  at  large  N. 
This  invariant  function  G  Is  a  linear  combination 
of  the  invariant  functions  g  and  k  introduced  by 
Fei#enbaam.'  Note  that  G(0)»G(1)*=0,  and  G(l/3) 

« 1.  The  Invariant  function  G( y)  is  shown  in  Fig. 

10. 

IV.  FOLDING  AND  THE  LYAPUNOV  EXPONENT 

A  useful  measure  of  the  degree  of  mining  of  a 
dynamical  system  is  provided  by  the  Lyapunov  ex¬ 
ponent  For  a  one -dimensional  chaotic  system, 
the  Lyapunov  exponent  is  simply  related  to  the 
folding  phenomena.  We  consider  a  simple  one- 
dimensional  mapping 

***"/(*)■  (4.1) 

The  chaotic  reglon  ls  characterised  by  the  fact 
that  the  separation  of  two  neighboring  points  (*,,  y„) 
will  increase  exponentially  an  we  iterate  (4.1)  a 
large  number  of  tlmp.  For  an,  infinitiemal  separ¬ 
ation,  we  refer  to  toe  average  logarithmic  in- 


as  ths Lyapunov  exponent  The  exponent  p.  ia  pos¬ 
itive  (a>0i  in  toe  obastto  region.  Let  ua  start 
from  a  small  neighborhood.  After  many  iterations, 
toe  neighhpxhoed  will  enlarge  and  eventaally  cover 
a  oompnstapaoe,  Furlher  Mnratoooa  will  amp  this 
,  1  utojact  space.  .Into,  ttsaKi  Eiataown  that  there 
<WWM*«MA  Mnqtfr  bmetl  np  »(*)  saeh  torn  PM*r 
pr®vi*SAtoe  peeMhlUfy  measure  tor  finding  an 
arbitrary  point  inside  toe  region  4*  after  many 
itarattona  of  Bq.  (4.1).  The  probability  function 
4>{x)  is  invariant  under  toe  iteration  and  hence  is 
the  eigenfunction  with  unit  eigenvalue  of  toe  linear 
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y)+(y)dy<H*), 


Fold 

_ _ 

*  PE 


A(*.y)»8f*-/(y)J.  (4.4) 

We  have  learned  from  Lemma  8  of  Sec.  IX  that  the 
support  of  p  may  consist  of  N  disjoint  regions. 
The  iteration  of  £(x,y)  {or  (4.1)J  wiU  map  points 
from  one  region  to  another,  sad  bring  Die  point 
back  to  the  original  region  after  N  iterations. 
Thus,  if  we  consider  the  iteration  defined  by 


F(x)»/'W-/(/*-‘U)),  (4.8) 

we  encounter  iterations  which  mff  a  region  back 
onto  Itself.  Since  the  Lyapunov  exponent  ji  associ¬ 
ated  with  /  is  related  to  the  exponent  nK  associated 
with  F  via 

(4.7) 

we  only  need  to  consider  F  and  n*  involving  simp¬ 
ler  mappings.  A  typical  distribution  function  is 
shown  in  Fig.  11.  Let  AB  be  a  region  which  maps 
into  itself  under  F.  The  end  point  B  is  mapped 
onto  point  A,  end  point  A  and  an  interior  point  D 
are  both  mapped  onto  £  [F(A)»  £(£)»£],  and  point 
P  is  mapped  onto  point  £.  This  mapping  may  be 
achieved  by  folding  PD  on  PA  and  then  stretching 
BP  to  the  original  length  AB  (see  Fig.  18).  For 
more  complicated  mappings,  it  may  be  necessary 
to  fold  the  line  segment  more  than  once.  All  the 


►  p  (b) 


/  \ 


FIG.  18.  uTT**f  Involved  In  a  chaotic  One  region: 

(a)  Adding  of  AN  about  point  P  gives  lb).  Stretching  BP 
completes  the  mapping  shown  in  ip). 

points  In  segment  AE  come  from  BD  and  the  map¬ 
ping  is  1  to  1.  We  shall  call  this  segment  AE  a s 
nonoverlapping.  The  points  in  segment  EB  come 
from  both  AP  and  PD  and  the  mapping  is  2  to  1. 

We  call  this  segment  EB  as  overlapping. 

We  denote  the  density  funetkHi  g  within  segments 
AP,  PD,  and  DB  as  <pAP(x),  aiid  <PDt(x), 

respectively,  obeying 

a  for*e(AJ>) 

|0>  *  outside  (AP)  (4.8) 

etc.  Then,  we  have 

KK,Pd»~  ^a*  •  (4.9) 

**♦**•(♦ ,)**,  (4.10) 

K’*nm  <*»>».  (4.11) 


*»-(*, >*■  +  (*.)«».  (4.12) 

Equations  (4.10)  and  (4.11)  may  be  taken  as  the 
definitions  of  d,  and  P,.  In  the  following,  we  shall 
express  the  Lyapunov  exponent  as  a  function  at 
p,  Pi,  had  g,.  The  d,  and  A,  for  the  p  of  Fig.  11 
are  ahown  in  Figs.  13  sad  14. 

S  ip  more  convenient  to  denote  the  separation  of 
two  neighbor  lag  points  as  P{x)dx  rather  then  as 
dx.  To  Understand  the  definition  of  m  using  p, 
eeaeigor  n  particular  point  x,  and  tta  neighborbood 
dx,.  Map  It  a  sufficient  number  of  times  N  such 
that  it  returns  to  its  original  aeigbborfcood.  If  tbs 
new  length  is  dx„  the  magnification  is 


The  anlnrgi 
factor  p{x) 


FB.  li.  Mteot 

4(a)  ter  k»s.n. 
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FK3. 18.  Density  function  0,  (M  text). 


ping  region  due  to  the  conservation  of  the  prob¬ 
ability  under  the  mapping  x=F{y),  giving 


dUVfx*  <p(y)dy 

(4.13) 

or 

...  _  i>(x)dx  . 

(4.14) 

However,  in  the  overlapping  region,  *  can  come 
from  two  or  more  different  points,  such  as  y, 
and  y„ 

(4.15) 

with 


i>ilx)dx‘  dfaVfy, ,  (4.16) 

♦,(*)***  ♦(*,)**.  (4.17) 

The  magnification  factors  doe  to  the  mapping. 
I,-*  and  y,-x  are  no  longer  one,  and  are  given 
by 


(4.18) 

(4.19) 


The  geometrical  average  of  the  magnification  fac¬ 
tor  m  approaches 


me 


ir(n-(7,f. 


and  the  Lyapunov  exponent  becomes 


Me*  lam -lim  -  yilnmCv.) . 

m~>  »  i 


(4.20) 


To  compute  the  limiting  average  magnification,  we 
note  that  the  weighting  factor  at  aqy  point  y  in  pro¬ 
portional  to  its  density  d(yWy.  Thus,  we  can  re¬ 
place  E,  in  (4.20)  by  an  integral  over  +(y)dy  and 
"by  /dyd( y),  giving 

l^./ssy  »{?>»»"{?> 

KatyHy) 

In  the  numerator  of  the  above  equation,  we  can 
separate  the  y  Integration  into  three  regions:  Re¬ 
gion  DB  which  maps  into  the  nonoverlapping  region 
(AS)  of  x  gives  no  contribution  since  lam*  0.  The 
other  y’s  are  separated  into  y,  and  *  regions  (A P 
and  PD),  and  they  both  map  into  the  overlapping 
region  (SB)  of  x.  The  Lyapunov  exponent  becomes 

p,  -  +  /dy,»(y.)lnm(y,) 

The  expression  has  a  nicer  form  when  expressed 
in  terms  of  x.  Using 


no.M.  Density  Amotion*,  (me  teat). 


we  finally  obtain 
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as 


/Bsdxp(x) 


The  expression  in  the  last  large  parenthesis  of  the 
numerator  of  (4.31)  has  the  form  of  an  entropy  as¬ 
sociated  with  the  mixing  of  taro  independent  sta¬ 
tistical  systems.  We  can  generalise  Eq.  (4.31) 
trivially  to  include  the  folding  at  three  or  more  re¬ 
gions. 

For  the  folding  at  two  regions,  we  introduce  an 
overlapping  factor  r  by 


It  is  easy  to  see  that  r=0  If  either  p,/p  or  p,/p  is 
zero,  and  r  has  a  maximum,  r  |  =  1 .  at  p,/0 
=  4>,/<t>  =  i.  In  terms  of  r  and  p,  we  hare 

y-n=  j  ix  p(xMx)ln2^/^  dx  p(x) .  (4.33) 

To  compute  m,  (or  p),  we  need  detailed  informa¬ 
tion  about  p,,  p,,  and  p  as  functions  of  x.  How¬ 
ever,  we  can  set  a  rather  good  upper  bound  os  u» 
it  we  know  the  integrals  of  <t>l  and  ps  in  the  over¬ 
lapping  regions.  Let 

M*m  *)•  M*~  (4.34) 


and 


M"  JL  ***^x)'  (4.35) 

we  then  define  a  mean  overlapping  factor  as 


M,  .  M,+M,  M, 

Ra) »^b;h>^rL  v+V 


One  can  show  that 


(4.38) 

(4.37) 


For  the  particular  coupling  leading  to  the  distribu¬ 
tion  funetfowa  of  Figs.  9-11,  we  obtain  m-0.36 
and  the  upper  bound  would  give  n  *  0.43. 

Wb  have  st udled  numerically  the  ««|»is<n*hv 
tutor  sag  ton  | 'ftgmtt  supeueut  for  tht  nsliw  eel 
teastfon  0  defined  rip  Bq.  (3.7)  by  sampling  300  000 
iterations  uadsr  G.  Tbs  averaga  magntflcatfoa 


m  0t  foe  Lyapunov  exponent  pc,  and  foe  average 
overlapping  factor  R0  are 

m0s  1.9873, 

|i0»lnm0=  0.6867, 


Ra  ■  0,9907 , 

respectively.  The  large  magnification  fat*  3)  and 
overlapping  factor  (H  *  1)  are  not  surprising.  Nu¬ 
merically,  foe  universal  function  G(y)  is  vary  close 
to  foe  quadratic  function  /(y)=4y(l  -y).  The  eigen¬ 
function  corresponding  to  the  latter  mapping  is  p 
=  l/[*(l  ~x)fn.  The  iteretioo  of  fly )  gives  rise  to 
the  maximal  magnification  and  overlapping  factor 
m=2,  p»ln3,  and It= 1,  exactly. 

Knowing  foe  Lyapunov  exponent  jx0  tor  G(y),  we 
can  compute  foe  Lyapunov  exponent  p(x)  for  the 
original  /(x)  at  the  2N-N  transition  point  for  large 
N.  la  particular,  for  the  transition  at 

largo  h,  we  have 

p(X.)  =  3-V0.  (4.28) 

Using  the  asymptotic  formula  Eq.  (3.6),  we  obtain 

(4.29) 


(4.30) 


(4.31) 

The  (X  -  Xj  power  dependence  has  been  recastly 
obtained  by  Huberman  and  Rndnick.*  We  are  also 
able  to  determine  the  prefactor  in  addition  to  the 
power  dependence  at  these  transition  points. 

V.  BEHAVIOR  OF  EIGENVALUES 
NEAR  THE  CHAOTIC  TRANSITIONS 

We  learned  from  Sec.  H  that  there  are  two  eigen¬ 
values  q=a  1  with  |q | - 1  in  foe  chaotic  two-region. 
On  the  other  hand,  there  is  only  one  eigenvalue 
q  =  1  with  |q  |  *  1  In  foe  chaotic  one  region.  Time, 
when  foe  system  makes  a  2  to  1  chaotic  transition, 
foe  q=  -1  eigenvalue  becomes  an  eigenvalue  of 

l»»l  cl- 

ln  foe  chaotic  two  region,  foe  eigenfunction  p, 
of  q=  1  Is  real  and  positive.  The  support  of  p0 
contains  two  separate  regions  I  and  □  as  In  Fig.  4. 
The  sigenfunction  p4  associated  with  q=  .1  also 
has  foe  seme  support  We  may  choose  p,  as 


with 


f  =  7-7-  =0.449  806967, 
1210 


Hence,  we  have 

^=0.6867^*-^)'. 
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*1“ 


P0,  la  region  I 
-  <f>0,  in  region  H . 


(5.1) 


However,  as  one  enters  tbs  chaotic  one  region, 
the  previous  two  regions  I  and  n  overlap  (see  Fig. 
5).  To  the  first  approximation,  p0  in  the  overlap¬ 
ping  region  may  be  viewed  as  the  sum  of  p„’s  in  I 
and  n,  while  p,  in  the  overlapping  region  is  the 
difference.  The  interference  of  p’s  in  the  overlap¬ 
ping  region  is  responsible  for  reducing  |q(to  less 
than  1.  To  understand  qualitatively  how  q  behaves 
near  the  transition,  we  consider  the  expression 

/d,d,HU,y)P.(y)(^+^g)’ 

-(2 ♦  an)  (5.2) 


which  implies 

/  dxdyK(x,  y)<f>0(y)( 
1+8*  — 


(5.S) 


for  s  in  the  overlapping  region.  One  can  show  that 
the  only  important  region  of  integration  in  the  Ihs 
of  Bq.  (5.2)  is  for  x  in  the  overlapping  region. 

With  *  in  the  overlapping  region  (OR),  we  may 
have  a  erode  estimate  of  1  +  q  as, 


(5.7) 


In  deriving  (5.7),  we  have  also  approximated 
|  P(x)  (VpeUr)  in  the  denominator  of  Bq.  (5.2)  by 
Po»  Equation  (5.7)  Indicates  that  1+q  la  propor¬ 
tional  to  the  integrated  density  of  p0  in  the  over¬ 
lapping  region.  Note  that  the  sine  of  the  overlap¬ 
ping  region  6  is  proportional  to  6X»  X  -  X,  and  that 
p0(x)  is  dominated  by  square-root  singularities 
l/(  I*  -  x,  |  f1,  in  the  overlapping  region.  Hence, 
we  have 


0,  X< A,. 


\>\ 


(5.8) 


In  (5.2)  and  (5.3),  we  have  made  use  of  the  fact 
that  q  and  p  are  real.  When  X  is  smaller  than  Hie 
transition  eigenvalue  A,,  the  regions  I  and  H  are 
disjoint.  In  this  case,  we  have 

and  hence 

t+q-0,  A<  X,.  (5.5) 

When  X>\  regions  I  and  H  overlap  and  the  left- 
hand  side  (the)  of  Bq.  (5.2)  no  longer  vanishes.  We 
(tad  that 


fU.  is.  G»p*m  aftawtioa  of  X-Xt  nearthe  «  k>  1 


The  constant  of  proportionality  depends  on  the  de¬ 
tailed  structure  of  p0(x)  and  on  the  rate  of  change 
of  the  overlapping  region.  However,  we  know  that 
in  the  vicinity  of  X» X,,  (l+q)P  is  aero  for  x<  x, 
and  increases  linearly  with  X  for  x>  \.  In  Fig. 

15,  we  plot  (1  +  tif  ae  a  function  of  X.  The  transi¬ 
tion  and  the  linear  dependence  in  X  are  Indeed  ver¬ 
ified  numerically. 

We  may  generalise  the  method  to  study  the 
change  of  eigenvalue  q  near  a  2  -1  chaotic  transi¬ 
tion.  In  the  chaotic  three  region,  the  eigenvalues 
with  (n |  =  1  are  l,q»“J(-’  +/S),  andqj.  The 
♦it*!*  ijstste  contain#  thrieeeperate pieces  I,  H, 
and  m  as  shown  in  Fig.  5.  The  %  (or  1?,*)  state  is 
obtained  from  p,  by  assigning  three  different  phases 
l.n»,n?  (l.qf.q,)  to  regions  I,  H,  andHI.  Inttae 
chaotic  con  region,  only  tho  eigenvalue  n  =  l  re- 
malas  to  be  of  maguitude  1.  The  other  two  eigen¬ 
value#  q,  q*  now  here  magaitude  leas  than  one 
(t.a.,  |q|<  1).  We  can  also  study  the  rate  of  |q  | 
approaching  one  near  the  traaeittoa  region  x  =X(. 
In  analogy  to  Bq.  (5.3),  we  have 


(5.9) 

For  x<  x,,  we  are  in  the  chaotic  three  region  and 
both  sides  of  Bq.  (5.9)  raatsh  Meattcaliy.  For  X 
>  X,  we  are  in  the  single  ohaatie  region.  Ae  we 

here  shown  is  Sec.  IT,  the  gepe  between  reglnee 


2  VI** 


fdxdyK(x,y)+,( y) 
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I  and  n  and  regions  H  and  HI  develop  nonvanishing 
4>„.  The  magaliwle  ol  S,  in  the  gaps  is  proportion¬ 
al  to  the  incoming  /tax  (source)  via  the  gate  JOB. 
The  of  the  source  is  given  by 

Scores-  /  dx ♦,(*)«  V^5a  (5.10) 

with  9A»  A  -  \.  Thus,  <*>„<*)  in  the  gap  region  is 
also  proportional  to  /Sx.  Now,  we  look  at  the 
numerator  on  the  right-hand  side  (rhs)  of  Eq.  (5.9). 
The  important  contribution  to  the  numerator  comes 
from  the  integration  region  where  y  is  in  the  gap 
region  DF.  The  requirement  of  K(x,  y)+0  implies 
that  x  is  in  the  gaps  or  in  region  H.  Then,  the 
phase  correlation  factor  is 


We  then  have  the  estimate 

-0(1)  f  dy p0(y)  / fdx+Q(x)<x'/~$\.  (5.11) 

m 

Near  q-%,  we  can  rewrite  the  Hut  at  (9.11)  as 

2-ijJq-wi*=i-oq*+  Ino-til’ 

=  2(1-  |q|)+  |q0  —  q I*  —  (l  —  Inl)3- 

(5.12) 

Ignoring  terms  second  order  In  smallness,  we 
have 

ua) 

as  before.  We  have  studied  (1  -  }ij  |)*  as  a  function 
of  X  numerically  and  indeed  observed  the  above  de¬ 
pendence. 
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APTENDU  A:  NUMERICAL  METHOD 

To  find  Dm  eigenvalues  and  eigenfunctions  of  the 
operator  X  given  in  Eq.  (2.7)  We  need  a  relatively 
straightforward,  hat  possibly  inefficient  method. 
Basis  functions  for  the  spues  were  chosen  to  be 

f  tjAm  | 

\o,  otherwise 


where  the  x.  were  a  prescribed  set  of  points  on 
(0, 1).  K  is  then  expanded  in  this  basis  set.  The 
problem  that  arises  Is  that  a  large  matrix  ts 
needed  (on  the  order  at  1000  x  1000)  for  reasonable 
accuracy.  The  reason  for  this  is  essentially  that 
discussed  in  Sec.  H;  that  is,  the  width  of  K  con¬ 
trols  the  number  of  bifurcations  that  are  possible. 
A  better  choice  of  basis  functions  might  alleviate 
this  problem  somewhat  Fortunately,  most  of  the 
elements  at  K  are  aero,  in  fact,  the  number  of 
nonsero  elements  is  about  twice  the  linear  dimen¬ 
sion  of  K. 

We  are  now  faced  with  finding  the  eigenvalues 
and  eigenvectors  of  a  large  sparse  matrix.  Since 
we  are  only  interested  in  the  largest  eigenvalues 
it  is  appropriate  to  use  an  iterative  technique.  To 
find  the  equilibrium-density  distribution  straight¬ 
forward  iteration  of  the  discretized  version  of  Eq. 
(Al)  works  fairly  well,  the  only  complication  be¬ 
ing  in  sorting  out  <f>0  from  the  other  eigenvectors 
With  eigenvalues  of  modulus  1.  Our  pictures  of  <f>0 
used  this  method. 

The  eigenvalue  (and  associated  eigenvector) 
problem  is  more  complicated.  We  used  a  version 
of  the  Lanczos  method.  See,  for  example,  Wilkin¬ 
son''  for  a  complete  description.  The  idea  is  to 
start  from  an  arbitrary  left  state  L,  aad  an  ar¬ 
bitrary  right  state  Rt  and  Iterate: 

—  fljUj  +  bg  Jtg , 

LgK^agLg  +  SgbgLg.  (A2) 

At  the  nth  step  we  have, 

RR*  -  o»R*+ 1,2,4  +  ®,N2fci , 

1,2*®,!,+  +  *  (A3) 

The  newly  generated  R,*  is  chosen  to  be  orthogonal 
to  all  previous  generated  similarly  tor  L.  We 
choose  (Ln  |A„)  =  1.  This  completely  determines 
the  coefficients  a.  and  S.  is  Juat  a  sign  *  1.  In 
this  scheme  it  is  easily  seen  that  a  tridtagouRl 
matrix  for  K  is  generated.  We  typically  used  a 
30  x  30  matrix.  Actually  it  is  important  to  use 
several  sizes  say  29-30  and  only  the  stable  eigen¬ 
values  are  relevant.  That  is,  new  eigenvalues 
an  introduced  at  each  stage  aad  they  typically 
change  as  the  matrix  si se  is  incresssd.  Tbs  eigen¬ 
vectors  can  easily  be  calculated  this  way,  but  we 
have  not  actually  thaw  that.  Tbs  absolute  accuracy 
of  predicting  transitions  was  of  the  order  of  sever¬ 
al  percent,  tap*  the  relative  accuracy  la  touch  bat¬ 
ter.  This  aaableeua  to  oosmate  the  sealinc  of 
eigenvalues  soar  transition  qpttt#  aesacatety. 


(Al) 
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AffENDIXB:  ALTEKNATTVE  FHOOF  OF  LEMMA  3 

In  thin  appendix  ire  present  an  alternative  proof 
of  Lemma  3  that  is  valid  for  any  number  of  dimen¬ 
sions.  We  start  with  a  proposition  that  ia  proved 
in  Ref.  6.  K  states  that  ander  fairly  general  eon- 
ditions  an  operator  satisfy inc  Bqs.  (2.2)  and  (2.3) 
has  a  cyclic  peripheral  spectrum.  That  is  to  say 
that  if  there  are  N  eigenvalues  with  absolute  value 
one,  those  eigenvalues  are  all  of  the  Mb  roots  of 
unity.  If  there  is  a  degeneracy  the  spectrum  most 
be  divided  into  classes  satisfying  the  above  con¬ 
dition.  We  will  assume  that  there  is  no  degenera¬ 
cy.  Thaw  are  are  given 

•  • » dr-i  > 

each  that 

with 

«*«**'.  (Bl) 

Front  Lemma  2  we  know  that  we  can  choose 

d/*)-«"'wdo(x)»  j«l,...,W-l.  (B2) 

We  fOrther  choose  the  phases  such  that  at  a  par¬ 
ticular  point  *««  with  p((a)#0, 

0,(s)-O,  J*1 . AT- 1.  (B3) 

We  first  prove  a  Lemma:  If,  for  K,  satisfying 
Kqs.  (2.2)  and  (2.3)  there  exist  p0, . . . ,  such 
that 

f/Cix.yWyHy-^Jx)  (B4) 


and 


proved. 

We  now  prove  the  important  theorem  which  says 
that  the  aupport  of  the  eigenfanettoas  conaiet  of 
disjoint  ragiona  and  the  phases  of  sach  eigenfunc¬ 
tion  are  constant  within  those  regions. 

Theorem.  There  exist  N nonnagative  functions 
x/x)  having  disjoint  aupport  such  that 

*X,«XW  (B7) 

and  the  x,  serve  as  a  basis  for  expansion  of  the 
eigenfunctions  having  modulus  1. 

Proof.  Construct 

x*-y  |l  ♦«**  +  ...+  ****-»  Id, , 

TUmJj  ll  +  w**1  +  «*«*»♦...  |#o, 

X»“^ |l +  «*«*** +  »*•**+. . .  Ido, 


+  +  Ido* 

It  is  trivial  to  verify  from  Kqs.  (Bl)  and  (B2)  that 
without  the  absolute  value  signs,  the  x,  satisfy  Eq. 
(B4).  From  the  above  lemma  we  see  that  they 
satisfy  Bq.  (B7).  We  further  note  that  at  x«a 
where  0/a)*0, 

Xo(a)-do^0, 

X<(«)*0,  i*  0. 

(This  follows  from  the  fact  that  the  sum  of  the  Wth 
roots  of  unity  is  aero.)  Not  construct 


drU)*A>(x), 

then 

/jc(x,y)  |d/y)  \dy-  |d*,U)  |  *  (B5) 

proof.  Taking  toe  absolute  value  of  Bq.  (Bd),  we 
obtaia 

/ir(s,y)|di(y)|dy»|dM(*)|.  (B«) 

Integrating  Bq.  (B8)  gives  the  cyclic  art  of  Inequal¬ 
ities 

/|d/*)|d*»/  |d*,(*)|d*. 

Because  of  the  eyHc  mtoare  toe  equality  atga  must 
hoM,  wtoch  implies  Bq.  (Mfc  naeept  pertnqw  on 


to”(Xo+Xa+*.*+XsJ, 

«i*(xo+5Xt+J3rXb  +  .*.),  (BB) 

-SL—  »— n. a ,  .at, 

otCe  inwB  ranenona  iiimy 

so  we  conclude  tost 

Mow  at  * ■«  wa  have 

t,“d».  1,  •  •  •  -*  1  • 

Thus 

C,~l,  jrn  o . Jf~J. 
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as 

Taking  Ow  absolute  value  at  Bqp.  (B8)  we  obtain 

|tj|i 

Because  the  coefficients  v>  point  in  different  dir¬ 


ections  in  the  camples  plane,  the  functions  x,  moat 
have  disjoint  support.  Finally  the  construction  of 
demonstrate  the  constancy  of  their  phases  (or'*) 
in  the  Mi  region. 
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Thi»  paper  Xiidice  the  iterative  properties  of  the  two-penmeter  faitily  of  nape  x. t ,  »  1  +  ax‘  +  bx'.  These  mips 
can  have  either  one  or  three  extrema.  Multipie  extrema  lead  to  a  oompiex  region  of  iterative  stability  and  to  new 
typee  of  iterative  behavior.  In  particular,  the  Ftipmhaum  critical  Unc  (locue  of  2*  -cycle  aoeumulatiooe)  preeent  for 
bs 0  terminates  at  two  tricritial  paints  T  and  T  with  ronartinatee  (a» -(0,  -  1.59490)  and  ( -  2.81403, 
1.40701),  respectively.  There  is  also  a  “dual"  critical  line  terminating  in  two  additional  tricritical  points 
fm ,  ( —  3.119 10, 2.543  71)  and  f’  —  (0.955  il,  -  1.149  81).  Behavior  near  the  tricritical  points  is  controlled  by  a 
new  Used  point  fT  at  the  fhnctkmal  recursion  relatioa/, +,fx)  -  [l//,(l)y°V/.(l))-  This  (bed  point  has  two  relevant 
directions  and,  therefore,  behavior  near  a  tricritical  point  depends  on  two  new  universal  numbers  6ln  «=  7.284  49 
and  <5^  —  2.857  13.  Crossover  and  scaling  behavior  near  the  tricritical  points  are  explored. 


L  INTRODUCTION 

Iterative  properties  *„♦,  =  /(*„)  ot  one-dimension* 
si  continuous  maps  provide  What  is  probably  the 
simplest  examples  at  nonlinear  dissipative  dynam¬ 
ical  behavior  and  have  attracted  considerable  at¬ 
tention  in  recent  years.1'*  Changing  a  single  con¬ 
trol  parameter  in  such  maps  may  lead  from  stable 
fixed-point  behavior  via  a  sequence  of  bifurcations 
to  chaotic  or  quasicbaotic  behavior.  Feigenbaum* 
showed  that  the  limiting  form  of  such  bifurcation 
sequences  exhibited  scaling  properties  character¬ 
ised  by  universal  numbers,  in  strict  analogy  to 
critical  behavior  In  statistical  mechanics,  and  in¬ 
deed,  that  the  renormalisation-group  ideas  devel¬ 
oped  In  the  statistical-mechanics  context4  could  be 
fruitfully  carried  over.  Such  trifurcation  sequences 
are  also  seen  in  higher -dimensional  models1'*'* 
and  in  experiments.*  They  represent  one  possible 
kind  of  “transition  to  chaos."  Feigenbaum’s  uni¬ 
versal  numbers  have,  in  fact,  cropped  up  in  Borne 
higher -dimensional  systems.7 

The  prototypical  one-dimensional  map  is  the  so- 
called  logistic  map, 

/Cx)-Ax(l-x),  (l.l) 

which  for  0  <  x  <  4  maps  the  unit  interval  [0, 1]  onto 
itself.  Equation  (1.1)  is  equivalent  under  change  of 
origin  and  rescaling  to  the  symmetric  quadratic 

map,* 

fW-l+ex*.  (1.2) 

In  particular,  the  main  bifurcation  sequence  has  Its 
limit  at  a,- -1.401 18.  There  is1  a  simple  sequence 
of  2"  cycles  for  a,  <  a  <  i  and  a  corresponding  in¬ 
verse-bifurcation  sequence  of  more  complicated 
2*  bands*  for  -2  <a<  a„.  The  vicinity  of  the  point 
a,  is  characterised  by  a  sealing  behavior.*  For  ex¬ 
ample,  the  position  a,  of  the  superstate  2*  cycle 
lor  a>a,  and  large  a  varies  as 
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0,-0,=  A/ 6},  (1.3) 

where  8, »  4.689 20  Is  one  of  Feigenbaum’s  univer¬ 
sal  “critial  exponents”  and  the  nonuniversal  ampli¬ 
tude  A  =  9. 349.  (The  actual  exponents  In  statistical 
mechanics  are  really  related  to  Ins,.)  So  far,  we 
have  discussed  only  the  logistic  map  (1.2).  It  is 
known,  however,  that  the  behavior  of  the  logistic 
map  provides  a  precise  model  for  a  large  class10 
of  single-peak  functions:  The  ordering  of  a  cycles11 
and  the  critical  behavior*’1*  are  universal;  only 
the  nonuniversal  constants  a,  and  A  vary  from 
map  to  map. 

D  is  the  purpose  of  this  paper1*  to  study  the  be¬ 
havior  of  a  simple  two-parameter  class  of  maps 
which  includes  (1.2),  the  symmetric  quartic  maps, 

/(x)*l+«f’+6x4.  d.4) 

The  distinctive  feature  of  the  maps  (1.4)  is  that 
for  ad  <  0  they  have  three  extrema,  at 

x„=0,  x,  =  *  V  —a/ib ,  (1.5) 

rather  than  one  as  (1.2).  This  multiple-peak  or 
multiple-valley  structure  makes  possible  a  pletho¬ 
ra  of  now,  interesting,  and  complicated  iterative 
behavior.  It  is  the  purpose  at  this  paper  to  sur¬ 
vey  in  a  partially  empirical  and  certainly  nonrig¬ 
orous  way  some  of  the  new  behavior  in  the  context 
of  the  two-parameter  family  (1.4). 

Section  n  of  the  paper  treats  general,  noncrltlcal 
properties.  In  the  one -extremum  regions  (oft  >  0) 
behavior  like  that  of  the  logistic  map  is  observed, 
only  the  accumulation  point  a,  is  drawn  out  into  a 
“critical  line.”  This  critical  line  terminates  in 
the  regions  ab  <0,  where  it  is  "pinched  off”  at 
tricritical  points  t  and  T  by  the  existence  of  be¬ 
havior  mere  complicated  than  a  monotonic  cas¬ 
cade  of  2" -cycle  regions.  When  mb  <  0,  it  la  quite 
possible  for  two  stable  cycles  (but  apparently  no 
more  than  two)  to  coexist.  Super  stable  cycles 
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moat  pass  through  aw  of  th*  atiuw  (l.SJu  For 
special  nloei  of  (0,6)  a  single  "doubly  stable" 
cycle  may  use  two  of  them.  Trleritlcal  points  are 
associated  with  x-  limits  of  doubly  stable  I* 
cycles. 

We  develop  a  special  term  at  topological  coctju*  . 
gacy,  called  duality.  Dual  points  P  and  P  have 
the  same  iterative  properties.  The  dual  of  the 
critical  line  TT1  is  another  critical  line  ff*  is 
the  a,A  plane. 

We  determine  the  region  of  the  a,b  plane  over 
which  the  function  (1.4)  is  Iteratively  stable,  Le„ 
possesses  a  set  of  points  *  of  aonsero  measure 
with  bounded  itineraries. 

Section  in  is  devoted  to  the  special  properties 
of  the  map  (1.4)  near  the  endpoints  T  and  V  of 
the  critical  line.  Feigenbaum'  showed  that  the 
points  on  the  critical  Lae  were  aseovlated  with 
a  fixed  point  of  a  certain  functional  recursion  re¬ 
lation.  The  trleritlcal  pointer  *nd  V  **•  Associ¬ 
ated  with  a  different  fixed  point  of  the  same  func¬ 
tional  recursion  relation.  At  this  new  fixed  point 
there  are  two  relevant  directions  and,  therefor#, 
two  new  universal  numbers. 

We  discuss  the  scaling  properties  of  the  regions 
near  the  trlcrtticai  points  and  the  crossover  be¬ 
havior  between  critical  and  trleritlcal  regimes. 
The  fixed  points  associated  with  rtf'  are  different 
from  those  associated  with  Tf  tad  related  to  them 
by  topological  eonjugacy. 

n.  GENERAL  PROPERTIES  UNDER  ITERATION 
OF  THE  SYMMETRIC  QUUmC  IMP 

This  section  describes  certain  general,  noncri- 
tical  properties  of  the  map  (1.4). 

A.  HdtehMriw 

The  possible  shspss  of  the  symmetric  quartlc 
map  (1.4)  depend  on  ten  relative  signs  and  magni- 
todss  of  tbs  paramatsrs  #  and  6,  as  shown  In  Fig. 
1.  The  curves  have  ana  local  sstrsmam  la  quad¬ 
rants  I  and  in  and  throe  in  quadrants  n  and  IV. 

The  number  of  intersections  of  ymfU)  with  the 
straight  line  y  -x  may  be  aero,  two,  or  fear.  The 
regions  of  each  behavior  are  Indicated  is  Fig.  1 
and  the  corresponding  curves  are  sketched  la  Fig. 

a. 

Rian* tries  xsx1,x,(*), *,(*), ...  of  an  initial 
point »  under  iteration  (1.4)  are  of  several  types.1 
WO  refer  to  as  t  oarndxf  those  fttetmtes  such  that 
Btesraries  sush  that  |xJ-« 
as  a-  -  are  mbamtkd.  Thu  set  of  x  having 
heundsd  mnsrartus  totes  jtyvhhsw  dssisteof/. 

ours,  than  f  Is  ttarmttaaly  stable,  otherwise,  it 
is  mrwUvah  amhtUt.  We  wish  to  explore  ia  the 
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FK3. 1.  Four  quadrants  Of  ths  t,b  plane,  indicating 
the  number  of  intersections  of  y-l+«x*+*x*  with  j-x. 
Detail  sheys  region  spar  tee  origin.  Nets  that  there  is 
s  small  region  of  quadrant  I  where  there  ere  two  inter- 
seeOosi  and  an  eveti  Smaller  region  of  quadrant  IV 
where  there  are  four.  Corresponding  curves  are 
sketched  ta  ng.  2. 


remainder  of  this  section  the  character  and  extent 
of  the  iterative  domain  of  the  map  (1.4)  for  various 
values  of  the  parameters  0  and  b. 

We  distinguish  the  following  types  at  bounded 
itineraries:  (a)  fixed  poiuts  and  other  s  cycles, 

(b)  itineraries  which,  as  x-«,  approach  (fere  at¬ 
tracted  by)  a  fixed  point  or  other  w  cycle,  and  (c) 
other  bounded  itineraries,  which  we  call  chaotic.14 
The  set  of  x  cycles  corresponds  to  the  set  of  nota¬ 
tions  of 

x -/'»>(*),  (2.1) 

where  we  hove  denoted /<">■/•/#/••••  e/  and 
/•*<*)■/<*<*)>•  This  set  is  countable  and,  there¬ 
fore,  contributes  in  and  of  itself  aero  measure  to 
the  iterative  domain  of  f.  The  x  cycle  xx,xt,  •  •  • , 
x.  ,xn+i  wx,  Is  stable ,  if  a  finite  region  of  nearby 
x’s  is  attracted  to  It,  otherwise,  it  is  xxstebfs. 
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Hi  particular,  it  it  stable  (unstable)  l 

(2-2) 

WkM 

*«,  <2S) 

tb«  n  cycle  to  super  stable }*  A  function  f  to  itera¬ 
tively  stable  if  and  only  if  it  has  one  nr  more 
stable  h  cycles  and/or  a  domain  of  x’a  at  nonsero 
measure  with  chaotic  itineraries. 

When/  to  Iteratively  stable,  sre  shall  wish  to 
.MsHngmsh  the  somber  of  different  types  of 
(hounded)  limiting  behavior  occurring  in  the 
iterative  domain  of  /.  If  all  x's  except  for  a  set 
at  measure  sero  have  the  same  limiting  behavior 
(e.g.,  if  all  are  attracted  to  a  single  stable  n  cycle 
or  if  all  have  chaotic  itineraries  using  the  same 
domain),  then /  to  unimodal if  two  distinct  types 
of  limiting  behavior  occur  with  nonsero  measure, 
then/  to  bimodal,  and  so  forth. 

B.  Tapotofkal  coe)ssKy  duHty 

U  is  useful  before  proceeding  further  to  establish 
two  important  symmetry  properties,  let  A  be  an 
arbitrary  locally  invertible  map  and  define 

geA‘‘e/»A,  <>-4) 

so  /  *  A  egoA"1  Maps /  and#  related  by  (2.4)  are 
called  topolog  'rlly  conjugate.1'*  Itineraries  of 
/  and  g  are  mated  to  one  another  in  a  simple  ooe- 
to-cne  manner;  any  bounded  (unbounded)  itinerary 
of  /  maps  into  a  bounded  (unbounded)  Itinerary  of 
g  and  conversely.  Map  g  is  iteratively  stable  if 
and  only  if  /  to.  Similarly,  n  cycles  map  to  s  cy¬ 
cles,  chaotic  itineraries  to  chaotic  itineraries, 
etc. 

For  general  A  the  functional  forms  of  the  topolog¬ 
ically  conjugate  maps  f  and  g  will  be  quite  differ¬ 
ent.  Two  instances  when  this  to  not  so  are  trans¬ 
lations, 

kfkt,x)mx -a,  A?(o,x)*x+e,  (2.8) 

and  seals  changes, 

kt<o,x)  max,  *il(n,x)*x/a ,  (2.8) 

both  of  which  tok»  nth  degree  polynomials  into  nth 

degree  polynomials.  For  example,  toe  general 
Symmetric  quartic  polynomial  ip  #0) 

F(x>-C«An**BF*  U-7) 

is  tognhigtoslly  conjugate  under  scale  change  to 
the  standard,  normalised  Arm  (1.4,  with  «« AC, 
h-JC’ne 

/•AftP)  eFehf  (C) ,  (*-•) 
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since 

F(*)'c[l  +AC(!)’+BC’(§)‘}  (2.9) 

There  ia  a  special  form  at  topological  conjugacy 
which  is  exceedingly  valuable  in  studying  the  itera¬ 
tive  properties  of  the  quartic  map  (1.4).  Note  that 

fix) =  1  -  ^  +  *’)“ =Gl(G,(*)) ,  (2.10) 

with 

G^Hl-^+b*2 

and 

G.W^+x’.  (2.11) 

Because  Gj*G,“GIl*  (G,*G,)  •  G„  /  to  topologi¬ 
cally  conjugate  to  the  map  Gt-G ,  and,  therefore, 
through  (2.8)  to  the  normalized  map 

/=*;*(</)•  G,* 0,0 A, (rf),  d^-~r  (l -g y . 

(2.12) 

The  map/  is  called  the  dual  at  f  and  has  the  farm 
/(x)  =  l+fix2+bx4, 
with 

a  =  (4b  -«*)(8ob  +  (4b  -a2)2]  /82b2, 

b-{8ob+(4b-c.2)2lV(8b)4.  (1,13) 

DuaUty  is  reciprocal,  (7)=/.  AU  of  the  topologi¬ 
cal  correspondences  mentioned  after  (2.4)  hold 
between/  and  its  dual.  For  example,  if  /  has  an 
n  cycle  using  the  point  x,  then  f  has  a  correspond¬ 
ing  n  cycle  using  the  point 

x»*-‘(d).CaW  =  (^+x2)/d.  (2.14) 

m  particular,  it  follows  from  (2.14)  that,  if/has 
an  n  cycle  wring  a  side  extremum  (*, =±  V-o/Sfc, 
ob  <  0),  then  the  corresponding  n  cycle  of  /  uses 
the  central  extremum  (f,  *0)  and,  conversely,17 
if  the  n  cycle  of  /  uses  the  central  extremum  xo=0, 
then  that  of  /  ueea  one  of  the  side  extrema  [x0 
-±  (-d/2S)1/2J. 

The  duality  relations  (2.13)  will  be  useful  in 
Sees.  nC  and  HD  below.  The  points  on  the  curve 

80b +(4 b-«*)*»0  (2.15) 

are  all  to*  dual  to  the  origin.  Memalnlng  a  axle 
points  0*0,  « #0)  map  to  infinity  under  duality, 
whU*  remaining  b-axto  potato  (s'O.b  *0)  map  to 
the  pnrnbola  b  *•  2/4.  in  addition,  there  to  a  lino 
of  self -duality,  (a,5)-fc,b),  whtoh  eaa  bn  parmet- 
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remark  appUee  toD  ends  and  their  duals  D  and 
B.  This  occurs  because  the  origin  is  the  dual 
image  of  die  entire  line  (3.16). 

The  iterative  behavior  of  /  is  completely  known 
along  the  self -dual  line,  since  g(x)  in  (2.17)  is  the 
well-studied1  quadratic  map  (1.2).  The  map/  is 
iteratively  stable  for  all  -2  <c  <i  and  otherwise 
unstable.  At  c  values  for  which  g  has  a  stable  2* 
cycle,  xux„ .... /  has  two  stable  s  cycles 
xux„ . . .  and  *„*«, . ...  and  is,  therefore,  M- 
modaL  Ate  values  for  which  ghas  a  stable 
(2s  +  l)-cycle,  /  has  a  stable  (3*+l)  cycle  and 
is  unimodaL  The  main  Ufurcatton  sequence  of 
g  ends  at1  c  '  -1.401 16. . . ,  which  translates  into 
to, »)•- (-1.67813,  0.17766). 


FKJ.  3.  Mapping  of  apecial  lines  under  Aullty  (2.13). 
The  origin  maps  to  the  line  (2.15).  The  i  axis  maps  to 
the  parabola  4-aV*.  Both  the  self-dual  line  and  the 
inverse  origin  have  breaches  In  Quadrant  IV,  aa  ahown 
In  the  detail. 


rised  ss 

«  =  2c,(l+c),  6=c*(l+c)s,  (2.16) 

corresponding  to 

/=*;1(l+c)ag<«e*,<l+C>,  (2.17) 

With 

gi jc)  =  1+cx*. 

These  special  lines  are  shown  in  Fig.  3.  The 
mapping  under  duality  of  some  sample  regions  of 
the  a,6  plane  is  tftown  in  Fig.  4.  Notice  that  re¬ 
gions  A  and  C  of  Fig.  4  are  adjacent,  while  their 
duals  A  and  6  meet  only  at  the  origin.  A  similar 


C.  Stable,  sepentabis.  md  doebtjr  (table  a  cycica 

(Ref.  15) 

S  is  straightforward  to  solve  the  a -cycle  condi¬ 
tion  (2.1)  numerically  and  then  to  check  stability 
(2.2).  The  condition  (2.2)  for  a  eeperstable  a  cycle 
requires  that  one  or  more  of  the  potato  xltx„ . . . , 
x,  of  the  cycle  be  an  extremum  of  /.  The  central 
extremum  *o*0  may  always  appear;  however,  at 
most  one  of  the  two  side  extrema  can  appear  In 
any  given  cycle,  since  the  evenness  of  /  guaran- 
teoe  that  the  Itineraries  of  xt*±  V-a/tt  (s6<0) 
coincide  alter  one  step.  Thus,  ssch  super  stable 
cycle  can  be  classified  4,  -,  0,  0*,  or  0-,  ac¬ 
cording  as  it  uses  the  indies  tod  extrema.  The 
last  two  cases  denote  a  cycles  (a*  2)  which  use 
both  the  central  extremum  and  one  aide  extremum 
and  are,  therefore,  dombly  stable.  Doubly  stable 
cycles  am,  indeed,  particularly  stable.  When 
(1.4)  has  mors  than  one  extremum,  all  are  quad¬ 
ratic.  Thas,  if  x  to  a  member  of  a  doubly  stable 
m  cycle, 

/<•>(*+ to) «*40(«**),  a-4  (3.18) 

Instead  <f»*l  (ordinary  eoper stable  cycle).  The 
argument  following  (3.14)  dhows  that  the  dual  of  a 
doubly  stable  cycle  Is  also  doubly  stable,  provided 
only  that  tbs  dual  points  50  and  x,  are  distinct.  The 
exceptional  esse,  i0«x,“  0,  occurs  only  for  6*0, 
i.e.,  for 

/(x)  *  1 +6x4 ,  (2.19) 

in  which  cap*/  has  s  single  fturtie  extremum. 

Any  cycle  tatted  uses  this  quartic  extremum  sat¬ 
isfies  (2.16)  and  will  also  bt  called  doubly  stable. 

Suppose  that  for  a  particular  sit'  of  parameters 
fc„6  j/  has  a  doubly  stable  cycle.  K  to  clear  that 
Itineraries  bejongjitf  to  all  points*  sufficiently 
close  to  | ay  pf  tbs  extrema  x,  dm  «*traetod 

to  suck  a  cycle.  It  turns  out  smpiriealfe  (see 
See.  DD)  that/ to  always  aimodaL  If /la  doubly 
stable  for  («W**)>  then/  is  aalmndal  for  (s,6)  la 
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some  nt  ox-hood  of  (a,,  b0).  Doubly  stable  points 
will  play  an  important  role  in  the  trlcritical  be¬ 
havior  which  terminates  the  Feigenbaum  critical 
line. 

We  illustrate  some  of  these  notations  by  simple 
examples.  Regions  of  stable  and  superstable  fixed 
points  (1  cycles)  are  displayed  in  Tigs.  5  and  6. 
Since /( 0)-l  by  our  choice  of  normalization  (1.4), 
superstable  fixed  points  can  only  occur  for  (a, 6) 
values  such  that 

*t(a,6)  ■*,=/(*♦)  OT  *-(«,&)■*-  =/(*-),  (2.20) 

and  the  various  branches  of  the  curve  in  Fig.  S 
have  been  labeled  accordingly.  The  condition 
(2.20)  reduces  to  (2.1S).  This  is  not  an  accident: 
Bad  we  chosen  the  more  general  normalization 
(2.1),  the  condition  for  a  super  stable  fixed  point 
associated  with  the  central  peak  would  have  been 
C  =0,  which  via  (2.8)  is  topologically  conjugate  to 
=  0.  The  discussion  following  (2.14)  then  shows 
that  the  xt  superstable  fixed  points  must  occur  for 
ia,b)  dual  to  (0,0). 

The  regions  of  fixed-point  stability  shown  in 
Figs.  5  and  8  contain  the  superstable  line  and  are 
bounded  [see  (2.2)]  by  curves  determined  from 

Vjfs.h)  »*,*=/(*»), 

with  (2.21) 


Over  most  of  the  1-cycle  region,  /  has  only  one 
stable  fixed  point1*;  however,  there  is  a  narrow 
region  visible  in  Fig.  8  in  which  there  are  two  in¬ 
dependent  stable  fixed  points.  A  typical  example 
of  this  behavior  is  shown  in  Fig.  7:  The  domains 
of  attraction  of  the  two  stable  fixed  points  are  dis¬ 
joint  and,  between  them,  exhaust  the  Iterative  do- 


FlQ.  3.  Loot  of  stable  and  superstable  fixed  pointe  in 
tee, 3  plane.  The  superstore leone  coincides  with  the 
dual-erifto  earn  of  Fig.  S  end  its  brenehee  have  bees 
labeled  at.  or  s.  anoosteg  to  te  aeeoeteted  extremum. 
The  »iree«  lute)  «sgte  of  fixed- point  ateNilty  to- 
ehtrtee  te  auperstebte  loose.  The  upper  edge  of  the 

tangmey  ecodjtfane 


FIG.  6.  Detail  of  the  stable  fixed-point  region  shown 
in  Fig.  5.  In  the  very  narrow  region  of  intersection  of 
the  two  branches  of  stability  the  function /(x)  has  two 
Independent  stable  fixed  points. 


of  /  (exeept,  of  course,  for  a  set  of  measure 
xero,  die  unstable  fixed  points).  The  map  /  is, 
therefore,  bimodal. 

The  superstable  2  cycles  are  shown  in  Fig.  8. 

The  xa  super  stable  2  cycle  satisfies 

0=xo*f(a,(xJ  =/(,,(0)  =  1  +u  +6  .  (2.22) 

The  xt  superstable  2  cycles  can  be  found  in  analogy 
with  (2.20);  however,  it  ia  much  easier  to  invoke 
duality  and  simply  compute  the  dual  of  (2.22). 

Intersections  of  superstable  »  cycle  curves  (e.g.. 
Fig.  8)  require  comment.  The  set  of  curves  be¬ 
longing  to  a  single  extremum  cannot  Intersect  ex¬ 
cept  at  the  origin.  Similarly,  curves  belonging  to 
xt  cannot  intersect  those  belonging  to  x.,  since  x, 
and  x.  have  the  same  itinerary  after  one  step.  x0 
curves  can  intersect  x,  curves,  as  illustrated  by 
Fig.  8.  These  intersections  are  of  two  quite  dlf- 


FIG.  7.  Iterative  plot  of  the  function /(*)«  1  -  0.76xs 
+  0.03 x*.  Thia  function  tea  two  (table  fixed  pointe  at 
xA  »  0,070  and  **-  -3.603.  The  Iterative  domain  of  /  is 
|*|  <6.470.  The  domain*  of  attraction  of  xA  and  x,  are 
Indicated  below  the  graph.  They  exteuat  the  iterative 
dnmeia  of/  (caoept  for  a  set  of  measure  aero)  and  inter¬ 
lace  with  each  other  as  shows. 


tote 
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FU.  8.  ShpsretsbU  1  cycles.  The  osutrsl-extramum 
(*•)  curve  ie  Just  the  lines +4+1-0.  Ksch  branch  le 
labeled  according  to  the  peak  frotnwhloh  its  stability 
derives.  The  point  A  and  its  dual  A  are  doubly  stable. 
Cither  intersections  are  nondynamic al. 


ferent  types.  At  the  intersection  point  there  may 
be  two  Independent  super  stable  cycles,  each  with 
its  own  domain  a I  attraction.  At  such  points  /  is 
bi modal  and  the  curves  are  dynamically  indepen¬ 
dent.  Alternatively,  there  may  be  a  single  doubly 
stable  cycle  at  the  Intersection  point,  hi  this  case 
/  is  unimodal  and  the  curves  are  dynamically 
coupled.  InFig.  8  only  the  intersection  A  and  its 
dual  image  A  are  dynamical.  The  remaining  in¬ 
tersections  are  nondynamic aL 
Figure  9  shows  the  regions  of  stable  2-cycle 
behavior.  Note  the  bimodal  rsgir  ns  surrounding 
the  nondynamlcal  intersections  at  Fig.  8.  The 
region  of  1 -cycle  stability  shown  in  Fig.  S  nests 
snugly  against  the  region  of  2-cycle  stability,  as 
shown  in  Fig.  10.  In  particular,  following  the 


plan*. 


b 

^  .  a 

>  • 

3"  i  -i  •-I — i 

TO.  10.  Nesting  regions  of  1- and  8-eyols  stability. 
Where  suoh  bifurcations  proceed  to  completion,  they 
end  in  a  Pelgeobaum  critical  line.  The  main  critical 
Una  TV  and  its  deal  ff*  are  shown.  The  normal  bifur¬ 
cation  sequence  does  not  occur  beyond  these  trioriUoal 
points. 

negative  a  axle  away  from  the  origin,  one  sees 
Hie  beginning  of  the  main  bifurcation  sequence  of 
the  quadratic  map  (1.2).  Such  sequencing  is  also 
seen  in  other  regions  of  Fig.  10.  Whenever  such 
bifurcations  run  to  completion,  they  terminate 
in  a  Feigenbaum  critical  point.  Regions  at  nesting 
therefore  indicate  a  Feigenbaum  critical  line. 

Note,  however,  that  the  complicated  topology  of 
the  2-cycle  regions  already  requires  that  the  cri¬ 
tical  line  be  broken  into  segments.  We  have  cal¬ 
culated  the  position  of  the  main1*  critical  lines. 
They  are  included  for  convenience  in  Fig.  10  and 
will  be  discussed  in  Sec.  HL 
Figures  11  and  12  show  loci  of  superstable  4 
cycles  and  3  cycles.  Dynamical  intersections 
(which  will  play  a  role  in  Sec.  in)  are  indicated. 

The  general  picture  which  emerges  far  the  central- 
extremum  (**)  curves  is  one  or  more  diagonal 
lines  running  across  the  (a,b)  plane  from  upper 
left  to  lower  right  and  intersecting  the  negative  a 
axis  at  the  known  values  for  the  simple  quadratic 
map  (1.2)  and  in  the  standard  sequence.*0  (Note 
that  there  are  two  such  4  cycles,  one  associated 
with  bifurcation  of  the  2  cycle  and  the  other  in  the 
region  of  Inverse  bifurcations  on  the  far  side  of 
the  critical  line.)  In  addition,  there  are  a  large 
number  of  “finger  F’  in  quadrants  □  and  IV,  some 
very  narrow.  These  poke  in  towards  the  origin, 
crons  no  axes,  and  go  off  to  tsfinlly,  apparently 
along  the  lines  t* 9,  e+MO,  e+h+1-0,  and 
e  +6 + 2  *  0.  The  x,  earves  are  obtained  from  the 
central  sartmsum  carves  by  duality  (2.12).  Wher¬ 
ever  an  xt  carve  orsoseo  (2.1$  (nee  Fig.  I),  the 
correopowllHi  dual  poaees  throagh  the  origin, 
creating  the  loops  visible  in  quadrant  IV. 
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FK5. 11.  Superstable  4  cycles.  The  4  cycles  necessary 
for  the  blfurcattcn  sequences  leading  to  XT'  end  tf1 
(Fig.  10)  la  present.  Much  adc&tlonal  structure  Is  also 
present.  Several  of  the  4-cycle  lines  are  too  narrowly 
spaced  to  be  resolvable  on  this  graph. 


FIG.  13.  Illustration  of  the  bin  construction.  A  big 
box  sad  two  possible  small  boxes  are  shown.  If  the 
function  remains  inside  spy  one  (or  more)  of  the  boxes 
over  the  appropriate  domain,  then  it  is  guaranteed  to 
be  Iteratively  stable. 


D.  Iterative  stabgity 

We  now  return  to  the  question  posed  in  Sec.  HA 
of  determining  the  region  of  the  ta,b)  plane  over 
which  the  function  (1.4)  is  iteratively  stable.  !*. 
principle,  this  is  Just  the  union  over  it  of  tl  >«•-- 
gions  of  stability  (see,  e.g.,  Figs.  B  and  9,  all 
h  cycles  plus,  in  addition,  those  regions  where 
/  has  no  stable  it  cycles  but  is  iteratively  stable 
by  virtue  of  chaotic  itineraries.*1 


IT*  sad  ft*.  Otherwise,  3-cycle  structure  is  hrcetRy 
startler  te  2-  sad  4-oyeIe  etmetoms.  Several  of  the 
3-eye  le  Uses  am  toe  aarrwwly  epeaed  to  he  resolvable 
ea  tUs  graph. 


Maps  of  the  types  4  and  □„  (see  Figs.  1  and  2), 
for  which  there  are  no  solutions  of  x=/(x),  have 
unbounded  itineraries  for  all  x.  Their  iterative 
domain  in  the  empty  set  and  they  are  iteratively 
unstable.  At  the  boundary  between  (1,-i-ilg)  and 
(It+llJ,  >=/(*)  is  tangent  to  y =x,  which  is  Just 
the  condition  (2.21)  giving  the  upper  edge  of  the 
stable  1-cycle  region.  Just  beyond  this  boundary 
in  (1,+iy  there  are  two  fixed  points,  x„v  (see  Fig. 
2)  and  x,,*,  with /'OO >1  and/' (x„J<l.  B 
/'(xmi„)>-l  (as  is  true  throughout  I,),  then  any 
x  with  |x |  <  flows  to  xm,„  and  any  x  with  |x| 
>xma  has  an  unbounded  itinerary,  so  /  is  itera¬ 
tively  stable  and  unimodaL 
In  other  regions  of  the  ( a,b )  plane  the  situation 
is  more  complicated  and  we  must  develop  a  meth¬ 
odology  for  establishing  iterative  stability.  Figure 
13  illustrates  one  such  criterion,  which  we  call 
the  box  construction.  Let  Xmn5x„,  (a, 6)  denote 
the  solution  of  x  =  fix)  at  largest  absolute  value. 
Suppose 

l*nml*l/(*)|  V  (2.23) 

We  refer  to  this  as  the  "big-boa”  condition.  It  is 
easy  to  see  that  points  with  |  x  |  « I  Xm«|  have 
bounded  itineraries,  while  those  with  |x|>| x„.J 
have  unbounded  itineraries.  Similar  “small -box" 
constructions  can  be  done  separately  for  centeral 
and  isle  extrema,  aa  suggested  in  Fig.  13.  B  any 
or  all  of  die  extrema  of  /  are  tumble  in  the  above 
sense,  then /  is  iteratively  stable.  For  example, 
in  quadrant  m,  where/  has  a  single  peak  at  x,-0. 
the  box  condition  (2.23)  requires 

|Xnn|  */<0)=l, 


(2.34) 
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so  the  boundary  of  guaranteed  stability  is  the  line 
-l=*n»=\f(x«J=/(l)  =  l+fl+&, 
and  (2.24)  is  equivalent  to 

a*b>~  2.  (2.2S) 

When  two  peaks  are  independently  (small-)  box- 
able,  then  /  is  at  least  bimodal.  Systematic  ap¬ 
plication  of  the  box  conditions  to  all  peaks  and 
over  the  whole  (a,b)  plane  leads  to  a  region  S  of 
guaranteed  iterative  stability,  pictured  in  Fig.  14. 
This  rather  complicated  figure  goes  into  itself 
under  duality. 

Failure  of  the  box  conditions  may  or  may  not 
imply  iterative  instability.  For  example,  com¬ 
parison  of  Figs.  9  and  14  already  shows  that  the 
region  of  stable  2  cycles  extends  beyond  the  bound¬ 
ary  of  guaranteed  stability  into  regions  D,  and  D, 
of  quadrants  II  and  IV,  respectively.  We  develop 
this  situation  in  more  detail.  The  map  (1.4)  is 
infinitely  differentiable  and  it  is  easy  lj  verify 
that  its  Schwartzian  derivative 


erary  is  unbounded),  /  has  no  stable  a  cycles,  ta 
the  former  case/  is  iteratively  stable  and  uni- 
modal.”  in  die  latter  case  /  is  Iteratively  un¬ 
stable.**  Thus,  the  region  a  +  6<- 2  of  quadrant 
m  and  its  dual  image  in  quadrant  n  are  iteratively 
unstable. 

In  quadrants  Q  and  IV,  /  has  three  extrema  and 
does  not  satisfy  the  hypotheses  of  Singer’s  theo¬ 
rem.  We  have  proceeded  empirically.  The  result 
of  our  study  may  be  summarised  in  a  form  that 
applies  over  the  entire  (a,  6)  plane  by  the  following 
three  statements. 

(1)  When  all  extrema  have  unbounded  itineraries, 
then  f  is  iterativety  unstable. 

(2)  When  xo=0  has  a  bounded  itinerary  and  xt 
have  an  unbounded  itinerary  or  vice  versa,  then 
/  is  iteratively  stable  and  unimodaL 

(3)  When  both  x0  and  xt  have  bounded  itineraries, 
then  f  is  iteratively  stable.  If  the  two  itineraries 
are  asymptotically  distinct,  then  /  is  bimodal;  if 
they  are  asymptotically  convergent,  then/  is  uni- 
modal. 


S/(x)« 


/~(«)  a/ru)V 

/'(*)  ~2\ /'{*)/ 


(2.26) 


is  everywhere  negative.  In  quadrant  m,  where  f 
is  single  peaked,  it  satisfies  the  hypotheses  of  a 
theorem  due  to  Singer,****  which  concludes  that 
any  stable  *  cycle  attracts  the  peak.  Thus,  when 
(2.25)  is  satisfied  (l.e.,  when  the  peak  itinerary  is 
bounded),  /  has  at  most  a  single  stable  n  cycle, 
while,  when  (2.25)  fails  (i  when  the  peak  itln- 


•W.  14.  Rerettve  •Ubillty  via  0»  bos  ocaetnotton. 
the  oross-httahed  isgtaa  f  is  guaranteed  stable,  the 
refloat  V  and  their  deals  0  are  feeruaeed  wuAh. 
Stability  of  the  reflet  D  and  their  dtle  D  oatot  be 
deelded  by  the  rtmple  bos  eoestruettoe. 


B  all  extrema  lie  outside  the  big  box  of  Fig.  13, 
then  statement  (1)  applies  and  /  is  Iteratively  un¬ 
stable.  Such  regions  have  been  indicated  U  (un¬ 
stable)  in  Fig.  14.  If  all  extrema  lie  inside  the 
big  box  at  Fig.  13,  then  statement  (3)  applies  and 
/  is  iteratively  stable.  B  the  central  extremum 
lies  outside  the  big  box  and  the  x,  extrema  lie 
inside  or  vice  versa,  then  there  are  two  possibili¬ 
ties:  B  an  inside  extremum  is  small-bootable, 
then /  is  stable  and  unimodal  [statement  (2)],  if 
not,  then  the  situation  remains  in  doubt.  These 
doubtful  regions  are  denoted  D  in  Fig.  14.  We 
Investigated  them  numerically  by  the  following 
method:  We  fix  (a, 6)  ~nd  follow  the  itinerary  of 
the  “inside"  extremum  for  32  steps.  B  the  itin¬ 
erary  ever  reachea  the  region  |x|>|xm.x(a,6)|I 
then  it  is  unbounded  and/  is  iteratively  unstable. 

B  the  itinerary  remains  within  the  range  - 1  xmv| 

<  x  <  I  Xn»\l  »  then  we  declare  it  bounded  (for  prac¬ 
tical  purposes)  and/  is  iteratively  stable  and 
unimodal.  In  practice,  the  distinction  is  sharp. 
Figure  15  shows  the  results  of  a  scan  In  (a,b)  in 
region  D,  at  a  grid  Interval  of  0.015.  Similar 
structure  occurs  in  the  other  D  regions.  What 
appears  la  a  complex  set  of  spines  or  spiculxn 
radiating  outward  from  the  known-stable  region 
S.  The  major  spicules  are  readily  identifiable 
with  low- it  regions  at  a-cycle  stability,  which  we 
have  sees  in  Figs.  9,  13,  and  13.  Many  other 
cpieulea  certainly  exist  in  these  regions,  however, 
they  become  rapidly  more  narrow  as  a  increases 
and  slip  between  the  mesh  of  our  grid  seen.  Well- 
defined  regions  of  entirely  unstable  behavior  ere 
also  suggested  by  Fig.  15.  The  boundary  between 
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FIG.  IS.  Iterative  stability  In  the  doubtful  region  Di 
of  Fig.  14.  Regions  S  and  t?i  are  guaranteed  stable  and 
unstable,  respectively,  by  the  box  construction.  Points 
shown  were  determined  iteratively  stable  fay  a  200  x  200 
scan,  as  described  In  the  text.  The  major  spicules  are 
readily  identifiable  with  low-w  (a  ■  2,3)  a  cycles,  as 

indicated. 


stable  and  unstable  regions  appears  to  be  very 
Irregular  and  may  have  some  fractal  character. 
All  the  above  considerations  could  be  further  re¬ 
fined  by  studying  boot  constructions  for  the  func¬ 
tions  fM  for  »>1. 

IH.  CRITICAL  AND  TRICRITfOU.  BEHAVIOR 

A  The  Fdeeninum  critical  Use  aad  ita  deal 

The  nesting  1-  and  2-cycle  regions  Shown  in 
Fig.  10  are  the  first  stages  of  a  standard  bifurca¬ 
tion  cascade,  1,2, 4,8, ...  .  When  this  cascade 
proceeds  smoothly  to  completion.  It  terminates 
at  a  locus  of  accumulation  points,  which  we  shall 
call  a  critical  line.4  We  already  know  two  points 
on  the  critical  line:  fo,b)  =  (-1.401 16, 0)  on  the 
negative-*  (logistic)  axle  and  the  intersection  at 
the  critical  line  and  the  self -dual  line  given  at 
the  end  of  Sec.  HB.  There  hi  no  difficulty  In  find¬ 
ing  the  rest  of  the  critical  Use  numerically.  The 
most  efficient  method  is  stmpty  to  plot  the  appro¬ 
priate  loci  of  central-peak  super  stable  4  cycles 
(Fig.  11),  S  cycles,  etc.,  which  rapidly  converge 
to  the  Feigenbanm  critical  line  TT',  a,*  0,(6) 

[or  *,«&»],  Anna  to  Fig.  10.  Approaching  this 
line  from  the  upper  right  always  toads  through  the 
standard  cascade;  continuing  beyond  to  tbs  lower 
left  reveals  the  usual  asquenca  of  inverse-bifurca¬ 
tion  bands.  Asymptotically  close  to  tbs  critical 
line  tbs  position  of  tbs  super  stable  T  cycles  (or 
any  other  reference  feature  one  cares  to  choose) 
varies  [cX  (1.2)]  ss 

s.®)-s,(l)M0)/8; 


or 

6,(s)»6,W+C(s)/6;,  (8.1) 

where  the  critical  amplitudes  A  (6)  or  c( s)  vary 
smoothly  along  the  critical  line  but  5,  la  universal. 

The  new  feature  here  is  the  fact  that  the  critical 
line  terminates  at  the  special  points 

T»(0,  -1. 594 90),  r  -(-2.81403,  1.40701), 

(3.2) 

which  we  shall  refer  to  as  tricritical  points  to 
analogy  to  the  theory  of  thermodynamic  phase 
transitions.19  The  mechanism  tor  the  termina¬ 
tion  la  the  existence  of  new  n-cycle  structures 
which  pinch  off  the  orderly  sequencing  character¬ 
istic  of  the  critical  line.  A  variety  of  interrelated 
structures  exist,  which  we  shall  study  In  mare  de¬ 
tail  in  Sec.  me.  Far  the  moment  we  mention 
three  such  structures:  (1)  A  tricritical  point  is  the 
limit  of  a  sequence  of  doubty  stable  2* -cycle 
points.  For  example,  the  points  r  and  7*  are  the 
limit  of  the  dynamically  coupled  points  listed  in 
Table  I,  the  first  two  of  which  are  visible  In  Figs. 

8  and  11.  (li)  A  tricritical  point  la  the  limit  of  a 
sequence  of  dynamically  independent  7-cycle  in¬ 
tersections  (the  first  few  are  also  visible  in  Figs. 

8  and  11).  Aa  a  consequence  of  (1)  and  (it)  any 
neighborhood  of  a  tricritical  point  contains  an  in¬ 
finite  sequence  of  interlaced  unlmodal  regions  (as¬ 
sociated  with  the  doubly  stable  points)  and  btmodal 
regions  (associated  with  the  dynamically  indepen¬ 
dent  intersections),  (ill)  Finally,  a  tricritical 
point  is  also  characterized  as  a  point  at  which  the 
critical  line  Is  pinched  off  by  a  sequence  of  nar¬ 
row  fingers  (see  the  end  of  Sec.  II C)  of  high-order 
out-of-order  T  cycles.  Such  a  finger  is  already 
visible  In  the  4  cycles  (Fig.  8)  in  quadrant  IV, 
where  it  is  still  well  away  from  T  [the  Anger  tip 
only  reaches  (0.80,  -1.25);  however,  the  corre¬ 
sponding  32-cycle  finger  reaches  (0.03,  -1.82) 
and  higher -order  cycles  close  in  rapidly]. 

Thanks  to  Felgenbaum9  we  know  that  the  critical 


TABUS  I.  Sequences  of  doUbty  stable  f -cycle  points 
(a, 4)  converging  to  the  tricritical  potato  T  and  T1 ,  which 
terminate  the  Fetgenfasum  critical  tine.  The  corre¬ 
sponding  potato  for  t  aad  f  may  bn  obtained  by  using 
the  duality  relatione  (2.13). _ 


Cycle 

T 

r* 

1 

9. 

(-*,  1) 

4 

0,  -1AO* to) 

<-*.00*10,  144*11) 

• 

<*,  -140*26) 

(—*.79*7*.  14*786) 

M 

0.  -1AM 14) 

(—1411*6,  14*67*) 

to 

0,  -14*4*01 

<-*413  to,  14*8*4) 

m 

<*,  -1404*0) 

(-44MW,  lAtftfil) 
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line  is  Associated  with  a  fixed  point  of  the  function¬ 
al  recursion  relation, 

i.e., 

f  m*l  •  (S-3) 

The  factor  /„{ 1)  in  (3.3)  is  chosen  to  preserve 
normalization,  /„,,(())  =  1  provided  /„(0)  =  1.  At 
any  fixed  point  of  (3.3)  (and,  as  we  shall  see, 
there  are  many), 

/* -A,  («)•/*»  oV<«>,  (3.4) 

with 

«=1//*(1),  (3.5) 

i.e.,  f*  is  topologically  conjugate  to/-**  under 
rescaling  by  magnification  factor  a.  The  mapping 
(3.4)  takes  a  function /„=/*+ 8/  near/*  into  an¬ 
other  function =/ * + 8/'  near/*,  with 

8/'  =  £(6/)  +0(8/*),  (3.8) 

where  £  is  linear.  The  eigenvalues  A„  and  eigen¬ 
vectors  em  of  £  satisfy 

£(«.)”  A.e,  (3.7) 

and  provide  a  convenient  basis  for  describing  the 
region  of  function  space  near /*.  For  example, 
if 

m 

far  some  set  of  coefficients  am,  then 

(«•») 

m 

The  fixed  point  (3.4)  of  (3.3)  associated  with  the 
Felgenbanm  critical  line  has  the  form* 

/♦<*)  ■  1  +ci,)x**c</>x*+  ••• ,  (8.10) 

with 

c/’- -1.5275,  c?*- 0.1048,..., 

and 

a*  *2.50381.  (3.11) 

There  is  only  one  “relevant"  *•  eigenvalue  At«8, 
-4.56020.  AU  other  eigenvalues  an  “irrelevant, ” 

|  A„|  <  1  for  m  >  1,  so  the  corresponding  perturba¬ 
tions  (3.8)  shrink  under  tteration.  The  fixed  point 
ff  therefor  •  has  an  bd  lowing  critical  hyparsttrface 
of  csdtonansfflii  one.  The  Mgoubuum  critical  line 
to  the  intersection  of  this  critical  hyper  surface 
with  the  a,  8  pitas.  Any  function  /,  on  the  Fsigen- 
bnuat  critical  Has  (but  efeladHg  the  endpoints  T 
and  f)  flows  to//  under  iteration  (3.3).  [Nearby 


functions  flow  near  to  ff,  through  the  linear  region 
(3.8),  (3.8),  and  then  (3.1)  follows.]  In  particular, 
if  /„(*)■/(*),  then  iteration  under  (3.3)  gives 

/.  •  hiH/^HO))  •/«a*>  •  h4(/(,*>(0)) ,  (8.12) 

with  a 

/<**)(0)  =  q/.(l).  (3.13) 

So,  if  /,  is  on  the  Feigenbaum  critical  line,  then 

Um  (/„).=//.  (3.14) 

Equation  (3.14)  is,  in  fact,  a  convenient  way  of 
finding  the  fixed-point  function//.  In  addition, 
putting  x  *  1  in  (3.14)  and  using  (3.12)  and  the  de¬ 
finition  (3.5)  of  the  magnification  factor,  we  find 

ar  *  lim[/J**>(0)/y[,*M,(0)] ,  (3.15) 

which  is  numerically  convenient. 

The  dual  of  the  Feigenbaum  critical  line  Tf  is 
another  critical  line  ff,  shown  in  Fig.  10.  End¬ 
points  of  this  dual  line  are  tricritlcal  points  lo¬ 
cated  dual  to  (3.2), 

f  =  (-3.18880,  2.54371), 

f  =  (0.855  61,  -1.148  81) .  (,,16) 

The  dual  critical  line  passes  through  (0, 0)  be¬ 
cause  Tf  crosses  (2.15)  and  through  infinity  be¬ 
cause  Tf  crosses  the  a  axis.  The  critical  lines 
TT  und  ff  intersect  at  the  self-dual  line;  how¬ 
ever,  they  are  dynamically  independent:  Feigen¬ 
baum  critical  behavior  uses  the  central  extremum, 
while  its  dud  uses  the  xt  extrema,  in  accordance 
with  the  discussion  after  (2.14).  The  continuity 
of  the  duality  relations  (2.13)  guarantees  that  (3.1) 
holds  in  the  vicinity  of  the  dual  critical  line.  Does 
the  dual  line  lie  on  the  critical  hypersuHhce  of  the 
fixed  point//?  The  answer  is  no.  Indeed,  the 
iteration  procedure  (3.2)~(3.4)  does  not  even  lead 
to  a  function  with  finite  coefficients  like  (3.10). 

The  reason  is  simple  enough.  The  2"  cycles  in 
the  bifurcation  behavior  of  ff  uae  one  of  the  side 
extrema  xt  instead  of  the  central  extremum.  How¬ 
ever,  the  iteration  (3.3)  continually  magnifies  the 
function  around  the  origin,  poshing  the  ectiue  part 
of  the  function  near  the  *,  extrema  off  to  Infinity. 
To  achieve  a  finite  fixed  point,  it  would  be  neces¬ 
sary  to  translate  the  active  peek  to  the  origin  be¬ 
fore  applying  (3.3),  (3.13),  and  (3.14).  The  Infinite 
fixed  point  which  coslroi*  ff  is  topologically  con¬ 
jugate  to// in  da#  ssase. 

»  .  • 

E  tdnMtetey:  Eapcma* 

The  tricrltioal  pohd  T  is  Iseated  ns  the  titett  of 
the  set  at  doubly  stable  8*  cycles  a  tasted  on  the 
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negative -6  axis.  Locating  T*  would  ajipear  to  re¬ 
quire  a  two -dime  nsiooal  search.  Luckily,  this  is 
not  so,  as  we  now  explain.  For  each  n>  I  there 
are  several  dosbiy  stable  points  in  qsadrast  n 
(see  See.  HID).  Bach  such  cycle  uses  the  extrema 
x0  and  x,;  however,  the  different  doubly  stable 
points  use  diem  in  different  orders.  Hie  sequence 
of  doubly  stable  2*  cycles  which  converges  to  T 
turns  out  to  use  the  two  extrema  in  the  order  x„ 

x+, - Normalisation  requires  /(*o)  mfW  “  1 

=x+td,6),  so  a  =  -2b,  as  s— Sap!  If  led  hi  Table  L 
thus,  the  doubly  stable  points  leading  to  T'  all 
lie  on  the  line 

/(x)»l-2 ex*+cx4  (2.17) 

and  can  be  found  by  a  one -dimensional  search. 

We  turn  now  to  trier  ttical  behavior  and  start 
with  an  analysis  of  the  point  T.  Functions  on  the 
6  axis  have  the  form 

/(x)  =  l+6x4  (1.1*) 

and  are  iteratively  stable  aad  unimodalfor  -2  «1 
<  {fg.  The  iterative  behavior  of  such  functions  as 
b  travels  down  the  negative -b  axis  Is  entirety  sim¬ 
ilar  to  that  of  the  functions  (1.2):  There  is  a  amin 
bifurcation  sequence  (the  beginnings  of  which  we 
have  already  seen  in  Figs.  10  and  11  aad  Thbte  0, 
which  has  its  limit  at  »r*  -1.804201 2*02282. . .  . 
Far  -2<b<br  there  is  a  corresponding  sage  esc* 
at  inverse-bifurcation  bands.  The  vteWKy^  of  *r 
is  characterised  by  scaling  behavior.  For  sram- 
ple,  the  positions  *.  of  the  doubly  stable  2"  eyeles 
for  6>br  and  large  n  vary  as 

6.-b,-B/(6y»r ,  (2.19) 

with  >7.284  *9  sad  £•  4.288,  In  preclra  paral¬ 
lelism  to  (1.3).  The  number  differs  from  B, 
because  (3.18)  has  a  quartic  nibsr  than  quadratic 
maximum*  but  is  otherwise  universal  in  the  same 
sense.  The  amplitude  B  is  uownatvereal. 

The  origin  of  the  scaling  behavior  (3.19)  is  the 
existence  at  a  new  fixed  point  of  (3.3),  which  we 
denote  [cf.  (3.10)] 

fffc)*l+el,V+eS,^*‘*%  (3.20) 

where/|  aarttflea  (3.4)  with megXEtr Mien  fast «r 

«,-l//|fl)»  (1  +c\*>+c\n+  •  •  •  f  13.31) 

There  are  rarer al  ways  of  finding  ths  trier itlcal 
fhmd-potnt  fraction.*’  Ora  of  thpsbnptsst  ways 
follows  ths  logic  of  (S.13j-(3.15):  whn  2  •bf, 
/(«)"/rWal+*fs4.  Thao 

n*  U«(/A  '  (3.18) 

a«d  ..*1  Ainr  *:'l  vh»’«T 


«,«llm{/p*»(0>//£"“»C0)l.  (3.33) 

Equation  (3.23)  can  be  used  to  find  the  coefficients 
ci?  as  follows:  Let 

/?’(*) = tfPW+CSV  ♦  •••,  (3.84) 

then  It  to  assy  to  develop  the  recursion  relation, 
Cj~*>  *  4b#t/^>(0)],C</> .  (3.25) 

However,  so 

end  inclusion  of  ths  additional  factors  in  (3.12) 
yields 

»■ 

cin -  Urn  4(46)**  ]j [/^>(0)1*.  (3.28) 

Derivation  at  similar  expressions  for  higher -order 
coefficients  to  straightforward.  Wo  find  from 
(2.23) 

•  r- -1.88020...  .  (3.27) 

ResnRsfor  C(4P  are  given  in  Table  I L 
Perturbation  about  the  fixed  point  (3.20)  deter¬ 
mines  scaling  behavior,  as  outlined  in  8ec.  IDA. 
There  are  two  types  of  perturbations,  6/-  5/“’ 
with 

8 /‘•>-a«x4+«tx*+---  (3.28) 

and 

5f<#)  =s,x,+a,x,+  •  •  * ,  (3.20) 

hnroiviag  powers  x“  wtth  even  or  odd  A,  respec¬ 
tively.  Because  fjktoa  polynomial  tax4,  the 
iteration  of  asuvsa  V  cannot  generate  any  odd 
component.  Humes,  in  terms  of  »($!>),  we 
hsve 


-r  vtr  v 

Vo 


TABLE  0.  Cnafflolsata  of  Um  trtorttteul  flxed-potat 

fcnaMon/  t(d*  l*«l*V44!rV* - These  eoto 

flPtaata  were  atoalaad  Ns  CT.2i|rad  its  htftor-owtar  i 
•tags.  Nate  Hat  toe  asm  ruts  (SJl)  ts  writ  satisfied. 

. -- -  . ---■•r- 

4n»  a  .sues 


SUpUi  im  i  I|||,JB^>  > 
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Although  0/u>  couples  to  both  and 
S/(<>  couples  only  to  tf'M.  The  reduciblllty  of 
(6/)  into  a  aemiblock  form  Implies  that  the  eigen¬ 
value  problems  (3.7)  for  the  even  and  odd  pertur¬ 
bations  do  not  mix  and  may  be  considered  inde¬ 
pendently.  Both  even  and  odd  sub  spaces  have,  of 
course,  an  infinite  number  of  eigenvalues  and 
eigenvectors;  however,  in  each  case  only  one 
eigenvalue  is  relevant  ( |  A )  >1)  and  the  rest  are 
irrelevant  ( |  A  |  <  1).  The  largest  even  eigenvalue 
controls  flows  in  the  8  (*4)  direction  near  tricriti¬ 
cality  and  is  Just  the  universal  number  intro¬ 
duced  in  (3.19)  and  determined  numerically  from 


superstable  S’-cycle  spacings.  Finding  it  and  the 
other  even  eigenvalues  by  solving  (3.7)  with  (3.38) 
is  feasible  bat  uninteresting.  The  largest  odd 
eigenvalues  we  denote  We  shall  show  below 
that  is  simply  related  to  the  magnification 
factor 

a<»»a*«2.88713,  (3.30) 

and  that  the  irrelevant  odd  eigenvalues  are  Just 
the  Inverse  odd  powers  of  (3.30),  (a^)"****1*.  To 
see  this,  take /„=/)+ 8/  and  iterate  (3.3)  once  to 
get 


/•«(*)  =/ fW+e/'W 


The  first  and  third  terms  in  the  brackets  are  al¬ 
ways  even;  only  the  second  tram  can  be  odd.  As 
we  mentioned  earlier,  we  only  need  to  work  In 
the  odd  k  apace  to  obtain  the  desired  eigenvalue. 
Thus,  to  linear  order,  we  have 

0fM'  <*)=«',**  +«',*•+*•* 

<»■«> 

Equation  (3.32)  has  the  structure 

< 'E^mUc^Ds.,  *,1-2,8,10 . (8.83) 

f 

The  eigeurfttae*  of  M  are  the  odd  eigftmlue*  of 
(3.7).  Because  the  coefficient  of  e/(*>  in  (I.  S3)  la 
a  polynomial  In  x\  o',  cannot  depend  on  a,  for 
i>»,  so  Jf  is  triangular  and  its  eigenvalues  are 
simply  lfu.  To  find  if*,  we  can  set  x  •  0  in  the 
coefficient  of  8 fM  to  get 

(S.MI 

However, ff  satisfies  (3.4), 


lotting  x—  0  lends  to 

wr<**/#wh 

so  finally, 

which  glvee  (3.S0)  sad  Hi 


dimension  equal  to  2.  This  hypersurface  evident¬ 
ly  Intersects  the  a,b  plane  at  T  and  controls  the 
scaling  behavior  there  (see  Sec.  ®C).  ft  is  hy  no 
means  obvious  that  the  trier ttical  point  T  at  the 
other  end  of  the  Felgenbaum  line  la  also  on  the 
trlcritical  hyper  surface  belonging  to/js.  We  have 
verified  numerically  that  it  Is  by  observing  that 
lim.^.0>Wf  [via  (3.23),  (3.28),  etc.].  It  fol¬ 
lows  that  behavior -near  T  i a  described  by  the 
eigenvalues  8$P,  and  a,.. 

We  tn#n  now  to  the  trtcrttlcal  points  terminating 
the  dual  critical  line  ff.  Duality  guarantees  that 
scaling  behavior  near  f  and  f  is  described  by  the 
universal  numbers  and  tip.  B  doss  not  follow, 
however,  that  f  and  f  flow  to  the  fixed  point  ff 
and.  Indeed,  they  do  not  (cf.  the  situation  with  the 
daal  critical  line).  Direct  iteration  In  analogy  to 
(3.22M3.28)  Mrows  that  #  and  f  flow  to  Mother 
trlerttfeal  fixed  point,  given  by 

(3.88) 


af-aj.  (3.29) 

Mote  that  (2.28)  is  Just  of  the  form  (2.4),  with  *(x) 
■  /*,  *'1(x)-x’,  and  ao/f  is  topologically  conju¬ 
gate  toff.  IMnse  amusing  nbitlSns  era  not  hard 
to  derive:  When  applied  to  points  an  the  ft  axU, 
fibs  duality  relab  win  (1.10)  and  (2.11)  read 

«!<»)«  1+Wx*,  «,&)-**,  (8.40) 

so 

fiM "/He) •Qjfidk&m  &  *  . 

(8,41) 

That,  is  aaafegy  with  (2.22) 


NfcnXjOrr'  - 
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Similarly, 


(3.42) 


/J«Um  (/**,,  (3.43) 

r** 

and  substitution  of  (3.41)  toads  to  (3.38). 

C.  TilciUinly :  CtMNWfldsndutoctioM 

Hie  trlcritical  fixed-point  discussion  of  the  pre¬ 
ceding  section  sets  the  stage  lor  an  examination  of 
the  structure  of  «  cycles,  etc.  (the  phase  diagram) 
near  the  points  T,  T,  and  their  duals.  B  ias 
property  at  the  iteration  (3.3)  that,  M  /,  has  an 
m  cycle,  xuxv. ..  ,x„  with  even  m,  then/,,, 

has  two  w/3  cycles,  *,//,(!),  . and 

x*/f,Q),xJfmQ), - V  the  as  cycle  of/,  is  sa- 

per stable,  then  both  file  m/t  cycles  of  /**  are 
super  stable  (/,M  is  at  tonst  Mmodal).  If  si  1* 
odd,  the/.,,  has  an  m  cycle,  *,//.(D,x,/y;<i), . . . , 
xJfJX),  xJJJXl, ....  Again,  /,  auger  stable  itn- 
plies/^  super toable. 

to  the  linear  rsgtos  (3. «)-{«.»)  near  a  fixed 
poll*,  theee  considerations  take  a  particularly 
concrete  farm.  Consider,  for  example,  a  tone  - 
ttdh /,  #tth  a  ryele  which  lies  naar/f  MB  hot 
on  the  tricrttical  hyper  earthen.  3ueh  a  twtttidn 
has  the  farm 

/.«/#♦  V. 

with 

8/W  «Stoi  W  +«-e.<*) ,  (3.44) 

where  •,(*)  end  «,(x)  are  the  eigenvectors belong- 
tog  to  sy*  and  8^*,  respectively.  The  iterate  func¬ 
tion  (which  has  two  3“‘‘  cycles)  la 

/*»*/#♦*% 

Wit#  (3.43) 

V*  W-3^M*W+»VW>)  i 
hs./frlM  etefttosway  frees  the  trtoritlaal  kypct- 

«V*  to  toe  sMMmo  [see  (3.33)  aed  (>.»•)]. 


terveuing  iterations.  These  scaling  variabUs  are 
a  natural  sot  of  coordinates  near  the  point  T.  They 
may  1m  taken  to  be 

ut=b -6r-atan0,  u, 1  a/cosO ,  (3.48) 

where  the  slope  at  the  critical  line  is 

tend- -1.2347,  (3.47) 

as  Illustrated  in  Fig.  18.  Cycle  structure  near  T 
has  much  regularity,  when  described  in  terms  of 
%  and  u~  For  example,  if  there  is  a  super  stable 
2*  cycle  at  then  there  Is  a  super  stable 

2**1  cycle  at  (8V^,8^n^,  a  auperstebie  V*  cy¬ 
cle  at  {(•9>)^ii,  (iV® )*“,).  etc.  peavfasd  only  that 
all  coordinates  remain  mull  enough  so  the  potato 
have  not  deported  stgeMtcnaSly  from  T  OLe.,  from 
the  trirrWcal  hypersarface). 

Thews  ccastthnttme  allow  a  scaling**  descrip¬ 
tion  ot  the  eroaacrer  from  critical  behavior  along 
the  Feigenbmtm  Ban  to  tricrttical  behavior  along 
the  h  axte  (region  2  of  the  scaling  variables  to 
Fig.  Id).  The  position  of  the  saperatahie 

2*  cycise  in  this  region  has  a  simple  scaling  form, 

ItpW'-HWru},  (3-48) 

which  in  valid  with  a  tingle  function  9  independent 
at  n,  presided  **>1  and  |«,  | ,  |«,|  eel  bid  for  ar¬ 
bitrary  vetoes  ottbe  scaled  comhinattona  (S^m. 
Equation  (3.48)  mast,  of  course,  incorporate  the 
known  behavior  on  the  b  axis  aad  along  the  critical 
line.  Ontheh  agls  a^O,  ao  (3.19) 
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N(0) -a.  (*.49) 

On  th«  other  hand,  at  fixed  «t#0  (3.48)  muet  agree 
with  (3.1)  aa  0,  so  at  large  %  S(x)~xy  and, 
comparing 

-C/8J ,  (8.60) 

we  deduce  that  the  critical  amplitude  C(a)  vanishes 
as  e—  0  as  s',  with 


which  we  have  verified  numerically. 

Such  seating  considerations  are  not  restricted  to 
region  3  (Fig.  18).  Indeed,  the  a -cycle  structure 
ia  the  entire  oeighborbood  of  T  scales  in  accor¬ 
dance  with  the  discussion  after  (3.47).  This  point 
is  iftastmted  by  Tigs.  17-30,  which  show  the 
structure  of  3*  cycles  for  **3,  3,  and  4  displayed 
in  terms  at  the  scaling*  combinations  (»“’)"%  and 
(Of'Tu,.  The  structure  displayed  should  become 
fixed  ess—  •».  For  a = 2  some  oooscallng  fea¬ 
tures  awe  observable  far  from  the  origin;  however, 
a  =  3  and  a=4  are  already  practically  indistin¬ 
guishable.  Note  that  at  fixed  djs,  large  k  probes 
email  s;  thus.  Pigs.  17-30  illustrate  graphically 
the  high-order  3*-cyele  structures  which  pinch  off 
the  critical  line  at  T.  A  similar  discussion  at  be- 
Nwur  near  T’,  f,  and  f"  could  be  made. 


In  this  peper  we  have  Confined  our  discussion 
to  the  main  Feigenbaum  critical  line  TV  and  its 
dual.  We  note  in  closing  that  tricrltical  struc¬ 
tures  exist  at  many  Other  points  in  the  a,b  plane. 


HO.  18.  Superstate  8  cycles  near  T  plotted  with 
scaling  variables  KtAr)>-  The  point  marked  by  the  arrow 
la  doubly  stride,  lids  figure  and  Fig.  18  oofneide  ex¬ 
cept  for  amail  uatnerloal  dMfomsoes. 


Let  ue  focus  on  the  characterisation  (Sec.  IIIA) 
of  a  tricrltical  point  as  the  it—  •»  Umit  of  a  set  of 
doubly  stable  2* -cycle  points .  Consider  only 
quadrant  IF.  The  3  cycles  have  a  single  doubly 
stable  point  (Fig.  8)  at  (0,-1).  The  4  cycles  have 
(Fig.  11)  two  suchpoints,  (0,-1.80893)  and 
(0.81179,-1.339*3).  The  B  cycles  have  four  such 
points,  etc.  This  sequence  IS  llustrated  in  Fig. 
31  through  order  a >3.  R  appears  that,  in  going 
from  *  to  a  + 1,  each  doubly  stable  point  splits, 
forming  as  »— •  a  whole  cascade  of  tricrltical 
points.  The  outer  sequence  leads  to  Tend  f\ 
Others  Ue  between  T  sad  f’ .  Bach  limiting  point 
has  an  infinite  number  of  others  in  its  nelgbboi  - 


afftt 
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fb.  20.  apwMbli  a  cyolea  fooutad  wMh  Ri  ex¬ 
treme  x,  vd  x.  (x.  fcaa  been  omitted  (or  eUnpUclty). 
Otherwise,  this  is  the  seam  m  Fig.  IS  esmpt  tor  Mm 
larger  Mala.  Note  the  addtttoual  doubly  states  print*. 


hood.  Figures  17-lS  show  two  Rats  at  doubly 
•taMo  J"-eycle  points  converging  to  t,  <ma  on  the 
«%=*0  (A)  axis  and  tbs  other  located  at 

*>  =  (66.4,10.2).  Figure  30,  which  la  on  a 
larger  scats,  shows  two  additional  doubly  stable 
points,  aad  it  seeans  likely  that  there  are  many 
more  at  larger  mines  ef  the  sealing  variables. 

Tbs  satf-slmllarity  observable  ta  the  sealing  plots 
suggest  that  tbs  trlerltieal  points  may  form  a  Can¬ 
tor  set.  At  es-ih  level  of  a,  we  may  label  a  Y* 
dmMy  stabts  cycle  msetioesd  above  as  right  (A)  or 
left  (A)  relative  to  its  parent  Wa  can  teen  specify 
a  trlerltieal  point  by  a  hiwwii  of  if  s  aad  fa 
Conversely,  any  such  a  R-L  sequence  specifies 
a  tricrttical  point.  For  tnatasce,  the  tricrttical 
point  T  is  R*  ■  lR,R, ...)  and  the  trlerltieal  point 
fi»t*.  In  Fig.  21,  wo  havo  also  plotted  the  lo¬ 
cations  of  two  additional  trlerltieal  points  T, 
•(0.M0I1.-1.7M0) aad  T,-(0. 9704,-1. *7*7). 
la  terma  of  the  R-L  notation,  T,  is  RV  and  T,  is 
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FW.  21.  Doubly  stable  *■  ovule*  in  the  region  between 
Taadf'.  Thn  paints  T  sad  f*  an  marted  by  crosses, 
la  this  regies  there  are  con  2  oyeis,  two  4  cycles,  four 
•  eyalos,  ate.  AU  the  s  -•  limit  patois  are  presumably 
trtoritioal. 


LA*.  We  have  not  at  this  stage  located  Ue  criti¬ 
cal  tinea  associated  with  these  tricrttical  points, 
nor  have  we  stptored  the  question  of  what  fixed 
point  they  flow  to  under  the  iteration  (3.S).  Pre¬ 
liminary  evidence  shows  other  doubly  stable  points 
top,  fba  •  eyeleo  and  above,  which  may  veil  belong 
to  additional  trlerltieal  sequences . 
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is  equlvaieat  toa-J,  nkda  X«4  ftayrnd  whiah  aU of 
tea  mtft  talarvata  awapt  a  sat  of  mauur*  aero  maps 
te|s)-*)  te  — dwta*  tee— k. 
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S.  J.  Chang.  it  Wortts,  and  i.  Wrigh t,  Xrioritioal 
potato  nd  Uforeattana  in  aquartto  map  D  pnm). 

"For  oar  purposes  Van  to  no  need  to  dtoBnpUta 
between  too  ergod to  and  nonperiodic  itineraries  of  Ref. 
t. 

‘*Cae  paa lr  laapa  with  this  property  are  called  "moat 
stable”  to  Bat.  U.  to  tola  paper,  we  follow  toe  conven¬ 
tion  «(  Bala.  V  and  S  and  refer  Van  as  super  a  table.  A 
cycle  la  called  doubly  stable  if  it  is  doubly  super  stable. 

“The  terra  “ant mortal"  la  used  to  a  similar  but  more 
reetrlotod  sense  to  Hat.  I. 

nThe  converse  la  proved  as  follows:  If  the  s  cycle  of  / 
uses  S|«0,  than  there  exists  a  y,  such  that  x, 
-GtfGjlyt)).  Application  of  (2.14)  to  y,  gives  a  y « which 
Is  a  side  extremum  of/. 

“A  ftmntlcn  /(»)  to  Iks  1-Uyole  region  mar  be  unimcdal 
or  bbnortal.  A a  we  shall  sea  below.  It  to  quite  poa- 
etoto  tee  a  stablas  apele  tonnaalat  wlto  a  stable  e» 
cycle  (*  »■*  or  a  ■*  m  )  or  with  ehsotto  trajectories. 

‘•There  are  nuegr  otoer  Ralgtahiiim  critical  Uses  to  toe 
phase  fRagnon.  The  “main"  orffteaft  Haas  shown  to 
Kg-  Id  have  a  complete  btfuroatioa  eeqnenoe  X,  2,4,8, 
...  .  The  other  critical  Uses  have  truncated  sequences 
2,4,8,  ... ,  4,8,16,  ... ,  and  so  forth, 

"Note,  for  example,  that  toe  lower  of  the  two  main 
diagonal  t  aynla  Haas  to  Kg.  11  crosses  the  negative 
e  ante  bopoade,.  It  corresponds  to  the  “ehnolto  4 
oycle. "  whlob  oocurs  in  tip)  region  of  Inverse  hUaroa- 
ttee. 

**K  the  quadratic  map  (I  A)  to  a  good  guide,  then  we  may 
expeot  that  thea-oyete  eaglhae  are  danse  (Baf.  2), 
eves  1*iriig*-  toe  ntoultta  reghtoe  have  Suite  measure 
fronjeotared  in  Baf.  t<  proven  by  it.  JatobeOn  tnipifi- 
lished)). 

°D.  Stager,  SUM  W.  Appl-  Mrtll.)  Bsv.  38,  240  (1878). 

"Whan  (2.2S)  la  satisfied  to  quadrant  HI,  She  whole 
aa  awra  eMhe  Iterative  domain  at/  to  apjataialy 
eahanaSad  bp  aHhta  a  atogla  atahle  a  spate  a*  a  atogla 


region  of  ohaotic  ltlnarariea.  In  principle,  there  are 
two  adrttttoaal  poealbfllUaa!  statultaneoue  exlatenoa 
of  two  or  more  duttaot  ehaotto  regions  or  ooextotance 
of  one  etabla  a  cynle  wlto  one  or  more  dtottaot  chaotic 
raghias  These  posathmttei  do  not  seem  to  ooour  to 
pipfttioe;  however,  to  the  beat  of  our  knowledge  they 
base  not  beta  rigorously  ruled  out,  even  for  toe 
sup. 

"Agate,  tears  exists  toe  logtoal  poaatbUUp  that/  pos¬ 
sesses  some  aonaero  region  of  ehaotto  itineraries.  So 
tar  as  we  know  tote  does  not  ooour;  however,  wa  are 
aware  of  no  rigorous  argument  to  enolnda  It. 

“to  the  toeoxy  of  pheae  trnrittone  a  trtorittoal  point 
maria  the  point  where  a  Una  of  saeoad-order  transi- 
t'an  |a  ortttoal  Una)  haconaaa  Brat  order.  B  is  char¬ 
acterized  bp  the  appaaranoa  of  a  new  fined  point  with 
aa  extra  relevant  variable.  See  R.  B.  Griffiths,  Kirs. 
Rev.  Lett.  34,  716  41870};  E.  K.  Btodai,  Mi.  28,  676 
'.1872);  K.  K.  Medal  and  I.  J.  Wagner,  Mi.»,  848 
(1972).  lie*  the  reader  be  warned,  however,  tost  the 
parallel  la  imperteot:  tatot  happens  In  toe  bifareaticn 
problem  beyond  Ihe  trlerittoal  paint  to  act  analogous 
to  first-order  behavior,  furthermore,  there  la  not 
hnaa  pus  there  would  be  la  the  phase  transition  prob¬ 
lem)  an  enlarged  parameter  apace  containing  three 
critical  lines  tad  a  triple  line  meeting  at  triorlttoaltty. 

“Wa  follow  hb*e  toe  pbaae-tranaltion  terminology  of, 
for  example.  Bit.  4.  ' 

"one  method  (Hof.  *)  la  to  substitute  (3.20)  Into  the 
fixed  potot  equation  DA)  end  time  tnmneto  to  some 
given  order.  A  tanpteJPWo-ooeSloieattrncexloa  . 
givqa  ci*1  —IM*,  *!* - (U78,  aad«$>*8.8l,  in 
HNMIlMfelt  IflMBUBl  rtnUl  (3JL9) 

"^eerty  example  of  tfoa  meteotf  Was  tee  analyate  Ip 

'  J.  T.BC  add  J.  O.  Utator,  Kpa.  Bav.  I d*L  Jg;  80S 
(MR,  of  tee  raagaetie  equation  of  state  of  CrBr,. 
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Abstract 


Highly  daaped  dyaaalcal  l)WMi  ait ft  lea  dinsnalooal  strange  attractors 
bass  tsaporal  fsaturss  which  cast  bs  described  by  ons  diaaasional  asps.  He 
dsaonstrate  that  tbsrs  should  aiPoheSpstial  fastness  dsscrlbable  by  ons 
dlasaslonal  asps  and  suggest  experlsentaliats  Investigate  certain  spatial 
csrrslatlons. 
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I 

In  this  paper  we  di8cu8s  soae  of  the  phenomenology  of  order  and  chaos  In 
highly  daaped  dynamical  systems,  with  emphasis  on  the  kinds  of  information 
that  can  be  obtained  experimentally.  Our  main  concern  is  with  systems  having 
a  low  dimensional  strange  attractor.  The  particular  examples  chosen 
correspond  to  an  attractor  with  fractal  dimension  very  close  to  two.  As  is 
well  known  from  the  work  of  Lorens  such  systems  have  important  temporal 
features  that  can  be  described  by  one  dimensional  maps.  We  wish  to  point  out 
that  there  will  also  be  spatial  correlations  describable  by  one  dimensional 
maps.  It  does  not  seem  to  be  widely  appreciated  how  much  information  can  be 
obtained  from  such  maps,  and  one  of  the  purposes  of  this  paper  is  to 
illustrate  the  kinds  of  information  that  can  be  obtained  and  to  suggest  that 
certain  experiments  could  be  performed  to  test  some  of  the  Ideas  presented 
here. 

We  will  use  the  Lorens  model ^  to  illustrate  our  points,  but  we  have  also 
done  computer  experiments  on  a  model  with  5  modes  to  verify  that  3  modes  was 
not  special.  Our  main  point  will  concern  spatial  correlations  in  chaotic 
behavior  implied  by  the  existence  of  low  dimensional  attractors.  To 
demonstrate  that  it  will  be  necessary  first  to  discuss  temporal  correlation* 
implied  by  the  existence  of  one  dimensional  maps. 

The  Lorens  equations  are  given  by: 

X  -  Y  -  OCX 

i  m  -  xz  +  x-eY 
i  "XT  -  e  82 


where  we  here  reecaled  the  variables  ee  discussed  in  ref.  2.  The  figures  ere 
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for  e  •  .02216,  o  -  10.,  and  B  -  8/3  which  correspond  to  r  *  l/(e2o)  -  203.64 
in  Loren*' notation*  The  particular  parameters  used  correspond  to  the 
chaotic  region  following  the  first  period  doubling  regies  as  e  is  increased 
frae  sero;  however  very  sieller  behavior  is  obtained  for  the  par centers 
studied  by  Lorens. 

The  tine  series  for  the  variable  X  is  shown  in  Figure  1.  Pictures  of  the 
three  dieenslonal  orbits  (see  refs.  1,  3-4)  show  that  the  notion  is  nearly 
confined  to  a  two  dieenslonal  surface  —  although  it  is  infinite  sheeted.  If 
the  attractor  were  in  fact  two  dieenslonal  we  could  predict  everything  about 
the  systee  froa  a  knowledge  of  any  two  variables  which  designate  the  position 
on  this  two-dleensional  surface.  One  choice  of  these  two  variables  could  he  X 
and  i.  He  construct  our  one  dieenslonal  asps  by  considering  various 
quantities  at  the  tine  when  X  has  an  extreeua,  (i.e.  i  ■  0),  although  any 
other  reference  point  would  suffice.  (8oae  choices  of  reference  points  give 
slapler  naps  than  others.)  Bsaaples  are  given  in  Figures  2  through  4.  Figure 
2  is  slailar  to  "Loren* '  a  Figure  4  and  shows  the  predictability  of  the  value 
of  X  at  its  next  naxlnua  froa  the  value  X  ■  at  the  present  naxlaua.  See 

Loren* 's1^  Figure  3  for  another  exaaple.  Note  that  there  is  a  nonuniqueneas 
that  we  will  discuss  later.  Figure  3(a)  shows  the  predictability  of  the  tine 
to  the  next  esxleue  of  X  froa  the  present  aaxlaua.  Figure  3(b)  shows  the 
predictability  of  the  value  of  Z  at  the  next  aaxlaua  of  X.  Figure  3(c)  shows 
the  veins  of  Z  at  the  tine  of  the  present  aaxlaua  of  X.  In  all  of  the 
exaaples,  the  present  value  of  X  at  its  aaxlaua  X  -  X^  serves  as  a 
predictor. 

Thus  we  expect  that  if  the  data  froa  a  single  probe  experiasnt  gives  a 
eae  dlasnsional  aap  such  as  that  shown  in  figure  2,  then  a  two  probe 
expwrlasnt  should  provide  a  aap  Such  as  that  shown  in  figure  3(c)  (or  figure 


3(b)).  To  dace  no  experlaentalietB  have  looked  at  euch  nape  and  we  auggaat 
that  they  do  ao. 

One  wight  aak  if  it  ia  possible  to  predict  at  tlaea  other  than  extrema 
and  the  anawer  la  yea.  Figure  4  a hows  the  prediction  of  Z  at  a  tine  At  -  20 
after  the  aaximua  of  X.  Other  values  of  At  give  elallar  curves.  It  la  of 
course  not  possible  to  take  At  arbitrarily  large.  There  are  two  reasons  for 
this,  one  of  which  Is  that  these  one  dimensional  naps  display  what  Is  called 
sensitive  dependence  on  initial  conditions  so  that  slight  differences  in  Imat 
are  greatly  magnified  at  large  At  giving  the  wap  a  complicated  structure  which 
Is  difficult  to  resolve.  The  second  Is  that  the  wap  Is  really  not  one 
dimensional,  but  the  lines  all  are  fractal  structures,  or  from  the  point  of 
view  of  an  experimentalist  they  are  thick  lines.  This  thickness  Introduces 
slight  errors  which  are  then  magnified  by  the  sensitive  dependence  on  Initial 
conditions. 

He  now  discuss  briefly  the  nonuniqueness.  Typically  there  will  be  a 
multivalued  problem  that  Is  easily  resolved  by  a  simple  rule.  In  Figure  2  we 
have  divided  the  curve  Into  four  regions  separated  by  *.  If  we  are  In  the 
region  where  we  have  to  select  either  the  curve  C  or  D  wi  do  io  by  the 
following  rule.  Look  at  the  previous  point.  If  It  is  was  on  D  or  A  use  C.  If 
It  was  on  C  or  B  use  D.  The  rule  Is  easily  obtained  by  watching  the 
generation  of  a  sequence  of  a  few  points.  He  have  looked  at  a  number  of 
systems  and  the  appropriate  rule  Is  usually  easy  to  construct. 

Bare  we  have  Ignored  the  thickness  of  the  lines.  If  we  resolve  them,  at 
the  next  level  each  segment  becomes  two  and  our  rule  becomes  more  complicated, 
but  still  tractable.  This  labelling  can  be  continued  to  any  depth  desired, 
but  la  practice  an  experimentalist  can  only  distinguish  a  faw  curves.  Mote 
that  the  previous  wine  is  not  needed,  only  an  integer  label.  This 
serves  to  label  the  other  curves  such  as  thoea  shown  la  Figures  3-4. 


label 


There  heve  been  several  attempts  to  dlacuaa  the  statistics*-^  of  systaas 
of  equations  such  as  the  Lorens  equations.  Although  we  have  nothing  new  to 
add  here.  It  should  be  pointed  out  that  the  amps  we  have  shown  inply  that 
there  are  strong  teaporal  correlations  of  soae  sort  In  such  systeas. 

Typically  the  autocorrelation  function  falls  off  aich  faster  In  time  than  the 
aeaory  of  the  asps.  Standard  statistical  descriptions  don't  correctly  treat 
this  feature.  Soae  different  type  of  description  is  needed  which  takes  Into 
account  the  existence  of  these  asps. 

Finally  we  turn  to  the  aaln  point  of  this  paper  —  the  possibility  of 
spatial  order  in  chaos.  Suppose  one  has  a  eystea  described  by  partial 
differential  aquations.  If  the  systea  Is  expanded  in  teres  of  soae  coaplete 
set  of  aodea,  and  the  nuaber  of  nodes  truncated,  ordinary  differential 
aquations  are  obtained.  Mow  if,  as  we  saw  above,  all  of  the  node  aaplltudes 
are  In  fact  given  in  terns  of  two  nodes  (plus  a  label)  at  any  tiae,  then  it 
follows  that  the  original  field  variables  are  also  given  in  terns  of  two  node 
aaplltudes.  In  fact  the  value  of  one  field,  say  teaperature,  at  one  point  at 
a  tiae  when  this  field  Is  an  extreaua  should  serve  as  a  predictor  for  field 
variables  at  the  sane  or  different  spstlal  points.  To  Illustrate  the  idea  let 
us  return  to  the  Lorens  equations.  The  teaperature  field  6(r,t)  Is  given  in 
terns  of  T  and  Z  and  the  streaa  function  u(r,t)  In  terns  of  X. 

(Kr,t)  -  Kj(/2  Y(t)  sin  vxgcos  kXj  +  ,Z(t)  sin  2«Xj) 

u(r,t)  ■  KgXCt)  sin  *Xj  sin  kx^ 

Where,  Kj,  Xj,  end  k  are  constants  sad  t  -  (Xj.x^Xj). 


In  Figure  5  we  show  the  uee  of  the  temperature  at  one  apatlal  point  at  the 
tiae  of  a  temperature  extremum  at  that  point  as  a  predictor  for  the 
temperature  at  a  different  spatial  point.  The  points  were  chosen  (at  random) 
so  that  the  predictor  is  .134  Y(t)  +  .6Z(t)  and  the  predicted  temperature 
is  82*  .5Z(t)  -  . 146Y(t).  The  temperature  at  its  extremum  at  one  spatial 
point  also  serves  as  a  predictor  for  the  velocity  field  at  any  other  point. 

We  suggest  that  it  should  be  possible  to  test  these  ideas  by  performing 
experiments  —  both  real  and  on  the  computer.  Of  course,  our  examples  have 
been  chosen  to  show  our  Ideas  in  the  most  favorable  light.  Real  systems  may 
have  a  more  complicated  sheet  structure  that  is  hard  to  untangle  —  but  if  one 
is  convinced  there  is  a  rule  and  if  the  maps  appear  to  be  one  dimensional  it 
should  be  possible. 

Specifically,  choose  any  two  quantities,  U(t)  and  V(t).  These  could  be, 
for  exasple,  the  values  of  some  field  at  two  spatially  distinct  points  at  the 
same  time.  Plot  U(t)  vs  V(t)  for  those  values  of  time  for  which  5(t)  ■  0.  If 
this  plot  shows  lines  as  in  Figures  2-5  the  dimension  of  the  attractor  is  near 
two.  If  the  plot  is  scattered  a  higher  dimension  is  indicated. 

In  systems  with  attractors  near  three  dimensions  it  should  be  possible  to 
hake  contour  plots  of  the  value  of  some  variable  at  00a  point  as  a  function  of 
the  values  at  two  other  points  (at  the  time  when  one  of  them  is  an 
extremum).  Note  hare  ws  mead  the  values  of  three  quantities  to  predict 
others.  The  plots  would  present  information  In  the  form  of  Lorens' s1^ 

Figure  3  except  that  there  is  the  additional  information  that  mother  variable 
has  prescribed  value  (i.e.  ft  ■  0). 

Given  a  map  such  as  figure  2,  we  can  estimate  something  about  how  large 
At  (see  figure  4)  can  be  before  temporal  information  is  lost.  At  present  we 
do  not  knew  hew  to  asks  the  corresponding  estimates  for  spatial  structure. 
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Experimental  loforeatlon  where  one  probe  Is  held  fixed  end  e  second  probe  is 
moved  to  obtain  date  at  several  different  spatial  points  mould  be  useful. 

At  this  point  we  would  like  to  mention  the  possibility  of  using  such  maps 
for  design  purposes.  By  changing  parameters  slightly  it  is  possible  to  see 
how  the  map  is  changing,  sad  by  using  a  graphical  Newton's  method 
experimentalists  could  predict  the  parameters  that  give  a  particular  desired 
behavior,  for  exsapla  a  periodic  orbit.  Us  have  used  this  scheme  successfully 
on  the  Lorens  equations  to  locate  periodic  orbits  with  particular  syanatry 
properties. 

lie  conclude  that  in  systems  with  low  dimensional  attractors  there  should 
be  highly  correlated  spatial  structures.  This  could  be  demonstrated  sad 
Investigated  by  asking  asps  (or  contour  plots)  similar  to  those  already  used 
for  displaying  temporal  correlations. 

Supported  in  pert  by  MR,  N00014-80-C-0840. 
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Figure  Ciftltu 

1  Stpnt  of  tte  tine  aeries  for  tho  variable  X  for  tho  Lorioi  equations  for 
tte  ptmitui  c  "  .02216*  o  ■  10.,  I  ■  8/3. 

2  Ntf  stowing  tte  value  of  X  «t  eta  none  seringa  of  X  (MX)  given  tta  vnlnn 
of  X  at  tte  prsNst  tolnn  of  X.  Sognontn  A,B,C,D  nrn  discussed  in  tta 
tost. 

3.  taps  showing:  «)  tta  tint  to  tta  oast  wartnue  of  X  (MIX)  b)  tta  value  of  Z 
at  MMX  and  e)  tte  wine  of  Z  at  tta  present  aarlaaa  of  X(PMX)  as  a 
function  of  X  at  fKX. 

4.  Map  stowing  Z  at  a  tins  At  «  20  after  FKX  as  a  function  of  X  at  MX. 

5.  Map  stowing  tta  tanparatura ,  62 ,  at  one  spatial  point  as  a  function  of  tta 
teaperature,  8j,  at  another  point,  given  that  8}  la  naxlnal. 
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